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Section 1 


REVIEW OF MAJOR ATC PROBLEMS 


The Continental Land Mass (COLM) Air Traffic Control (ATC) Study is essential- 
ly a vehicle for comparing various satellite systems and techniques relative to provid- 
ing CONUS-wide ATC services. In order to provide a meaningful framework for the 
study, a review of the ATC operational environment has been made to highlight the 
significant problem areas. Much of the material presented in Subsections 1. 1 and 1.2 
has been abstracted from an EBM ERAD Study, Reference 1. Source material includes 
various Federal Aviation Administration (FAA) and National Transportation Safety 
Board (NTSB) documents. References 2 through 9. 

1. 1 MID-AIR AND NEAR MID-AIR COLLISIONS 


The mid-air collision (MAC) record of US civil aviation over the past 15 years is 
shown in Table 1-1. The major item of concern, the mid-air collision between two air 
carrier (AC) aircraft, occurred only three times during this period. Two of these 
accidents involved pilot error and the third an anomaly of the control system. The bulk 
of the accidents involving air carriers are caused by collision with a general aviation 
(GA) aircraft. The statistics show an increasing trend during 1966 through 1969 in the 
number of air carrier/general aviation collisions. As expected, the majority of the 
mid-air collisions, 83%, are between two general aviation aircraft operating in a VFR 
environment. 

How do actual accidents compare with the number of accidents expected on the 
basis of aircraft populations? To facilitate this comparisions, data on flight hours and 
operations are summarized in Table 1-2. Since there are too few collisions in any given 
year from which to draw valid statistical conclusions, the total accidents occurring 
during the 15-year period are compared against the aircraft population ratios of 1968 
and the results shown in Figure 1-1. Collisions in the traffic pattern or during final 
approach to an uncontrolled airport have been removed because they are a special prob- 
lem related to airport traffic. All three collision categories involving air carriers are 
lower than predicted since the air carriers operate almost exclusively under instrument 
flight rules. General aviation which operates under visual rules approximately 88% of 
the time is higher than predicted. 
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Table 1-1 

MID-AIR COLLISION RECORD, US CIVIL AVIATION 1 




























Table 1-2 


AIR TRAFFIC ACTIVITY, 1968 


Flight Hours by Flight Rule (millions) 




IFR 

VFR 

Total 


Air Carrier 

7.2 

— 

7.2 


General Aviation 

2.9 

20.4 

23.3 


Military 

4.9 

2.9 

7.8 


Total 

15.0 

23.3 

38.3 


Flight Hours by Aircraft Class (millions) 




Single Engine 

Multi-Engine 

Prop 

Turbo 


Piston 

Piston 

Jet 

Jet 

Air Carrier 

0.3 

0.1 

1.0 

5.8 

General Aviation 

19.0 

3.6 

0.4 

.3 

Military 

0.9 

2.2 

0.5 

4.2 

Total 

20.2 

5.9 

1.9 

10.3 

Operations (millions) 




At airports without FAA towers 




57.7 

At airports with FAA towers — 

VFR 

14.6 



At airports with FAA towers — 

IFR 

40.5 




Subtotal 



55.1 



Total 


112.8 
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Figure 1-1. Actual vs Predicted Mid-Air Collisions by User 
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In January 1968, a special FAA group undertook a study of near mid-air collisions. 
The near miss situation is at best a subjective judgment. Aircraft may pass in close 
proximity to each other, particularly in a terminal area, without the incident being con 
sidered a hazardous condition. A near miss is defined as a situation in which one air- 
craft passed within close proximity, 500 ft or less, and created an unsafe condition. 
Near miss incidents were further divided into hazardous and nonhazardous categories. 
Hazardous near misses are defined as those where 

r '' 

1) Time was not available for an evasive maneuver and avoidance of an 

accident was strictly by chance - ■ . . > 

2) Collision might have resulted if evasive action had not been taken by 
. either pilot. 

In the past, near misses frequently were not reported by the pilot in order to avoid 
action under the FAA enforcement program. To obtain as many reports as possible, 
plus an objective analysis of this situation, immunity to prosecution was granted by the 
FAA administrator to all pilots or controllers involved in a near mid-air collision. 
Results of the study are shown in Figure 1-2. Of the 2230 reported incidents, 1102 
were judged to be nonhazardous and are not considered in this study. Of the 1128 haz- 
ardous incidents, 36% occurred in the en route area and 64% in terminal areas. In 
Figure 1-3 these incidents are compared to the expected value based on the number of 
aircraft operations for each flight rule category. Predictions for the en route area are 
based on total flight hours. Correlation between actual and predicted values for the 
VFR versus VFR case is extremely good. A higher degree of interference between IFR 
and VFR aircraft is experienced than would be predicted. The effect of the control 
system in reducing incidents is clearly shown in the IFR versus IFR case. Returning 
to Figure 1-2, controller error is responsible for four times as many IFR versus IFR 
incidents as pilot error. 

Sixty-four percent of the hazardous incidents occurred in the terminal area. The 
majority of the terminal incidents (87%) occurred at the 322 airports equipped with an 
FAA tower, even though total operations are almost equally divided between tower and 
nontower airports. There are two major reasons for the low percentage of near mid- 
air collisions (NMACs) at uncontrolled airports. First, the private pilot is reluctant 
to report a near miss because of the inconvenience and paperwork involved. Secondly, 
the private pilot has much higher tolerance for close proximity passing since this is a 
frequent occurrence in and around the light plane airport during VFR operations. 

On the next level of Figure 1-2, terminal incidents are classified by flight rule, 
and in Figure 1-3 are compared to the expected value. For the terminal area, the 
number of operations is used as the basis for the expected number of incidents rather 
than flight hours. IFR versus IFR incidents are only 3%, less than half that which 
would be predicted. Incidents between IFR and VFR aircraft account for 63% of the 
total, more than double the number expected. A high degree of interference between 
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Incidents (%) 


ifr/ifr: ifr/vfr vfr/vfr ifr/ifr ifr/vfr vfr/vfr 

En route Terminal 



Figure 1-3. Near Mid-Air Collisions — Involvement by 
Flight Rule En Route and Terminal vs Predicted 
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aircraft operating at major commercial hubs and nearby general aviation airports is 
indicated. VFR/VFR incidents accounted for only 34% of the total, approximately half 
the number expected. The terminal ATC task is generally thought to be more difficult 
than the en route task since airspace is restricted, aircraft traffic densities higher, 
and aircraft altitude is frequently changing. Yet, the percentage of total incidents oc- 
curring between two IFR aircraft in the terminal area is one-quarter that of the en 
route airspace. 

1. 1. 1 GEOGRAPHIC DISTRIBUTION 

A geographic distribution of mid-air collisions for the years 1966 through 1968 is 
shown in Figure 1-4. Although they tend to group in California, Florida, Ohio, and 
New England, collisions are also distributed over low traffic areas. The problem is 
better defined by a plot of the near mid-air collisions shown in Figure 1-5. Here the 
incidents are clearly shown clustering around major hubs such as New York, Los 
Angeles, San Francisco, and Washington. Smaller clusters are shown over the train- 
ing areas in South Florida, Pensacola, and Phoenix. Incidents occurring away from 
population centers fall to a great degree on airways and, more specifically, on inter- 
sections of airways. 

An idea of the concentration of NMACs in the area of major hubs may be obtained 
from Table 1-3. Total operations at airports within a 50 mile radius of the hub have 
been estimated*. NMAC incidents occurring in the same area are taken from the 1968 
report and the NMAC rate per million operations is determined. The hubs are then 
listed by incident rate in decreasing order. It is readily apparent that factors other 
than total number of operations in the area determine the incident rate. Phoenix and 
Pensacola have relatively low traffic counts but a high proportion of military training 
flights. The incident rates for New York, Philadelphia and Washington are similar and 
quite high. San Diego, San Francisco, and Los Angeles possess comparable geographic 
features and have similar but lower incident rates. No explanation is offered for the very 
low rates found at Dallas, Miami, Minneapolis, Columbus, and New Orleans. 

Distance from terminal at which incidents occurred was also investigated. Number 
of actual and near mid-air collisions is plotted as a function of distance from the terminal 
in Figure 1-6. Since many more near mid-air collisions occur than actual mid-air colli- 
sions, a relative frequency scale is utilized. Actual mid-air collisions peak roughly 2 
miles from the airport while the near mid-air collisions peak slightly more than 3 miles 
from the terminal. This difference is partially explained in Figure 1-7 which shows the 
location of 20 mid-air collisions which occurred in 1988 during the landing phase. As 
the airport is approached, density of traffic increases and the pilot is alert, so the num- 
ber of incidents reported is high. The majority of accidents, however, occur on or near 
final approach. Four factors appear to be responsible 

1) When on the base leg and final approach, the pilot is preoccupied with 
landing procedures. 

* Operations at tower equipped airports in the area multiplied by the ratio of total 
operations to operations at tower equipped airports. 
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Figure 1-4. Mid-Air Collisions 1966 through 1968 



Figure 1-5. Special Near Mid-Air Collision Study, 1968 
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Table 1-3 


NMAC RATE FOR MAJOR HUBS 



Operations 1 
(millions) 

NMAC 

Rate 

Per Million 

Phoenix 

0.7 

12 

17.1 

Pensacola 

0.3 

5 

16. 6 

New York 

3.2 

53 

16.6 

Philadelphia 

0.9 

15 

16.6 

Washington 

1.5 

24 

16.0 

Denver 

1.1 

17 

15.5 

San Diego 

1.5 

19 

12.6 

San Francisco 

3.2 

38 

11.9 

Los Angeles 

6.3 

74 

11.7 

i 

Cleveland 

0.6 

6 

.11.6 

Chicago 

1.9 

21 

11.0 

Boston 

0.6 

6 

10.0 

Norfolk 

0.7 

6 

8.6 

Kansas City 

1.3 

11 

8.5 

Seattle 

1.4 

11 

7.8 

Dallas 

1.7 

10 

5.9 

Miami 

1.8 

10 

5.5 

Minneapolis 

1.7 

9 

5.3 

Columbus 

1.3 

6 

4.6 

New Orleans 

0.9 

3 

3.3 

- 

■^Estimated local and transient operations in 50 mile radius area not 
including military operations. 
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2) The airspace through which the flow of aircraft must pass is enormously 
compressed as the end of the runway is approached. 

3) Visibility .of the aircraft is poor in the downward direction. 

4) Due to the angle of bank, aircraft turning a short base leg have poor 


visibility of any aircraft already on final. 



Approximate Position of 20 Mid-Air Collisions that Occurred in the Down-Wind, 
Base- Leg and Final Approach. 


Source: Mid-Air Collisions in US Civil Aviation-1968, NTSB, July 1969 
Figure 1-7. Mid-Air Collisions During Landing Phase 

1.1.2 OPERATIONAL FACTORS 

As expected, the phase of operation significantly affects the number of mid-air 
collisions. Total collisions for the period 1966 to 1968 are plotted by phase of opera- 
tion in Figure 1-8, with each aircraft plotted separately. The predominance of acci- 
dents in the final approach and traffic pattern have already been discussed. One would 
intuitively assume that climbing and descending aircraft are a major factor in mid-air 
collisions. This does not appear to be the case since the bulk of the nontraffic pattern 
accidents occurred with both aircraft under normal cruise conditions. 

Angle of closure is also a significant factor. In Figure 1-9 near miss incidents 
are compared against actual collisions for both en route and terminal air space. As 
before, there are many more NMAC incidents than actual collisions so data is plotted 
on the basis of percentage of total. For the en route area, the majority of the NMAC 
incidents occur from ahead or from crossing traffic. Actual collisions follow an in- 
verse pattern with the majority of the collisions taking place between overtaking traffic. 
The same pattern is evident but even more pronounced for the terminal area. Poor 
cockpit visibility adds to the visibility problems created by weather and sun glare. 
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Figure 1-8. Mid-Air_Collisions by Phase of Operation, 1966 to 1968 
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Figure 1-9. Angle of Closure 
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Closing speed between the two aircraft involved in near miss also seems to be a 
major factor in mid-air collision situations. Unfortunately, data is available only on 
NMAC incidents. Closing speeds between the conflicting aircraft are plotted in Figure 
1-10. As expected, speeds are lowest in the terminal area and increase through the 
transition area, and further increase in low altitude en route and are greatest in high 
altitude en route. Two factors are significant 

1) Most low altitude incidents occur in the vicinity' of 300 knots, a speed 
range where see-and-avoid is estimated (Reference 10) to be 95% 
effective. 

2) Above 10, 000 feet the majority of incidents occur above 400 knots 
where the estimate of the effectiveness of see-and-avoid has fallen 
to 50%. 



Figure 1-10. Closing Speed at Time of Near Mid-Air Collision 

_v 

The role of radar advisories in NMAC incidents is shown in Figure 1-11. In the 
terminal area 85% of the incidents occurred within radar coverage, and in 60% of the 
incidents one or both aircraft was in contact with an ATC facility. Radar advisories 
were issued in only 21% of the cases due to controller workload and other reasons. 
Similar figures are shown for en route airspace. Here radar advisories alerted the 
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Figure 1-11. Relationship between NMAC and Radar Advisories 
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pilot to the near miss situation less than 8% of the time. It appears that the usefulness 
of radar advisory information is limited by the control system’s ability to issue advisor- 
ies and communicate with aircraft rather than lack of radar coverage. 

1.1.3 ALTITUDE DISTRIBUTION 

The final and perhaps most important aspect of the mid-air collision is the distri- 
bution of these incidents by altitude. Figures 1-12 and 1-13 are altitude distribution 
of near- miss incidents taken from the FAA NMAC report. Incidents in the terminal 
areas occur at the lower altitudes and IFR/IFR incidents are almost negligible. Inci- 
dents between VFR aircraft occur lower than those between mixed traffic. En route 
incidents between two controlled aircraft are principally a phenomenon of the higher 
altitudes peaking at 15,000 feet. Although these distributions are as one would expect, 
the main question is how they compare with aircraft population at each altitude. Official 
figures on total aircraft population as a function of altitude do not exist so data from a 
number of sources (References 6 and 7) has been brought together into a composit esti- 
mate. The number of flight hours at each altitude has been used for comparison purposes 
rather than number of flights. Estimated distributions are shown in Figures 1-14 and 
1-15. Figure 1-14 is a cumulative plot of the total hours by class of user. General 
aviation is obviously the dominant factor with the bulk of its hours occuring below 
10, 000 feet. Figure 1-15 is that of IFR and VFR flight hours. Note that VFR opera- 
tions extend all the way to 24, 000 ft and that both IFR and VFR hours peak below 6000 ft. 



Figure 1-12. Type of Flight by Altitude Distribution (Terminal - 719 Incidents) 
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Number of Incidents 

Source: Near Midair Collision Report of 196$ FAA, July 1969 

Figure 1-13. Type of Flight by Altitude Distribution ( En Route - 409 Incidents ) 

One would intuitively expect altitude transitioning of jet traffic to be the cause of 
the high incident rate at intermediate altitudes. Yet, earlier charts, Figures 1-8 and 
1-9, show that incidents occur between aircraft in level flight. There appears to be 
several reasons for this paradox with altitude transition an indirect cause. Altitude 
maneuvers are complex but relatively safe for controlled aircraft (IFR) since they are 
closely watched by one or more controllers. Thus, during these maneuvers the con- 
troller's attention is somewhat diverted from other aircraft under this control. Atten- 
tion is also diverted to aircraft at other altitudes. Since the bulk of the air traffic is 
located below 9000 feet the controller is concentrating on providing separation in these 
higher density areas. The third factor is speed. At the time of this study, aircraft 
below 10, 000 feet were restricted to 250 knots. Many aircraft routinely descend to the 
10,000 feet level at much higher velocities. These higher speeds create problems both 
in time^yailable for_the .controller to detect and react to4he=closing situation, and in 
time to plan and coordinate handoffs. The majority of the low altitude controller's ac- 
tivity involves lowspeed traffic and he apparently does not react well to the high-speed 
situation. 
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Figure 1-14. Total Hours Cumulative Distribution by Altitude and User 
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Figure 1-15. Total Hours Cumulative Distribution by Altitude and Flight Rule 
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1.1.4 CONCLUSIONS 


1) This analysis of MAC and NMAC statistics indicates that the CONUS 
airspace can be divided into five regions, each of which has unique 
problems and requirements: 

a) High altitude en route 

b) Major air carrier terminals ' 

c) Low altitude en route 

d ,) High density general aviation airspace 
e) Nontower general aviation terminal airspace 

2) The principal hazard in the terminal area is the mixture of controlled 
(IFR) and uncontrolled (VFR) or high and low performance aircraft. 

3) The present airway navigation system concentrates the traffic into a 
small percentage of the airspace significantly increasing the collision 
hazard, especially at the airway intersection points. 

4) Visual flight rules are not an effective method of separation of aircraft 
at speeds in excess of 250 knots and at altitudes in excess of 10, 000 ft. 

5) A large percentage of the near mid-air collisions take place within 
radar coverage; but the number of radar advisories issued is small. 

Lack of communication between the pilot and the control facility 
reduces the potential effectiveness of radar advisories. 

6) Restricted cockpit visibility is the major cause of mid-air collisions 
at general aviation airports. 

7) In a controlled system low activity, and low stress periods are 1 not 
necessarily the safest since pilot and controller inattention becomes 
a factor. 

8) As a possible corollary to item 7, in the en route phase of flight, 
controller error is the principle source of near mid-air collisions 
involving IFR/IFR aircraft. 

1.2 REQUIREMENTS FOR ATC COMMUNICATION TO AIRCRAFT 

One of the current near critical problems in ATC is in the area of ground to air- 
craft communications. There is a capacity problem because of the shortage of available 
VHF channels and high channel utilization under peak loading conditions. There is also 
a workload problem for both controller and pilot imposed by the need to exchange essen- 
tial information. The RTCA, the FAA, and others have recognized that an automatic 
data link system for ATC can provide significant relief in both the capacity and work- 
load areas. The purpose of this subsection is to develop data rate requirements for 
such a system. 
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1. 2. 1 TERMINAL AREA 

Terminal ATC is generally concerned with either the first 15 minutes of flight or 
the last 15 minutes depending on whether a departure or arrival operation is required. 

In Table 1-4, a set of ATC commands for terminal operations is presented. An estimate 
of the frequency of command issuance is given (Reference 11) along with an estimate of 
the message length. The message length includes the overhead items listed in Table 1-5. 
An 8-bit character was selected to be compatible with the RTCA standard (Document DO- 
136) and ARINC/AEEC design. More efficient coding could be employed from the stand- 
point of the rate of information exchange, but at a corresponding loss in error checking 
capability. A missed approach sequence is described in Table 1-6. 

Referring to Table 1-4, there are several assumptions implicit in the frequency with 
which various commands are given. These assumptions are 

1) That a flow control capability will exist which will tend to reduce the 
frequency of holds 

2) That track assignments will consist of coded routing like SIDs and 
STARS 

3) That communications channel switching will be primarily a ground 
function. 

The product of the command message length and the frequency of command occur- 
rence yields the weighted message length shown in Table 1-7. The totals shown in the 
arrival and departure columns are an estimate of the total number of bits of data which 
must be transmitted per aircraft. 

1.2. 2 EN ROUTE AREA 

The en route ATC is concerned with the control of aircraft from departure terminal 
handoff to arrival terminal handoff. An estimate of ATC en route requirements essenti- 
ally follows the same procedure used for terminal requirements. The en route ATC 
commands are shown in Table 1-8, with message length and frequency of occurrence 
(Reference 12). The frequency of command occurrence for route, altitude and speed 
changes is based on pfeak day statistics for the Chicago ARTCC, December 1968. The 
weighted message lengths are shown in Table 1-9. 

1.2.3 DATA RATES 

"TKsH dataratFsreq^^ fimctidnOf traff i^ loading, in particular the peak load- 

ing. . The data rate computations are shown in Table 1-10 for both the en route and termi- 
nal areas. The en route estimates are based on the ATCAC (Reference 13) load of 2500 
IFR aircraft. The traffic distribution is based on the Chicago ARTCC statistics. One 
way data rate estimates are 1785 bps. Assuming that aircraft acknowledgement is accom- 
plished by parroting the message to the ground, the two way data rate will be 3570 bps. 
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The terminal data rates are based on the ATCAC (Reference 13) Los Angeles 
terminal basin peak load for 1990 of 1160 aircraft/hour arriving and departing from 
several airports. The one way data rate is 1513 bps and the two way rate is 3026 bps. 

1.3 REFERENCES 

1 0 "Control Subsystem for the f Next Generation 1 Air Traffic Control System", 
IBM IRAD Study, July 15, 1971. 

2. Briefs of Accidents Involving Mid-Air Collisions, U. S. Civil Aviation 
1966, National Transportation Safety Board, June 1968. 

3. Briefs of Accidents Involving Mid-Air Collisions, U. S. General Aviation, 
National Transportation Safety Board, 1967. 

4. Mid-Air Collision in U. S. Civil Aviation 1968, National Transportation 
Safety Board, July 1969. 

5. Near Mid-Air Collision Report of 1968, Federal Aviation Administration, 

July 1969. 

6. En route IFR Air Traffic Survey Peak Day FY 1969, Federal Aviation 
Administration. 

7. Traffic Patterns for IFR and VFR Aviation, Calendar Year 1969, Federal 
Aviation Administration. 

8. FAA Air Traffic Activity FY 1970, Federal Aviation Administration, 

August 1970. 

9. Air Traffic Control Advisory Committee Vol I, Department of Transportation, 
December 1969. 

10. Reparation of Air Traffic by Visual Means: an Estimate of the Effectiveness 
of See and Avoid", Proceedings of the IEEE, March 1970. 

11. Fred Holland, Mitre, WP 8266, "Computer Sizing of Terminal Area Command 
and Control for the ATC Advisory Committee, " April 11, 1969. 

12. Fred Holland, Mitre, WP 8272, "Computer Sizing of En route Command and 
Control for the ATC Advisory Committee, " April 22, 1969. 

13. f1 Report of Department of Transportation Air Traffic Control Advisory 
Committee, " December 1969. 


1-21 



TERMINAL UP LINK COMMANDS 
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* Includes message overhead (see Table 1-5). 

** Assumes the implementation of some form of the discrete frequency concept in which the ground 
system performs most of the frequency switching. 

*** Probability of a missed approach is assumed to be 0.01 (see Table 1-6 for details). 


















Table 1-5 


DATA LINK MESSAGE OVERHEAD 


Component 

Characters 

Bits 

Synchronization 

2 

16 

Start of Header 

1 

8 

Aircraft Address 

7 

56 

Ground Station Address 

4 

32 

Mode 

1 

8 

Start of Text 

1 

8 

End of Message 

1 

8 

Block Check 

1 

8 

TOTAL 

18 

144 


Table 1-6 

UP LINK COMMANDS FOR AIRCRAFT EXECUTING 
A MISSED APPROACH SEQUENCE 


Command 

Number 

Message Length* 

Bits /Sequence 

Char 

Bits 

Execute Missed Approach 

1 

19 

152 

152 

Take Up Heading 

8 

20 

160 

1280 

Change Speed 

4 

20 

160 

640 

Change Altitude 

2 

20 

160 

320 

Enter Hold 

1 

35 

280 

280 

Leave Hold 

1 

35 

280 

280 

Intercept Landing System 

1 

24 

192 

192 

Total 

18 

— 

— 

3144 

* Includes message overhead (see Table 1-5). 
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Table 1-7 

TERMINAL UP LINK MESSAGES 


Commands 

Weighted Message Length (bits) | 

Arrivals 

Departures 

. Enter Hold 

70 

- * 

Leave Hold 

70 

- 

Track Assignment 

192 

192 

Change Speed 

640 

320 

Take Up Heading 

690 

320 

Change Frequency 

176 

176 

Intercept Landing System 

192 

- 

Change Altitude 

640 

320 

Position for Takeoff 

- 

224 

Cleared for Takeoff 

- 

176 

Missed Approach Sequence 

31 

- 

Total 

2971 

1728 


Table 1-8 

EN ROUTE UP LINK COMMANDS 


Commands 

Message Length* 

Frequency /Flight Category 

Chars 

Bite 

Arrivals 

Departures 

Overs 

Withins 

Enter Hold 

35 

280 

- 

- 

- 

• 25 

Leave Hold 

35 

280 

- 

- 

- 

.25 

Route Assignment 

78 

624 

- 

1.0 


1.0 

Change Route ** 

38 

304 

.2 

.3 

.2 

• 4 

Change Speed** 

20 

160 

.8 

— . 

- 

• 2 

Take Up Heading 

20 

160 

3.0 


- 

6.0 

Change Frequency*** 

22 

176 

2.0 

mM 

- . 

2.0 

Change Altitude** 

20 

160 

3.9 


.8 

3.1 

Total 



9.9 

7.8 

1.0 

13.20 


♦ Includes message overhead (see Table l-5) ' — ' "" — 

** Command frequency Is based on peak day statistics for Chicago ARTCC, December 1968 

*** Assumes the implementation of some type of discrete frequency concept in which the ground 
system performs most of the frequency switching. 
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Table 1-9 


EN ROUTE UP LINK MESSAGES 


Commands 

Weighted Message Length (bits) | 

Arrivals 

Departures 

Overs 

MJU.II.fcM 

Enter Hold 

— 

— 



70 

Leave Hold 

- 

- 

- 

70 

Route Assignment 

- 

624 

- 

624 

Change Route 

61 

91 

61 

122 

Change Speed 

128 

- 

- 

32 

Take Up Heading 

480 

480 

- 

960 

Change Frequency 

352 

352 

- 

352 

Change Altitude 

625 

240 

128 

496 

Total 

1646 

1787 

189 

2726 


Table 1-10 
DATA RATES 



Arrivals 

Departures 

Overs 

Withins 

En route (assuming 2500 IFR Aircraft) 

Message Length/Aircraft 
(bits) 

1646 

1787 


189 

2726 

Traffic Distribution (%) 

34.9 

24.8 


16.3 

24.0 

Traffic Count 

873 

620 


407 

600 

Control Life (minutes) 

33.1 

40.2 


59.3 

46.8 

Aire raft /Hour 

1580 

925 


412 

769 

One Way Data Rate 
(bits/s) 

722 

459 


22 

582 

Total One Way Rate = 1785 bits/s 



Arrivals 

Departures 

Terminal 

Message Length/Aircraft 
(bits) 

2971 


1728 


Aircraft/Hour 

1160 


1160 


One Way Data Rate 
(bits/sec) 

956 


557 


Total One Way Rate =1513 bits/s 
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COMMUNICATION TECHNIQUES 
FOR SATELLITE ATC SURVEILLANCE 


2.1 INTRODUCTION 

The communications connected with a COLM Air Traffic Control System using 
satellites are diverse, covering a wide range of operational requirements. A pre- 
liminary categorization which can be used to focus the problem and direct attention 
toward specific topics which are peculiar to satellite systems., and which have not 
been systematically treated in the pas t is the following: 

1) Voice communications for advisories and control 

2) Communications specific to surveillance position location 

3) Communications specific to navigation position location 

4) Digital data communications other than those required for navigation and 
surveillance (e.g. , company data transmission, weather advisories, etc. ) 

5) Communications related to satellite control, (e.g. , telemetry and house- 
keeping, ephemeris determination, and satellite electronics timing and syn- 
chronization) 

The voice communication problem is not peculiar to satellite ATC systems, 
although the inefficient use of bandwidth, which is characteristic of voice, makes it 
even more unattractive for satellite relay (with its premium on bandwidth) than for 
line-of-sight communications. The trend, in, present planning, is to replace voice 
transmissions with formatted binary data links. We have correspondingly assumed 
that such voice links as remain in the satellite-based system will not be relayed over 
satellite links, and hence do not constitute a problem appropriate for detailed con- 
sideration in this study. 

The satellite control communications of item 5 take place over links between the 
satellites and fixed earth stations. Such links are characterized by high-gain antennas 
at the ground, and are well understood, having been the subject of attention in many 
contacts over the past 10 or 15 years. They are, therefore, also not considered in 
detail. 

Items 2, 3, and 4 involve communication paths between the ground control center 
and the airborne users which use the satellite repeater as a relay point. The 
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ground-to-satellite segments of these links are not considered critical for the reasons 
cited above. The satellite-to-user links, on the other hand, are characterized by 
low-gain nondirectional aircraft antennas, and by a large number of users. 

Thus, although the data rate from the control center to each user may be quite 
low, the satellite-to-user downlink must carry a high data rate, with correspondingly 
stringent requirements on satellite prime power. For example, 10-bit/second user 
implies, for 10 5 users, a one megabit per second downlink data rate. In the other 
direction, the return of data from user to control center, which cannot be avoided for 
surveillance, imposes a severe multiple -access problem, with all users requiring 
simultaneous access to the satellite repeater. 

From an operational point of view, the most attractive feature of a satellite-based 
ATC system is its ability to transcend line-of-sight limitations, and to provide sur- 
veillance coverage for all users without low-altitude blind spots. Because of the 
potentially low cost to the user of equipment to provide surveillance, the surveillance 
aspect is considered to be more critical than the area navigation capability that satel- 
lites can provide. Since navigation is inherently a high user cost item, it will proba- 
bly be limited to large commercial and military users in any practical system. 

Consequently, in this section we concentrate primarily upon the surveillance 
communication problem, and examine in detail the satellite-to-user link. To do this 
in a systematic way, we first postulate a number of satellite configurations which are 
sufficiently general to encompass almost all conceivable surveillance and navigation 
systems. 

These configurations are shown in Figure 2-1, and are briefly described as 
follows: 

1) Configuration I - A number of satellites transmit synchronized ranging wave- 
forms to the users. Users with the appropriate equipment can use these 
signals for area navigation, but, in general, the users will make time-of- 
arrival difference measurements and store this data temporarily. According 
to some multiple-access and polling procedure, each user transmits this 
data through a single relay satellite to the control center, where the measure- 
ments are converted to a three-dimensional position fix by hyperbolic compu- 
tation. 

2) Configuration II - Each user autonomously transmits a periodic ranging 
waveform through at least four satellites to the control center. Time-of- 
arrival difference measurements are made for each user, and a position fix 
is computed hyperbolically. 

3) Configuration m - This differs from Configuration II only in that the ranging 
waveforms are elicited by an interrogation waveform which polls each user 
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at a specified time. The position fix is then computed by trilateration using 
the two-way range measurements. 

4) Configuration IV - This configuration is characterized by the use of a single 
satellite equipped with crossed-boom interferometer. Suitable waveforms 
are transmitted by the users either autonomously or by interrogation, and 
measurements of phase differences of the signals transponded by the individu- 
al interferometer antennas permits a two-dimehsional position to be deter- 
mined. The third dimension is supplied by the user, who includes his pres- 
sure or radar altimeter reading as part of his transmission. 

In analyzing candidate approaches to the communications aspects of each of these 
configurations we take a top-down approach which emphasizes the multiple-access and 
polling problems rather than the details of specific modulation methods. It is the 
decisions made at these higher levels of system design that have the greatest effect 
upon the efficiency of the system in terms of bandwidth and power. Specifically, we 
recognize four levels of design which are described in the following subsection. 

2. 1. 1 ELEMENTS OF SATELLITE COMMUNICATIONS DESIGN 

The elements of satellite communications design for ATC surveillance and navi- 
gation applications can be viewed as consisting of four separate but related levels of 
design considerations. These considerations are (1) first level, polling; (2) sec- 
ond level, multiple access; (3) third level, modulation; and (4) fourth level, ranging. 

1) First Level - Polling Alternatives 

In synthesizing satellite ATC surveillance systems, the most fundamental 
design decision involves how the operation of polling users is accomplished. 
Basically there are three alternatives, namely, 

a) Autonomous - Asynchronous user generation of position messages are 
exclusive of external control signals or precise timing standards. 

b) Elicited - User position messages are generated in response to central 
control command or to a central control common timing signal. Group 
polling may be accomplished, for example, by airspace areas or by 
user aircraft performance categories, etc. 

c) Hybrid; Individual, Group, Common Timing Control - Hybrid techniques 
utilize advantageous features of both autonomous and elicited polling. 

That is, such criteria as aircraft performance, flight profile, and air- 
space density facets may be utilized to determine transitional points 
between autonomous and elicited techniques to facilitate an optimum 
position update rate. Under processor control, the polling technique 
could be altered in real time to emphasize surveillance of high risk air- 
space areas. 
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2) Second Level - Multiple Access Alternatives 

These alternatives involve the methodology of obtaining signal orthogonality 
to permit surveillance of many users within the required interval. These 


alternatives* are 



a) TDMA 

TDMA/CDMA 

TDMA/FDMA/ CDMA 

b) FDMA 

TDMA/FDMA 


c) CDMA 

CDMA/FDMA 



d) All the above with SDMA 

3) Third Level - Modulation 

The third level of investigation, modulation, essentially impacts communica- 
tion link performance in terms of detection sensitivity, error rate, and infor- 
mation bandwidth. As such, the analysis for a given application is straight 
forward. However, when combined with conjunctive ranging and reporting 
configurations, additional modulation requirements exist which impact trade- 
off selection. 

4) Fourth Level - Ranging 

In those configurations where user-initiated signals are utilized by the control 
center to determine position, or, as in Configuration I, where satellite rang- 
ing waveforms are used for position location measurements, the impact of 
ranging signals must be considered with communications requirements. 
Ranging signal structures considered include 

a) Sidetone 

b) BENOR** 

c) Pseudo-Noise 

d) Pulse 

e) Pulse Compression 

1. Analog 

2. Digital 

2. 1. 2 PRINCIPAL SYSTEM PARAMETERS AND ASSUMPTIONS 

In the analyses performed in this section, the basic approach is parametric in 
order to allc^th^ widest posslbre application of the results. Where specific calcula- 
tions are made, the numerical assumptions listed on the next page are used. 


♦TDMA 

FDMA 

CDMA 

SDMA 

♦♦BINOR 


Time Division Multiple Access 
Frequency Division Multiple Access 
Code Division Multiple Access 

Space Division Multiple Access (multiple antenna beams) 
Binary Optimum Ranging 
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1) Number of system users (max), N, 

2) Update interval, T , 

3) Total message length (Configuration I), T., 

4) Desired data quality, P^, 

5) Signal-to-noise (thermal-no CDMA) 

6) Signal-to-noise (thermal with CDMA) 

7) Signal-to-Clutter (CDMA) 

8) Data transmission efficiency 


1 second 
100 bits 

' 1(T 5 bit error rate (BE R) 
13 dB for 10" 5 BER 
20 dB for 10~ 5 BER 
13 dB for 10 -5 BER 
1 bit/Hz 


Aircraft-to-satellite link budget is as follows: 


1 ) 

2 ) 

3) 

4) 


User antenna gain 
Satellite antenna gain 
Path loss (@ 1.6 GHz) 
Nq KT (T=910°K) 

P /N 

R 0 W 


0 dB 

24 dB (CONUS coverage) 
187 dB 
-199 dB 

36 dB/ Hz for each watt 


2. 1. 3 PRINCIPAL DESIGN PERFORMANCE PARAMETERS 


The design performance parameters employed for trade-off comparison relate 
to either the use of basic resources such as bandwidth and spectrum, or to practical 
considerations which ultimately limit feasibility. These principal performance 
parameters are , 

1) Bandwidth efficiency (users per megahertz) 

2) Power/BER (channel quality per watt) 

3) Signal process gain required 

4) Number of spacecraft required by approach 

5) Number of spaceborne antenna per spacecraft 

6) Ground center complexity 

7) Requirements for return ground to user link (via satellite) 

8) Return link requirements, e.g. , power, speed, BER. 
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2.1.4 OUTLINE OF SECTION 2 

The next two parts of this section deal with topics which have application to all 
configurations. Subsection 2.2 presents a trade-off analysis of candidate waveforms 
for ranging. Subsection 2. 3 considers the problem of sequential polling and the design 
of available polling waveforms for each of the satellite configurations. In subsections 

2.4 through 2. 7 the multiple access problem is considered in detail for each of the four 
configurations. Each multiple-access candidate is analyzed using a unified approach 
as far as possible. The results are presented in a way which is convenient for com- 
parison of candidates. Subsection 2. 8 consists of a summary which includes a com- 
parison of the approaches analyzed and a preliminary selection of a number of candi- 
dates which appear suitable for further consideration. 

2.2 RANGING WAVEFORM ANALYSIS 


This subsection presents a comparative analysis of a number of ranging waveform 
designs. The analysis is applicable to all the configurations, although it is directed 
primarily toward the requirements of Configuration I, in which navigation signals are 
processed by the user for later retransmission. 

2.2.1 COMPARISON CRITERIA 

The comparison is made on the basis of the following four criteria: 

1) Ranging power utilization 

2) Bandwidth occupancy 

3) Attainable fix rate 

4) Cost and complexity of user hardware. 

The principal operational criteria are ranging accuracy and fix rate, which are 
linked by the noise bandwidth of the phase-locked loop used for phase measurement. 
The required loop bandwidth is in turn determined by the received power-to-noise 
density ratio, the loop center frequency, and the system dynamics. In the following 
analyses, ranging accuracy is normalized to the ranging signal bandwidth, which is 
assumed to be determined by the chip rate in binary codes and the upper tone fre- 
quency in sidetone signals, and found as a function of Pr/N 0 , the signal power-to- 
noise density ratio at the receiver. The attainable fix rate is taken to be the recipro- 
cal of the ranging phase-locked loop bandwidth. 

The cost of the hardware required by the user in each of the candidate techniques 
is influenced by several factors. Signals must be received from approximately four 
satellites for each fix, and these signals must access the user’s receiver according to 
some multiple-access discipline. Frequency-division, code-division, and time- 
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division multiple access methods are considered where applicable, and the implica- 
tions of each method are determined in terms of number and complexity of the receiv- 
ers required, bandwidth utilization, and fix rate. 

2.2.2 BINARY CODES 

2. 2. 2.1 Continuous PN 

The basis of this waveform is a PN sequence which is transmitted continuously 
by each satellite. In FDM systems each satellite may use the same code, while in 
code division each satellite will use a different code. For what follows, we initially 
assume an FDM system. 

The fine ranging measurement is made by measuring the phase of the receiver 
delay-lock loop clock relative to an on-board clock. Practical restrictions limit the 
accuracy of phase measurements to about 0. 01 cycle, which must correspond to the 
best desired rms error, which is assumed to be about 50 ft. This indicates a mini- 
mum code chip rate of 400 kHz. The range ambiguity requirements dictate a code 
repetition rate not exceeding 100 Hz. The limits for both accuracy and ambiguity 
resolution are exceeded by the use of an 8191 bit code, for example, a maximal length 
sequence from a 13 stage shift register. Another well-known code of this length is 
the JPL code (see Reference 1), which is synthesized from three sequences of 7, 19, 
and 31 bits, respectively. Typical chip length and sequence repetition rates for an 
8192 chip code are 500k chips/second and 61 Hz, respectively. 

Continuous transmission of a PN code sequence provides sufficient information 
for the user to acquire the signal, synchronize with the code, and perform ranging 
measurements. Using the code parameters of the preceding paragraph, a represen- 
tation of the code characteristics in time and frequency is shown in Figure 2-2. The 
spectrum shown is of the binary PN code itself. In practice, the baseband code will 
be filtered to remove the high frequency components, which will result in de-emphasis 
of the spectrum lobes above the first null. 

The transmission of data in conjunction with the PN waveform can be handled in ^ 
several ways. The binary data can be modulated directly onto the PN sequence by 
code phase reversal at the appropriate bit rate. The data, which is narrowband, can 
also be placed on a separate subcarrier located above the baseband spectrum of the 
PN signal. The low data rates involved in the satellite transmission make the data 
transmission a fairly uncritical item in the waveform design. Quantitative support 
for this statement is obtained from the link characteristics given in subsection 2. 1.1, 
which shows a received signal -to-noise density ratio of 36 dB/Hz for each watt of 
satellite transmitted power. If 13 dB is assumed to be sufficient E/Nq for a satisfac- 
tory error rate, then -21 dBW /bit/s must be transmitted by the satellite. For the 
present assumption of no more than 10 bits/s, only -11 dBW of satellite power need 
be used for data transmission. 
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fit) 


N = Code Length 


l J I I t 

0 { NA 2NA 

b. PN Code Time Profile 

Figure 2-2. Code Characteristics in Time and Frequency 

Although acquisition of the code signal is only done by each user upon take-off (or 
if synchronization is lost due to noise) it is noted that the use of a composite code such 
as the JPL code permits rapid acquisition, with each code component being acquired 
in turn. For example, an analysis given in Reference 1 gives an acquisition time of 
1. 0 second with a P r /N 0 of 42 dB/Hz. This figure assumes a 99, 9% required proba- 
bility of acquisition. 

2. 2, 2. 1.1 Range Measurement 

The range measurement is performed by measuring the phase of the clock gener- 
ated in the receiver delay-lock loop, while coarse range is determined by the offset 
between thqreceivedcode and^the_,referencecode. gene rat edby.the loop. -ForFDMA 
(i,e. , a separate frequency channel per satellite), we assume coherent demodulation 
of the PN sequence by a delay -lock loop such as that shown in Figure 2-3. The pres- 
ence of several signals in the same channel precludes the use of coherent demodula- 
tion; thus, for CDMA, a noncoherent receiver is required such as shown in Figure 2-4. 
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CDMA performance is discussed in the next subsection. The ranging accuracy is 
determined by the mean-square phase jitter at the YCO output. In a coherent receiver, 
two factors contribute to this jitter: (1) the thermal noise at the input, assumed to be 
Gaussian and white across the IF bandwidth, and (2) the self-noise generated by the 
correlation process, which results from the fact that under locked conditions the input 
is time-displaced from both reference code generator outputs by one-half a chip length. 



Figure 2-3. Coherent Delay-Lock Loop for FDMA Satellite Signals 
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Figure 2-4. Envelope-Correlator Delay-Lock Loop for CDMA Satellite Signals 


Gill (Reference 2) analyzed a delay-lock loop of the type shown in Figure 2-3, and 
derived an expression for the mean-square error resulting from the thermal noise. 
The error normalized to the chip length A (and hence to the RF bandwidth) is given by 
11/2 


A 


B N 
n o 


2P 


(2-1) 


where B_ 


= loop noise equivalent bandwidth 


P /N = received signal -to-noise density ratio, 
r o 

The normalized self-noise of the correlation process is given by (see Reference 3): 


n 

■£ 


AB 


1 + 24 ( 


A 
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(2-2) 


where, e is-the delay measurement errors Under in-lock conditions, c - 0, and 
the following term can be written: 



(2-3) 
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Under in-lock conditions, it is reasonable to assume statistical independence between 
the self-noise and the other error contributions. Therefore, the total normalized 
error becomes 


1/2 



An effect which must be considered in evaluating error performance by means of 
Equation 2-4 is the delay -lock loop threshold. This arises from the shape of the loop 
discriminant curve shown in Figure 2-5, which shows the correlator output as a 
function of delay error. For errors lying in the interval (-A/2, A/2) the discrimi- 
nant is linear. Errors exceeding this range, however, cause the loop to lose lock. 
The chip length should be chosen such that the mean square error is less than about 
0.3 A (Gill, Reference 1). For a minimum range error of 50 ft, this condition is 
satisfied for A > 170 x 10~^ seconds. 


Correlator 

Output 



Figure 2-5. Delay-Lock Loop Discriminant 


Equation 2-4 is plotted in Figure 2-6 for A = 10 seconds, and two values of loop 
noise bandwidth. The performance degradation at high signal -to-noise density ratios 
arises from the correlator self-noise. The threshold of the delay -locked loop is also 
shown at cr = 0.3 A. 
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Figure 2-6. Ranging Accuracy PN-FDMA 
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2. 2. 2. 1.2 Fix Rate 


The fix rate which can be achieved is determined directly by the rate at which 
independent samples can be taken of the loop delay estimate. This depends upon the 
loop bandwidth, and is assumed to be its reciprocal. The curves of Figure 2-6, 
therefore, correspond to fix rates of 1 and 10 per second. In any practical applica- 
tion, a tradeoff must be performed between desired fix rate and error performance. 

2.2. 2. 1.3 Multiple-access Techniques 

As indicated above, the multiple satellite signals can access the user receiver 
either by frequency division or code division. The use of code division will affect 
performance by increasing the interference level. An essential tradeoff to be consid- 
ered is between RF bandwidth and range error. A secondary tradeoff arises from 
consideration of user costs. Each frequency occupied by a separate signal implies 
additional receiver front-end complexity with its attendant cost in space, weight, and 
power. It should be noted that a separate delay-lock loop is required for each satel- 
lite signal regardless of the multiple-access technique used. 

2. 2. 2. 1. 4 CDMA Performance Analysis 

In the foregoing analysis of PN-FDMA ranging, a coherent receiver is assumed. 

In CDMA signaling, the presence of multiple carriers in the same frequency band 
makes coherent reception infeasible. The ranging PN signal must be recovered by 
means of square-law envelope-detector delay-lock loop such as shown in Figure 2-4. 
The normalized error performance of such a system for a single input signal is given 
by (Gill, Reference 2): 


jj 

A 


B N B„B N 2 
n o ^ 1 n o 

2P + 12 + 


AB 
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(2-5) 


where Bj is the bandpass filter bandwidth, assumed to be 10 kHz to account for dop- 
pler uncertainty. This expression differs from that of Equation 2-4 in the addition of 
the middle term. This term arises from the action of the envelope detector, which 
has a signal suppression effect at low signal-to-noise ratios. 

In analyzing the effect of code-division multiple -access upon error performance, 
one can proceed by assuming that the interfering signals are constant-amplitude 
sinusoids of random, uniformly distributed phase angles. The resultant interference 
will be assumed Gaussian and white across the RF bandwidth. Then, if the power in 
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each signal is equal to P r at the receiver, and assuming three interfering signals, the 
effective noise power density becomes 


3P 


N EFF N o + B 


RF 


(2-6) 


The error performance for the code division system is then found by replacing N 0 by 
NgFF * n Equation 2-5. If we let B^ F = l/A , the result is 


a_ 

A 


Note, that for fewer than three interfering signals, the Gaussian assumption can no 
longer be made realistically, and the analysis involves a detailed investigation of the 
cross correlation properties of the codes. In general, no closed form analytical 
procedures are available, however, Gold (Reference 4) and others have completed 
significant work which allows reasonable performance estimates to be made on the 
basis of computer analysis. Equation 2-7 is plotted in Figure 2-7. The deterioration 
of performance from the FDMA case at high received signal power is shown by com- 
parison with the dashed curves, which are identical to those in Figure 2-6. At low 
signal -to-noise ratios, a performance degradation is noted arising from the effect of 
the envelope detector. 


J1-(N + 3P A)+^2.(N + 3P A) 2 + — 
r P 


1/3 


(2-7) 


2. 2. 2. 1.5 Summary 


The principal feature of PN ranging is the ability of several signals to share the 
Same frequency band. In view of the scarcity of the spectrum resource and the inher- 
ently wideband nature of ranging signals, this is a significant asset. 


Ranging performance in noisy environments is dependent upon the presence in the 
signal of high-frequency "fine structure" while the resolution of ambiguities requires 
the inclusion of low-frequency components. These two requirements place a lower 
bound on the code chip rate and a lower bound on the code sequence length. Short 
codes which have lengths close to the lower bound are generally easier to generate 
and to acquire. Longer codes may be desirable in some applications where privacy 
is a factor and the longer acquisition time can be tolerated. Composite or separate 
codes could be used in a. system. which accommodates both militarv and^civilian users 
with the same satellites. 

For a given satellite transmitted power, the ranging accuracy and fix-rate 
achieved are determined essentially by the delay-lock receiver bandwidth. This, in 
turn, is determined primarily by the Doppler rate that is expected to be encountered. 
Reasonable aircraft performance assumptions lead to practical loop bandwidths of 1 
to 10 Hz. 
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Normalized rms Error 



Figure 2-7. Ranging Accuracy PN-CDMA and PN-FDMA 
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2.2. 2.2 BINOR Code 

In order to reduce the acquisition time required for long binary PN codes, Stiffler 
(Reference 5) has developed a code referred to as the binary optimum ranging (BINOR) 
code, which is distinguished by the characteristic that the acquisition time required 
in the receiver requires on the order of log 2 N trials for a code of length N, whereas 
a PN receiver normally can require up to N trials. The code is synthesized by ma- 
jority logic from a set of log2 N square waves of harmonically related frequencies. 
Acquisition is accomplished by noting the polarity of each of log 2 N successive partial 
correlations with the component square waves. The spectrum and time profiles of the 
BINOR code are essentially equivalent to those of PN-FDMA as fiven in Figure 2-2. 

A typical receiver structure is shown in Figure 2-8. 



Figure 2-8. Binary Code Receiver 

The transmission of binary data using the BINOR waveform cannot be done by code 
phase reversal because of the acquisition procedure, which is sensitive to code polari- 
ty. Data transmission must be by separate carrier adjacent to same frequency band. 
The same calculations of power required for data transmission as were made for PN 
also hold in this case. 
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2 . 2 . 2 . 2 . 1 Phase Measurement 


Once the code has been acquired by the user, the phase measurement is made on 
the lowest frequency square wave, which is derived from the recovered bit clock by a 
divider chain. The phase jitter on the low-frequency square wave is thus equal to the 
bit clock recovery phaselocked loop. The normalized mean square delay error result- 
ing from this jitter is given by the expression , 



1 _ 
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B N 
n o 


2P P 2 
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where B n is the equivalent noise bandwidth of the bit-clock phase-locked loop, 
and P is the correlation coefficient between the clock and the received signal, which 
is given by 



f(log 2 N) t r 


(2-9) 


Since, from Equation 2-9, the correlation coefficient is dependent upon log 2 N, the 
ranging accuracy is a function of the code length N. Hence, a lower bound on accuracy 
is imposed by the ambiguity resolution. Assuming, for equivalence with the PN code, 
a chip length of 1 /i.s and a code length of 8192 chips; then 13 square waves are 
required for code synthesis. Equation 2-8 is plotted in Figure 2-9 for log 2 N = 13 
and typical values of loop noise bandwidth. 


,2. 2. 2.2.2 Fix Rate 


As in the case of PN ranging, the attainable fix rate is dependent upon the recov- 
ery loop bandwidth. 


2. 2. 2. 2. 3 Multiple-access Methods 


Due to the uniqueness of the BINOR code, code-division multiple-access cannot 
be used. Frequency separation of the signals is the only viable method. Separate 
code receivers are required for each satellite, and sufficient bandwidth for four 
simultaneous transmissions must be available. 

2.2. 2.2. 4 Summary 

The BINOR code differs from PN codes primarily in its regular structure which 
permits rapid acquisition. The regular structure, however, also implies a loss of 
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the desirable autocorrelation characteristics which permit the use of a delay-locked 
loop receiver for PN codes. Although a DLL could be used to receive the BINOR 
code, its performance is not easy to predict. The foregoing analysis assumes that a 
phase-locked loop, locked to the bit clock, is the principal instrument of phase 
measurement. This implementation causes the performance to be dependent upon 
code length, and results in a loss of about 2 dB in power requirement over a PN-FDMA 
code with equivalent parameters. 

Besides rapid acquisition (about 0. 2 s for the code described above), the BINOR 
code has the advantage in that the synchronizing and recovery circuitry, being prima- 
rily digital, is relatively inexpensive. This claim (Reference 1) has not been investi- 
gated in detail. 

It is concluded that for the accuracy and ambiguity requirements of the present 
problem the BINOR code achieves ranging accuracy only slightly inferior to that of 
PN-FDMA, The acquisition performance of BINOR is superior to PN, being about 0.2s 

at P /N = 42 dB/Hz. 

r o 

2.2.3 SIDETONES 

In a sidetone ranging scheme, a set of sinusoidal tones is transmitted. The 
frequency of the highest frequency tone is determined by the range accuracy required, 
assuming the same limitation on phase measurement resolution as in the case of 
binary codes. The lowest frequency is set by the ambiguity resolution required. The 
intermediate tones are used to resolve the ambiguity in each successively higher tone. 
The problem of optimum tone frequency selection has been discussed by several inves- 
tigators, e.g. , Reference 1 and Reference 6. A typical sidetone spectrum is shown 
in Figure 2-10. 

2. 2. 3.1 Phase Measurement 


The mean-square error in the time delay measurement made with a tone ranging 
system is determined by the upper tone frequency, f , the number of tones, and the 
received signal -to -noise ratio, and is given by m 
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where n is the number of tones. The assumption that each tone contains l/n of the 
total power is implicit. Normalized to the half-period of the upper tone, T^, this is 



( 2 - 11 ) 
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Normalized rms Error 




Figure 2-10, Typical Sidetone Spectrum 

This normalization permits a comparison with the corresponding normalized error for 
the binary code given by Equations 2-4 and 2-7. Equation 2-11 is plotted in Figure 2-11 
for n=5. If we assume a parallel between binary code clock rate and fine ranging tone 
frequency, then we can conclude (in a general sense) that performance is equivalent 
per unit bandwidth to a digital code for a tone system using 10 tones. For fewer tones, 
the performance is relatively better. 

2.2. 3.2 Data Transmission and Spectrum Compaction 

Unlike the spectra of the binary codes, which contain discrete components every 
l/N Hz where N is the ambiguity resolution period, the spectrum of a sidetone wave- 
form consists of a small number of components and a large amount of unused spectral 
capacity. The availability of the unused capacity can be increased by frequency- 
translating the component tones to more closely-spaced spectrum positions. A syste- 
matic way of doing this is, for a sequence of tones fj_, f£, . . . , f transmitting the 
tones ^ -f n _i, f n -f n _ 2 > •••» *n ~V 

For example, if a sequence 0.4 kHz, 2 kHz, and 10 kHz is required, the trans- 
mitted set would be 10 kHz, 9. 6 kHz, and 8 kHz. Provided that f -f ^ this 

n n-JL l’ 

procedure results in a compaction of the overall spectrum. The required tones are 
synthesized in the receiver from the received set. 

The reduction in accuracy that results from spectrum compaction is found as 
follows. In the receiver, the tone f^ is found as the difference of the two received 
tones f n -fj and f n . If the phase error variance in the recovery of each of these tones 
is equal to cr<^, and the two errors are assumed independent, then the rms phase 
error is equal to V2 The rms ranging power utilization is therefore reduced by 

3 dB . 


2-22 



Normalized rms Delay Error 



Figure 2-11. Ranging Accuracy Sidetones 
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Data transmission in conjunction with the sidetone waveform must be placed on a 
separate carrier and positioned in the spectrum between the tones and the main car- 
rier. The reduction of power available for ranging due to the transmission of data is 
-24 dBW/bit/s, as computed in the preceding subsection. The 10 bits/s assumed for 
this system requires negligible power. 

2.2. 3. 3 Fix Rate and Multiple-Access 

Sidetone ranging fix rate is determined solely by the receiver PLL bandwidth for 
the high frequency tone. Frequency-division multiple-access is required. 

2.2. 3.4 Summary 

The ranging accuracy of a sidetone ranging waveform is primarily determined by 
the amount of power in the highest frequency tone. Thus, proliferation of tones re- 
duces performance. In this respect, sidetone ranging is similar to BINOR code rang- 
ing which loses accuracy with an increase in the number of synthesizing square waves. 
A sidetone waveform, however, does not require a component tone at every harmonic 
of the ambiguity resolution frequency as does the BINOR code; thus, it is inherently 
more flexible and efficient. Its performance is also an improvement on PN ranging. 

Reception and demodulation of a sidetone waveform involves a set of PLLs, one 
locked to each tone frequency. This set follows a coherent demodulator and carrier 
PLL. 


Sidetone ranging achieves a somewhat better utilization of ranging power than 
digital codes for the problem parameters under consideration. 

2.2.4 PULSE RANGING TECHNIQUES 

Pulse ranging systems differ from the CW methods discussed in the preceding 
subsections in several important respects. Ranging accuracy is determined by the 
resolution of the time-of-arrival of a narrow pulse. Pulse receivers inherently oper- 
ate in an incoherent mode, using matched filter detection rather than serial corre- 
lators or phase-locked loops. Thus, signal acquisition is not a factor in pulse rang- 
ing, the low duty cycle of the pulse signals makes time-division multiple-access a nat- 
ual operating mode, and peak power devices are appropriate. 

The principal constraining factor in the generation of pulse ranging signals on- 
board a satellite is the liihitation on the peak power that can be generated by satellite 
transmitting equipment. The present state-of-the-art in high-power devices for space- 
borne applications is about 500 W (with some laboratory bulk effect devices producing 
an order of magnitude more peak power) at L-band using a Gunn effect negative- 
resistance oscillator. The device efficiency under these conditions is about 30%. 
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In order to mitigate the peak power problems, proposals for satellite-generated 
pulse waveforms generally assume the use of pulse compression techniques to obtain 
short, high-energy detected pulses by filtering the longer duration, lower power 
pulses which are actually transmitted by the satellite. Well-known methods for imple- 
menting pulse compression signaling are by the use of short PN codes and by the use 
of linear frequency modulation. The following analysis is general, and applies to both 
techniques. For comparison, an uncompressed pulse ranging approach is also ana- 
lyzed. Figure 2-12 shows a PN pulse spectrum and profile. 
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b. PN Pulse Time Profile 

Figure 2-12. PN Pulse Spectrum and Profile 

The functions to be performed by the pulse ranging waveform are identical to 
those performed by the CW waveforms — range measurements, ambiguity resolution, 
and data transmission. The range accuracy achieved is directly related to the com- 
pressed pulse bandwidth and received pulse energy to noise ratio as discussed below. 
Ambiguity resolution is determined by pulse repetition rate, which must therefore be 
less than 100 Hz (for 10 7 ft maximum range difference). Data transmission can be 
accomplished in a number of ways. Optimum error rates can be achieved by using 
antipodal binary codes’ with dual matched-filter reception. 
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2 . 2 . 4. 1 Range Measurement 


Assume the rise time of the compressed pulse at the output of the matched filter 
is T b- If the time-of-arrival is measured at the leading and trailing pulse edges and 
the two readings are used to obtain an average measurement, then (Reference 6) the 
delay error is given by 
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where E is the pulse energy, and B is the IF bandwidth. If B = l/r b> then 
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This can be related to peak P r /N Q by noting that 
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where T is the duration of the uncompressed pulse. In pulse ranging, different fix 
rates are obtained by averaging the range measurement over different numbers of 
received pulses. For averaging over n pulses, assuming independence, the error 
variance is reduced by l/n. Then Equation 2-13 is 



Equation 2-15 is plotted in Figure 2-13 for T = r b = 1.0 /as, that is, no pulse com- 
pression effect. It is observed, for example, that for cr- 10“ V s (100 ft ranging 
error), and n'=^l, the required peak P r /N 0 is 74 dB/Hz. The assumed link provides 
34 dB/Hz/W of satellite power; hence, 40 dBW, or 10 kW of peak power must be 
transmitted by the satellite* 
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Figure 2-13. Ranging Accuracy Pulse 
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In Figure 2-14 an uncompressed pulse length of 10*^ seconds and 100 pulses per 
second are assumed. The normalized range error, a’/ T ^ ) , is shown as a function of 
both peak and average transmitter power. The compression provides 20 dB reduction 
in required peak power; thus, the same performance as in the above example can be 
achieved with 100 W peak assuming as before al.0/is compressed pulse length. 

The use of pulse compression provides the ability to differentiate (by signal or 
code design) between signals from different satellites, perhaps eliminating the need 
for staggering transmission between satellites. 

2.2.4. 2 Data Transmission 

The transmission of satellite data is conceptually straightforward for pulse sig- 
nals. Two signal codes are used with dual matched-filter reception. In the case of 
chirp signals, the binary signaling states can be represented by an increasing or 
decreasing chirp frequency, while binary codes can be modulated by binary phase 
reversal. The ranging signals and data transmission signals are the same. 

In the matched-filter reception of pulse signals, attention must be paid to both the 
probability of a missed pulse and to the probability of a false alarm. False alarm 
rates can be kept very low by using predictive pulse gating based upon a priori knowl- 
edge of the pulse repetition frequency. 

.2.2.4. 3 Multiple-access Techniques 


The low duty cycle of the pulse signal from each satellite makes time division the 
most attractive and bandwidth-efficient multiple-access scheme. The user equipment 
consists of a single pulse receiver with the necessary timing and data decoding cir- 
cuitry. 

2. 2. 4. 4 Summary 

Other than PN-CDMA, pulse ranging is the only system that permits operation in 
a single frequency channel. On an average power basis, it achieves ranging accuracy 
comparable to that of a sidetone waveform (see Figure 2-10). The inherently TDMA 
nature of pulse transmission leads to somewhat simpler user equipment since receiver 
duplication for multiple signals is not required. 

The principal drawback J;o pul selranging is the requirement for the generation of 
relatively high peak powers on board the satellite. Present approaches to the design 
of spacecraft transmitters emphasize the use of constant power output devices, such 
as TWTs. As mentioned previously, devices presently under development will meet 
the power requirements for pulse systems using compression ratios of at least 100:1, 
but the peak power requirements for uncompressed pulses currently are beyond the 
state-of-the-art. 


'1 


2-28 



Normalized rms Error 



2-29 


2.2.5 MULTIPATH PERFORMANCE 

This subsection briefly discusses the problem of multipath reception and its effects 
upon the ranging waveforms. 

The reception by an aircraft of a satellite signal is susceptible to interference 
from the signal reflected upward from the ground plane. The severity of the interfer- 
ence is determined primarily by the strength of the spurious signal, its dispersion in 
time and in frequency, and the difference in propagation time over the two paths. 

The strength of the reflected signal, its time and frequency dispersion are in turn 
dependent upon the nature of the reflecting surface as described by its reflectivity coef- 
ficient and its average roughness, and upon the wave incident angle. The reflectivity 
coefficient, which is a function of incident angle, determines spurious signal strength, 
while the surface roughness determines the dispersion of the reflected signal in time 
and frequency. To a first-order-of-approximation the reflectivity coefficient is depend- 
ent only upon the incident angle and the dielectric constant of the surface. The total 
reflected power is independent of roughness. Experimental measurements of reflectivi- 
ty coefficients (Reference 8) for asphalt and short grass at L-band give results which 
for low elevation angles (10° to 40°) are generally less than -10 dB. At higher angles 
the aircraft antenna, being hemispherical in coverage pattern, can be expected to pro- 
vide additional discrimination against reflected signals. The prediction of quantitative 
performance of any given configuration of aircraft antenna, wave polarization, satellite 
elevation angle, and surface material can generally be done only by computer simula- 
tion or experimental measurements. 

The effects of surface roughness appear as a dispersion in the frequency of the 
reflected wave which, at L-band, is generally negligible, and a random dispersion in 
the path length traveled by the reflected wave which appears as a random time disper- 
sion in the interfering signal. At L-band, and assuming a surface roughness which is 
just sufficient to differentiate it from a smooth surface (’’critical roughness”) this 
time dispersion is on the order of 5 ^s, which cannot be considered negligible with 
respect to the baseband signals under consideration. In evaluating the ranging wave- 
forms for a particular system application the effects of the time dispersion of the 
multipath signal must be investigated by simulation and experiment. Note that time 
dispersion due to roughness is of the same order of magnitude as the time delay due to 
path length difference. 

Since the analytical prediction of multipath performance as a function of waveform 
design, on a quantitative basis, is relatively intractable, all that can be done is to 
evaluate the waveform candidates against certain qualitative design principles which 
lead to good multipath rejection. 

Ranging measurements are based upon the cross-correlation of the received wave- 
form with a stored replica, and it is intuitively evident that waveforms having autocor- 
relation functions (approximately an impulse function) are the most desirable. Blasbalg 
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(Reference 8) has demonstrated this result rigorously by showing that optimum multi- 
path rejection is achieved by signals having the flat spectrum of white noise. Accord- 
ing to this principle, multipath performance is improved by narrowing the central peak 
of the autocorrelation function and reducing the sidelobes. Reducing the width of the 
peak is achieved by raising the clock rate for binary signals, reducing the compressed 
pulse width for pulse coding, and raising the upper tone frequency for tone ranging 
signals. As a rough rule, the width of the correlation peak should be small compared 
to the path length differential delay. The reduction of sidelobes in the correlation 
function implies the removal of periodic components in the baseband spectrum (i. e. 
flattening the spectrum). Both of these criteria argue in favor of the use of noiselike 
(PN) signals and against the use of signals with strong specular components such as 
sidetones and the BINOR code. 

Multipath performance in a satellite system is strongly dependent upon the ele- 
vation angle of the satellite. At low elevations the differential delay becomes small 
while the reflectivity coefficient increases. Thus the spurious signal strength 
increases under conditions in which the increase of signal bandwidth is least effective. 
This double deterioration may argue in favor of a multisatellite constellation in which 
at least four satellites will always be at sufficiently high elevations that the multipath 
is acceptable. 

2.2.6 VE LOCIT Y MEASUREMENT 

Signaling waveforms which permit coherent carrier recovery can be used for 
direct measurements of the radial velocity of the user relative to the satellite. Radi- 
al velocity measurements from three satellites permit the calculation of user velocity 
in x, y, and z. The measurement is implemented by determining the frequency of the 
carrier received from the satellite and calculating radial velocity from the equation 

v - f A (2-16) 

r d 

where f^ is the Doppler frequency. Knowledge of the frequency transmitted is 
assumed, and X is the wavelength of the transmitted frequency. The mean-square 
velocity error resulting from thermal noise at the receiver input is given by 
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(2-17) 
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where cr is the phase-error variance at the output of the carrier recovery phase- 
locked loop and is the loop bandwidth. The phase-error variance is given by 
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where P r is the carrier power. Substituting (2-18) in (2-17) gives the result 
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Under the assumptions that X= 0. 65 ft , P T /]$ 0 = 40 dB/Hz, and the desired rms 
velocity error is 1 ft/s, the required loop bandwidth is 
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= 123 Hz 


It is noted that only those waveforms with a continuous carrier;, i. e. PN-FDMA, 
PN-CDMA, BINOR, and sidetones; can be used for velocity measurement. 

2.2.7 SEPARATION OF FUNCTIONS 

In the foregoing analyses of candidate ranging waveforms, we implicitly assumed 
that the satellite transmission consists of the repetitive transmission of a single rang- 
ing waveform containing no periodicity of a longer period than that required for ambi- 
guity resolution. Carrier synchronization, code (or frequency) acquisition, and range 
measurements are all carried out on the same waveform. Improved performance can, 
however, be obtained if the functions of ranging and ambiguity resolution are carried 
out on different waveforms at different times. 

The fundamental expression for delay resolution, from which the expressions for 
ranging accuracy in the preceding subsections are directly or indirectly derived, is 



( 2 - 20 ) 
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where /3 is the effective signal bandwidth defined by 
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where F(tu)cis the Fourier transform of the signal. A single sinusoid at frequency cu 
has, from Equation (2-21), an effective bandwidth cu. In ranging measurements, 
therefore, the optimum use of a baseband channel of width cu c is obtained by transmit- 
ting a single sinusoid of frequency cu c at full power. This is equivalent to a sidetone 
ranging waveform containing only a single tone. The lower frequency components, 
which are included in all the waveforms considered here, are required only for ambi- 
guity resolution; their effect on ranging accuracy is to reduce it by reducing the effec- 
tive bandwidth of the signal. 


This point was recognized in the design of the TRW NAVSAT navigation waveform 
(Reference 1). This waveform is based upon the use of the BINOR code for acquisition 
and ambiguity resolution while range measurements are performed during periodic 
transmission of a bit-clock signal at full power. The ranging accuracy performance, 
thus obtained, is that of a sidetone system consisting of a single tone (n = 1 in Figure 
2-11). It is clear, however, that both PN-FDMA codes and sidetone ranging wave- 
forms are susceptible to the same approach. When viewed in this way, it is seen that 
the essential differences in the waveforms consist of differences in the method used 
for ambiguity resolution. 


With functions thus separated, each waveform consists of a sequence of code (or 
multitone signal) of sufficient length to permit both initial acquisition and resynchroni- 
zation of the ambiguity-resolving loops (DLL for PN codes, and PLL for sidetones). 

This segment is followed by a segment of pure bit clock (or high-frequency tone) during 
which the precision range measurements are made. With this waveform structure, 
the approaches yield equal performance in both ranging and ambiguity resolution. Com- 
parative decisions between them must be made on the basis of considerations other 
than these — primarily the cost and complexity of user equipment. We note that nei- 
ther PN-CDMA nor pulse waveforms are viable candidates for separation of functions; 
hence, since all other candidates occupy the same bandwidth, bandwidth is not a factor 
in the comparison. 


In discussing the separation in functions in time we note that the data trans- 
mission function can also be separated in time from ranging and ambiguity determina- 
tion, and transmitted during a short time segment when full carrier power can be 
devoted to it. In this way a much higher data rate can be supported than in the case 
with continuous transmission at low power, and data can be sent with a low duty cycle. 
This approach may simplify the receiver compared to the user of a separate carrier. 
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2.2.8 DISCUSSION AND CONCLUSIONS 

In this subsection, the normalized delay error curves are used as the basis for 
comparison of performance from an operational point of view. For convenience, these 
curves are reproduced in Figure 2-15. In addition to the performance curves for the 
five waveforms considered in detail, two additional curves are shown. One of these 
shows the ranging accuracy that is achieved by separating the waveform functions as 
discussed in the preceding subsection. It is equivalent to the performance that would 
be obtained by a single sinusoid at the highest baseband frequency. The second addi- 
tional curve gives the performance of a pulse ranging system operating with no com- 
pression and with no averaging over groups of pulses. Although this curve is plotted 
with the same normalized ordinate as the other curves, it holds only for a pulse length 
of KP^second. Separate labeling of the abscissa shows the peak power that must be 
transmitted from the satellite for each of the pulse waveforms. The average power 
abscissa holds for all the curves. 

It is seen that sidetones, PN-FDMA, and BINOR waveforms achieve roughly 
equivalent performance, with the differences between them arising from differences 
in effective bandwidth. Of the three, sidetones are the most efficient in the utilization 
of ranging power. This is an intuitively satisfying conclusion, since in a sidetone 
system power is allotted only at those specific points in the spectrum where it is 
required, whereas in the binary code waveforms power is diffused over the entire base- 
band. Although this feature contributes to the ranging efficiency of sidetone waveforms, 
it is, as discussed below, an disadvantage in other respects. 

The performance of the 100:1 pulse compression waveform is essentially equal 
to that of the five-tone sidetone waveform for equal average power. A 3-dB factor 
in its performance arises from the assumption that two delay measurements are made 
on each pulse. If only one measurement is made then the pulse performance curve 
becomes equal to that of PN-FDMA. 

PN-CDMA suffers in performance both at high and low values of transmitted 
power relative to the FDMA and TDMA waveforms. The deterioration in performance 
of this approach must be traded against the reduction in required RF bandwidth. 

The curves of Figure 2-15 can be viewed from several different perspectives. 

For example, if the rms error is fixed at the same value for each waveform the 
curves become those of required RF bandwidth (per satellite for FDMA and TDMA, 
and total for CDMA) as a function of transmitted power. Similarly, for a fixed band- 
width, the curves show rms error performance as a function of transmitted power. 

To determine the overall bandwidth, efficiency of each waveform the total band- 
width required by the particular multiple-access scheme must be taken into account. 

Of the waveforms analyzed, only in PN-CDMA and pulse signaling can the signals from 
all satellites be accommodated in a single RF channel. The remaining waveforms 
require FDMA and use one channel per satellite. This factor must be taken into 
account in using the curves of Figure 2-15 to determine relative bandwidth efficiency 
of the different waveforms. 
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Figure 2-15. Comparative Ranging Performance 


It is reasonable to conclude on the basis of Figure 2-15 that the ranging per- 
formance and bandwidth occupancy of sidetones PN-FDMA and BINOR are sufficiently 
close that a decision between them should be made on the basis of other considerations, 
such as user set complexity and susceptibility to multipath interference. 

The fix rates attainable by the candidate waveforms are determined by the 
receiver loop noise equivalent bandwidth (or averaging time for pulse waveforms). 

Since this parameter can be independently chosen it is concluded that attainable fix 
rate is not a significant trade-off factor in selecting between the approaches consid- 
ered. 


In contrast to the conclusions relating to the efficiency of utilization of power 
for ranging purposes, the discussion of subsection 2.2.5 points to the conclusion that 
waveforms with relatively flat spectra, such as binary, are less susceptible to multi- 
path interference than are waveforms with strong specular components. The difficulty 
of obtaining quantitative estimates of this effect prevents a trade-off being made at the 
level of this study, except to note the importance of investigating this area, particu- 
larly since the ranking of waveforms relative to multipath susceptibility appear to be 
opposite to the ranking for ranging performance. 

In considering the user set complexity required to implement receiving and pro- 
cessing circuitry for the waveforms considered, several observations can be made. 
The waveforms that require separate RF channels (PN-FDMA, sidetones, and BINOR 
code) require a wideband RF receiver and a separate demodulator and baseband pro- 
cessor for each satellite signal. The use of PN-CDMA places all the signals in a 
single RF channel, but separate processors are Still required for demodulation and 
processing for each satellite. Comparison of these approaches requires further 
design of the user set to permit a cost comparison to be made at a lower level of 
design detail. The only approach which seems to offer a clear advantage in user 
costs is the pulse system, in which all the satellite signals time-share a single 
channel and a single receiver. 

As a result of the investigation reported herein, it is concluded that none of the 
waveforms considered offers a clear and uniform advantage over the others relative 
to the comparison criteria taken into account. Selection of a waveform for a particu- 
lar system application must be done on the basis of evaluation against what may be 
called secondary criteria, many of which will be system-specific. The selection may 
be made on the basis of these secondary criteria with reasonable confidence that the 
system performance in ranging and position location will not be severely degraded 
regardless of the selection. ^ 
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2.3 SEQUENTIAL POLLING 


This Subsection discusses the mechanism of polling and response which lies at the 
heart of elicited-response surveillance systems. The discussion is concerned with 
specific polling procedures and the design of the interrogation waveform. 

Elicited-response systems depend upon the use of interrogation waveforms that 
key user responses, at appropriate times. The interrogation waveforms are trans- 
mitted on the satellite downlinks to groups of users who respond to the recognition of 
a specific point in the waveform. The basic problem in a sequential polling waveform 
design is to assure that the responses in a given frequency channel do not overlap at 
the repeaters of the relay satellites unless protection is provided by the use of code 
division multiple access. Stated another way, the problem is to sequence responses 
from groups of users using the same response waveform in the most efficient way. 

The questions to be answered are 

1) In what order and with what time spacing must the users be interrogated? 

2) What is the best signal design for the interrogation waveform? 

The answer to the first question depends upon whether a system with a single response 
path (Configuration I) or a system with multiple paths (Configuration III) is being con- 
sidered. In Configuration I, two approaches to range-ordered polling can be consid- 
ered. For Configuration III, the requirement to avoid signal overlaps simultaneously 
at several satellites imposes additional constraints which dictate a different approach 
to the problem. 

2. 3. 1 RANGE-ORDERED POLLING 


Consider the situation shown in Figure 2-16. A group of users is to be polled by 
an interrogation waveform such that their response messages, each of length r m , are 
sequenced through the repeater of the relay satellite. The number of users that can 
be polled in a time T u is clearly maximized by minimizing the guard time between mes- 
sages. If the control center has knowledge of the positions of all users, this can be 
done by placing the interrogation point for user j in the interrogation waveform at a time 
Ty where Tj is given by 
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( 2 - 22 ) 


in which r m is the message length, Arj is the range difference between user j and user 
j + 1, and c is the signal propagation speed. The situation is shown in Figure 2-17 in 
which the horizontal axis represents time and the vertical axis represents distance be- 
tween the satellite repeater and the users. A single frame of polling and response is 
shown, although in practice the process is continuous. Equation 2-22 can be under- 
stood by observing that if all the users are at the same range the second term is zero 


2-37 





and the interrogation points are spaced by the message length. The second term cor- 
rects for the propagation delay resulting from unequal ranges. It is noted that the or- 
der of the interrogation points is not the same as the order of the responses, and also 
that more than one user may have to respond at the same (or almost the same) point. 
This rolling method is referred to as "minimum guard time" range-ordered polling. 


The implementation of this minimum guard time approach introduces several dif- 
ficulties. First, the processor at the control center must compute the propagation de- 
lays for each user at each interrogation scan, and order the interrogation waveform to 
assure the correct sequential response. Second, the ordering of the interrogation 
waveform must be changed every scan to reflect changes in user positions regardless 
of whether the ordering by range has changed. 

This approach maximizes the number of users that can be polled sequentially. The 
number of users, N, that can be polled in a given time T u is limited only by the mes- 
sage length, and is given by 


T 

N = (2-23) 

m 


The implementation difficulties encountered in minimum guard time polling can be 
eased considerably at a slight sacrifice in capacity. Assume that the interrogations are 
ordered by user range, and are spaced at the transmitter at intervals equal to T m . 

Then for users at the same range the responses will be sequenced through the repeater 
with essentially zero guard time, as in the previous approach. The responses from 
users whose ranges differ by an amount Ar^ will be spaced at the satellite by a time in- 
terval equal to 2Ar/c. If rj is the range difference between the i’th user and the 
(i + l) f th user, then the total time, T u , required for the responses from N users is 
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If R is the range difference between the closest and furthest users, then* 


N 

X Ar. - AR (2-25) 

A 1 

and 

T _ 2_AR + Nt (2-26) 

c m 


2-39 



For specified T, Equation 2-26 limits the number of users that can be polled using 
this technique. N is given by 


2 AR 


N = 


m 


(2-27) 


This simplified range-ordered polling scheme represents a decrease of 2 AR/c r m in 
the number of users from the minimum guard time approach. An estimate of the mag- 
nitude of this reduction can be made using the typical figures of T u = 1 second, r m = 
1(T 5 second, and AR= 3000 miles (CONUS). Then, 

2 AR 2 x 3000 

c t 5 -5 

m 1.86 x 10 x 10 

= 3220 

which represents a loss of some 3% of the 10 5 users which could otherwise be polled. 


2. 3. 2 SEQUENTIAL POLLING FOR MULTIPLE SATELLITES 


The conceptual simplicity of the range-ordered polling method depends upon the 
assumption that the user signals return to the control center by way of a single satel- 
lite, and hence, this method of polling is confined to use only with a Configuration I 
system. Assuming only time orthogonality, in the multi satellite elicited polling associ 
ated with Configuration in, the problem is complicated by the fact that the response 
messages must be sequenced in such a way that no overlap occurs at any of the relay 
satellites. Clearly, if a set of users are range-ordered with respect to one satellite, 
the same range-ordering cannot hold for any other satellite. Another approach to the 
problem must be taken. 


Assume that a group of N users is distributed randomly over the total coverage 
area of interest, with the maximum possible distance between any pair being AR, and 
consider a low-elevation satellite as shown in Figure 2-18. Assume also that the in- 
terrogation signal is being transmitted by the same satellite. If User 1 and User 2 are 
interrogated successively, then the response messages will not overlap if the interro- 
gations are spaced by an interval t given by 
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Figure 2-18. Low-Elevation Satellite 

If no overlap is assured for the lowest elevation satellite, it is also assured for satel- 
lites at higher elevations. For a specified polling interval of T u seconds, the number 
of users that can be polled is 


T 

u 

~ 2 Alf^ 

c T m (2-29) 

A generalization of this idea is obtained by dividing the total coverage area. Ax, 
into n equal-area square cells of dimension A^/2, where A is the cell area. Within 
each cell, users may be interrogated in random order, with no overlap of messages 
provided that the successive interrogations are spaced by an interval given by Equation 
2-28. In this case, we have 



AR = (2A) 
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since the maximum possible separation of users within a cell is the length of the cell 
diagonal. In polling all users, the procedure is to interrogate all the users in a given 
cell, with spacing given by Equation 2-28 before passing to the next adjacent cell, ac- 
cording to a prespecified and fixed cell order. In determining the number of users that 
can be interrogated using this procedure. Equation 2-29 must be modified to reflect the 
fact that when the interrogation sequence moves from one cell to another, the maximum 
distance between the last user polled in one cell to the first user polled in the next cell 
in greater than AR. For simplicity, we ignore the fact that the true distance is the 
diagonal across two cells, and use 2 x (2A) 1 / 2 as an approximation. The interrogation 
interval in this case is 
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which exceeds the value given by Equation 2-28 by 2(2A)l/2/c. This additional amount 
must be allowed for each of the n cells. In a total polling interval of time T U . there 
will be N polling intervals of length r, and n cell-transition intervals of length 
2(2A)l/ 2 /c. Summing these intervals and equating to T u gives 
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Solving for N, 
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The question to be investigated is how the cell area. A, should be chosen to maximize 
N for a given Tu* T m» and A^. Noting that 


n = 


A ’ 


Equation 2-33 can be rewritten 
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Differentiating Equation 2-34 with respect to A 1 / 2 and equating the result to zero leads 
to 
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The solution for the positive real root of Equation 2-35 is 
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It is easy to verify that A 1 / 2 given by this expression corresponds to a maximum of N, 
which we denote by N m . Table 2-1 gives a tabulation of values of N m and A 1 / 2 for 
various values of At, T u , and r m* The number of cells, n, is also given. 
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Table 2-1 


1/2 

N , n, AND A 7 FOR VARIOUS VALUES OF A , T , AND T 


a t 

2 

(mil ) 

T 

u 

(s) 

T m 

(A^s) 

N 

m 

n 

A 1/2 

(mi) 

4 x 10 6 

10 

1 

26,742 

26,963 

12.2 



10. 

24,670 

25,682 

12.5 



100 

12,670 

18, 310 

14.8 


3 

1 

2,431 

2,439 

40.0 



10 

2,396 

2,427 

40.0 



100 - 

2,324 

2,391 

41.0 


1 

1 

270 

271 

12.2 



10 

269 

270 

12.2 



100 

257 

257 

12.5 

4 x 10 5 

1 

1 

2,674 

2,696 

12.2 



10 

2,440 

2,568 

12.5 



100 

1,267 

1,831 

14.8 


The value of A t of 4 x 10^ square miles corresponds to a CONUS coverage sys- 
tem, At of 4 x 10 5 square miles corresponds to a system in which the CONUS area is 
divided into ten sectors, within each of which the interrogation procedure is independ- 
ent of the others. The table entries for N m and n for At = 10 5 square miles are per 
sector. 

The number of users that can be interrogated using this procedure is quite low, 
particularly for small values of T. This is not surprising when it is noted that the 
channel duty cycle is low. That is, the message length is generally short compared 
to the interval between successive response messages. In designing a practical system, 
the capacity that is unused under a strict TMDA procedure is filled by users separated 
by code division. Various ways of doing this are analyzed later in Subsection 2. 5. 

As long as a user remains in the same cell, he can respond at the same point in 
the interrogation waveform. When he moves from cell to cell, his response point in 
the interrogation sequence must be changed. The frequency with which this occurs is 
important, as it affects the processing load imposed on the control center and also the 
data rate on the link from control center to user. The smallest cell dimension in Table 
2-1 corresponds to an update rate per user of about 3 minutes at an average speed of 
250 mi/h. The required downlink data rate is investigated in the next subsection. 
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It should be noted that the determination of N m makes no assumptions about the 
distribution of users, and the analysis is valid for any distribution, except the unlikely 
case in which all the users occupy the same cell. In this case, N m is larger, but it is 
obviously not of practical interest. 

2. 3. 3 DESIGN OF THE POLLING WAVEFORM 

The polling waveform must perform several functions. It must 

1) In two-way ranging measurements (Configuration HI), be designed with 
sufficient fine structure to permit the user to recognize and respond to 
its assigned interrogation point with a time accuracy determined by the 
required position error 

2) Transmit data for changing interrogation points when necessary 

3) Transmit data to selected users for collision avoidance 

4) Operate within the downlink power restrictions imposed by the low user 
antenna gain and satellite prime power limitations. 

The requirement of Item 1 is relaxed somewhat in Configuration I, in which the inter- 
rogation signal is not itself part of a ranging measurement. 

First consider the requirements of interrogation: there are two basic approaches 
that can be taken. First, the users can be polled by transmission of user addresses. 
When a user recognizes his address, he responds with the appropriate message. This 
approach has the advantage that changes in the polling sequence and can be accomplished 
simply by rearrangement of the addresses by the control center. No specific data trans- 
mission to the user is necessary to change response points. A second possibility is to 
transmit a common timing waveform which is received and monitored by each user. 

The waveform repeats every polling scan. Each user is assigned a point in the timing 
waveform to respond. An obvious candidate for the timing waveform is a long PN se- 
quence. Using this approach, the polling sequence must be changed by signaling each 
affected user to change his response point. The basic criterion which is used to evalu- 
ate the merits of these two possibilities is the data rate required on the satellite to 
user downlink. A secondary criterion is the processing load imposed on the control 
center to generate and update the polling sequence. 

2. 3. 3.1 Configuration I 

In Configuration I, it is easy to show that minimum guard time range-ordered poll- 
ing imposes an excessive bit rate requirement on the downlink. The sequencing re- 
quirement necessitates a waveform with a fine "granularity" of response points which 
can only be s ati sf i ed by a common timing waveform/ But a common timing waveform 
requires explicit update messages when the user response points change, which must 
be done every scan. Assume N users, and B bits per update message, and a scan time 
of T seconds. The downlink data rate R is 



For N = 10 5 , B = 10 2 bits, and T = 1 second, R = 10 7 bits/s. The postulated satellite 
link, with a characteristic of 36 dB/Hz/W, would require about 50 kW of transmitted 
satellite power to achieve a bit rate with a 13 dB E/N 0 . The minimum guard time 
approach can be dismissed on this basis alone. 

The use of a common timing waveform with the simplified range-ordered polling 
method can also be dismissed by the same argument, since the range -ordering will 
change from scan to scan. Consider the implementation of this method with discrete- 
address polling. The polling addresses must be transmitted with a spacing essentially 
equal to the response message length. Assuming M bits for each address and a re- 
sponse message length of r m seconds, the downlink data rate is 

R = — (2-38) 

r m 

which is independent both of the total number of users and the scan time. Equation 
2-38 is plotted in Figure 2-19 with M assumed to be 20 bits. It is noted that the rates 
are sufficiently high that the downlink power may be the controlling factor in establish- 
ing the polling rate using this method. 



Figure 2-19. Downlink Data Rate Discrete-Address Polling 
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2.3. 3.2 Configuration IE 

In the multisatellite elicited-response system, the rate at which users can be in- 
terrogated is limited by the geometric considerations discussed previously. The dis- 
crete-address data rate can be found directly from the scan time T u , and the number 
of users, M: 



u 


(2-39) 


For example* for N = 105, M = 20 bits, and T u = 1 second, R = 2 Mb/s. This rate re- 
quires 10 kW satellite power, and appears to be unworkable for this reason. 


To find the bit rate for the common timing waveform, the rate at which update 
messages must be sent must be found. This is equal to the rate at which users move 
between cells. Assuming an aircraft speed v and a cell dimension A*/2 # the average 
number of cell transitions per second, M, with uniform user density is given by 


M = 


Nv 

. 1/2 


(2-40) 


The data rate for M-bit update messages is then 
MNv 


R = 


A 1 ^ 2 x 3600 


bits /second 


(2-41) 


For M = 150, N =10^, v = 250 mi/h, and Al/2 = 25 miles, R - 41. 5 kb/s. The re- 
quired satellite power for this data rate is about 200 W. 

2. 2. 3. 3 Waveform Structure 

When a user moves from cell to cell, his interrogation point in the polling wave- 
form is moved to a point associated with the cell he is entering. This implies that 
such a point is available (that is, unoccupied), otherwise the change of the one interro- 
gation would require the change of many others. An approach that accomplishes this is 
to preassign segments of the waveform to each cell, with the duration of the segment 
determined by the maximum number of users expected in that cell. This is shown in 
Figure 2-20. The success of this approach depends upon providing sufficient capacity 
that the probability of a cell becoming saturated is suitably low. 
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Figure 2-20. Interrogation Waveform Structure 
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2.3.4 SUMMARY 

In this subsection we have discussed some of the aspects of the polling problem for 
elicited-response systems. Sequential polling methods for both single-satellite and 
multiple-satellite response systems were considered, and some of the implications of 
the polling procedure on the downlink interrogation signal design were investigated. 

. In the single-satellite (Configuration I) system, two approaches to range-ordered 
polling were considered. The first, which involves sequencing the user responses 
through the satellite repeater with essentially zero guard time, was found to require 
excessive downlink power in the interrogation signal, in addition to being undesirable 
because of the load it imposes on the ground processor. The second simplified ap- 
proach is feasible with discrete-address interrogation, and involves only a small loss 
in capacity. 

For multiple- satellite (Configuration IE) signaling range-ordered polling cannot 
be used, and a range-cell approach to polling was developed. It was shown that the cell 
arrangement can be chosen to maximize the number of users interrogated. In this case 
the common timing waveform is the preferred interrogation signal design. It is noted 
that the ground processing load imposed in Configuration I range-ordering can be re- 
duced by adopting a range-cell approach at some sacrifice in capacity. 

2.4 CONFIGURATION I ANALYSIS 

The objective of the following paragraphs is to describe the aircraft-to-satellite 
performance of pairs-of-polling and multiple-access schemes in terms of the per- 
formance parameters of Subsection 2. 1.3. The determination of polling/multiple ac- 
cess pairs is the most significant design decision to be made in selecting preferred 
ATC surveillance communications techniques. A uniform analytical procedure will be 
developed so that alternatives can be viewed and interpreted on an equal basis. Subse- 
quent subsections, will apply the result to Configurations II, IE, and IV. The ground- 
torsatellite links are not considered to involve practical technical limitations since 
large aperture ground antennas can be used. 

In Configuration I, the user derives his position information from a navigational 
system such as 

1) IMU 

2) Omega/Loran 

3) Navigational Satellites. 
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These navigational systems can be assumed to be independent of the ATC system. 
It is assumed that regardless of the navigational data source, the user makes suitable 
measurements on the received navigation signals, and the resultant data is formatted 
for transmission to the ground center. Investigation of Configuration I, then, revolves 
about purely communications link considerations, since ranging and position deter- 
mination is independent. For this analysis a performance matrix for each polling 
technique and multiple access technique was prepared. The analysis leading to the 
contents of these matrixes is presented in Appendix A for each multiple access tech- 
nique. 

2.4.1 ANALYSIS PROCEDURE 

To illustrate this unified analysis procedure the data sheet for elicited TDM A/ 
CDMA is presented in Table 2-2. Essentially the top of the data sheet describes the 
multiple access technique and applies realistic signal-to-noise factors. Subsequently 
the expressions for the prime performance parameters are presented and calculated 
as a function of number of system users and transitional points between access tech- 
niques. From the results of these analyses a normalized bandwidth factor was de- 
rived, calculated, and plotted for each polling-multiple access technique. These 
graphs illustrate the interrelated trends of key resource parameters of bandwidth, 
user transmitter power, and code process gain as a function of the number of system 
users. The resultant plot is shown in Figure 2-21 for elicited TDMA/CDMA. Similar 
graphical presentations for each technique are contained in Appendix A. 

Table 2-2 and Figure 2-21 will be described in some detail illustrating the content 
and use of the additional figures supplied in Appendix A. 

The description, elicited TDM/CDM' implies that user aircraft transmit only upon 
interrogation, and that responses are caused to fall into particular time slots relative 
to the satellite as illustrated in Figure 2-22. Furthermore, more than one user can 
occupy a particular time slot, via CDMA, and the number of users per slot, N s , is 
controlled (by assignment) and is given by the total number of system users divided by 
the number of time slots, N/R^. The first entry in Table 2-2 is the signal or code 
division process gain which is required to accommodate N total system users with Rt 
time slots (Rt varies from 10 to 1000). The process gain (PG) requirement was taken 
to be that PG which would result in signal -to-clutter (desired signal-to-other signals) 
ratio corresponding to a BER* of 10-5, The second entry, Tj, is transmit time per 
interrogation, during which M bits must be transmitted. This results in a real time 
data rate of M/Tj and a message bandwidth of Bj - D/Hz (assuming a 1-bit/Hz modula- 
tion efficiency), for the third and fourth table entries. Since this example is elicited, 
or ordered, the probability of a clear^channel does not apply and is entered as unity. 
The sixth entry, total bandwidth Bt> reflects the process gain requirements and the 
message bandwidth and is given as the product of the desired signal-to-clutter ratio, 
the maximum number of simultaneous users per time slot, and the detection bandwidth. 


*Bit error rate. 
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PARAMETRIC ANALYSIS OF ELICITED TDMA/CDM A 
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Normalized Bandwidth, -rjr T , for T = 1 s, and M = 100 bits 



1 10 10 2 10 3 10 4 10 5 10 6 

Number of Users (N) 


Figure 2-21. Normalized Bandwidth vs Number of Users 
for Elicited TDMA/CDMA 
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Figure 2-22. Time Profile for Elicited Polling, 
TDMA/CDMA Multiple Access Pair 


The seventh entry, required transmit power, is calculated on the basis that a 20-dB 
signal-to-thermal noise is necessary, where CDMA is used, so that the BER is prin- 
cipally a function of the signal-to-clutter ratio. The eighth entry is the total number of 
system users per megahertz, an indication of the relative spectrum utilization of each 
candidate. 


The next entry is the timing accuracy required by both the interrogation signal and 
the user response if less than one data bit overlap is to occur in the time division 
format. It is interesting to note that approaches which use elicited TDMA can impose 
position location accuracy requirements on the navigation system. For example, in 
Table 2-2 with 1000 time slots, the position of the aircraft must be known at the time 
of transmission to within 10^ ft in order to ensure proper arrival time at the satellite. 
If the ground control station is controlling the transmit time of the user response, the 
impact of the delay between transmission and the aircraft velocity on this requirement 
is evident. 


Finally, Table 2-2 illustrates the minimum number of parallel ground signal proc- 
cessors needed to accommodate the requirements of N users. Here it is evident that 
the use of TDMA reduces the number of processors from N to N divided by the number 
of time slots. 


Figure 2-21 summarizes the parametric variation of total (normalized) bandwidth 
occupancy, transmitter power, and processing gain as a function of the total number of 
system users. The figure shows the tradeoff between transmitter power and signal 
processing gain. For example with N=10 5 and 10 time slots, 10 W of transmitter 
power and 53 dB of signal processing gain is required. From a conservative point of 
view, this signal processing gain probably represents a borderline state-of-art speci- 
fication. With N=lo5 and 100 time slots the transmitter power is 750 W and the proc- 
ess gain is 43 dB, perhaps the best compromise for elicited TDMA/CDMA combina- 
tions since the timing requirements to support this degree of TDMA are minimal. In 
any case, for large N, the bandwidth occupancy can be seen to be independent of the 
mix of TDMA and CDMA. The real value of these curves can best be derived by over- 
lapping or comparison with other alternatives such as autonomous TDMA/CDMA, etc. 
Operating points for the popular proposed ATC systems can be located and tradeoff for 
parameter variation can be interpreted. 
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Similar analytical procedures and graphical presentations were made for a polling- 
multiple access pairs and are presented in the appendices. 

2.4.2 SYSTEM FIGURE OF MERIT 


The analysis in the foregoing subsection can be further generalized by eliminating 
the link assumptions and the assumption relative to modulation efficiency, and guard 
frequency and guard time allocations. As indicated above the principal resources are 
transmitter power, total occupied bandwidth, and signal processing gain. A preferred 
system approach would minimize the product of these three items, i. e. , Txp * B^PG. 
For convenience, a figure-of-merit describing this performance was defined as 


Q 


T 


L. G G 
fs ac sc 


KT 


* PG 
T 


where 


L p G G 
fs ac sc 

KT 


is the product of free space loss, aircraft antenna gain, and spacecraft antenna gain di- 
vided by KT, which is the Johnson noise in a 1 -Hz bandwidth. Table 2-3 illustrates the 
form which this product takes for polling /multiple access pairs. The entries are gen- 
eral enough to facilitate the comparison of any candidate approaches. . As an example, 
let us compare the figure-of-merit for autonomous FDMA and TDMA: 
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where 


To = time interval per bit 

Ci = constant for modulation efficiency 

Cp = constant for FDM guard channel space. 

Here, 
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is the detected signal-to-noise ratio per data bit, a design factor to be specified and 
one which affects both candidates equally. We observe that the 
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factor; i. e. , the number of bits per message, M, times the modulation efficiency (bandwidth 
in Hz divided by the number of bits per second) C, divided by the interval between messages, 
Tu; appears squared for both candidates. Consequently, the Q factor is more sensitive 
to changes in these design factors than to SNR. Finally it is most interesting to note 
that Q for autonomous TDMA is a function of 10^ and N squared. The 10^ factor 
results from a requirement for the probability of a clear transmission of 1-10 -5 . For 
105 users this means that on the average two users are overlapped all. the time and 
the ground processor must continually input data from two users. * In any event the 
A- TDMA scheme increases Q as the square of the number of users whereas the Q for 
A-FDMA varies only linearly with N and Cp, the frequency guard band. 

The following additional observations can be made from Table 2-3: 

1) Q varies as N 2 for all candidates excepting those which employ FDMA, 
the FDMA hybrids vary as N 2 /Rp where Rp is the number of FDMA 
slots. For pure FDMA, Q varies linearly with N. 

2) All candidates with CDMA include a factor (Pt^b/^c) which is the de- 
sired signal-to-clutter ratio (clutter is effective noise power due to all 
other code division signals) 

3) Autonomous TDMA without FDMA involves a factor approximately equal 
to the reciprocal of the probability of two transmissions overlapping. 

For large N, and the condition that only two signals overlap, this factor 
is equivalent to one additional N 2 multiplier. 

4) For the autonomous TDMA/CDMA hybrid if enough process gain is pro- 
vided to allow for three times the expected number of signal overlaps 
than the additional N 2 factor associated with autonomous TDMA is re- 
duced to nine times the ordinary CDMA (Pr T D/ N e) factor. 

The net conclusion from the above analytical procedures is that if the performance 
of satellite ATC surveillance systems are to be compared only on user ATC trans- 
ponder power, total system bandwidth, and processing gain then candidates involving 
FDMA are desirable. We observe that for a ’’tactically" oriented system; i. e. , those 
with many divergent low-power, low-antenna-gain users; that this conclusion is con- 
sistent with minimum user costs and the traditional command and control approaches. 

It is also observed that in trunkline, point to point systems, wherein the system meas- 
ure of performance is strictly capacity per megahertz, that an entirely different basis 
for comparison and selection of techniques is appropriate. Considering that transmit- 
ter power, total system bandwidth, and signal processing gain are limited resources 
constrained either by availability o r sta te-of-art, the above analysis offers a useful 
evaluation tool for comparing the many possible ATC surveillance schemes. The fore- 
going evaluation procedure can be used for all possible ATC satellite configurations 
and therefore provides a uniform basis for comparing system approaches which can be 
radically different. 


*Systems with smaller numbers of users (N) could conceivably tolerate overlap rates 
greater than 10~ 5 in which case the multiplicative penalty for A- TDMA would be re- 
duced. 
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Figure 2-23 graphically represents the results for the standard link assumption, 
fix intervals, Ci = C^p = Cp = 1, etc. , that have been assumed for this study. A de- 
tailed analysis of all of the candidates has been completed and documented in Appendix 
B. Based on these results, a preliminary elimination of Configuration I candidates can 
be made on the basis that an eventual growth to 10^ simultaneous users is necessary. 
The results of this elimination procedure and the rationale is illustrated in Table 2-4. 
The effects of space division multiplex, i. e. , the use of multiple antenna beams will 
be treated subsequently. 

2.4. 3 SPACE DIVISION IMPLEMENTATION IMPACT ON TECHNIQUE ANALYSIS 

The foregoing multiple-access/polling techniques can be implemented in conjunc- 
tion with space division satellite techniques. A space division implementation is feas- 
ible by configuring 10 to 20 antennas on a spacecraft. If each antenna pattern is de- 
signed to cover a subdivision of the COLM area such that a near uniform traffic density 
exists for each antenna beam, critical code, power, or bandwidth requirements of the 
foregoing techniques can be reduced. In accomplishing this implementation several 
operational problems exist. One problem is range ordering in an elicited TDMA tech- 
nique. Another problem is frequency, code, or time slot management between adjacent 
and overlapping beam patterns. The ground processor may have to have flight plans 
stored for each active aircraft, and direct aircraft to change frequency or code as they 
move from one beam zone to another. These facets of space division are not consid- 
ered in depth. Rather, the analysis is constrained to deriving the impact of a 10 beam 
space division system on each of the foregoing techniques. For this analysis, it is as- 
sumed that 105 users are uniformly distributed between the 10 beams. 

Table 2-5 summarizes the effect of using space division on most of the previously 
described candidates. The transmitter power reduction reflects not only reduced N 
factors but also an increase in spacecraft antenna gain of 10 dB. Both E-TDMA and 
A&E- FDMA/CDMA are subsequently added to the list of potential candidates. 

2.4.4 SUMMARY AND CONCLUSIONS 

As derived from figure-of-merit calculations, the two candidates retained for 
further trade-off analysis against other configurations are elicited TDMA/FDMA and 
autonomous or elicited FDMA. Both of these candidates provide the least user com- 
plexity in terms of power and implementation while utilizing bandwidth effectively. 
Although the TDMA/FDMA approach requires the ground control to user link for all 
surveillance replies, a reduction in system timing and channelization accuracy re- 
quirements relative to pure TDMA occurs. 
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Figure of Merit, Q = 


Autonomous TDMA 



Number of Users, N 


Figure 2-23. Comparison of Polling Multiple Access Technique 


2-56 



Table 2-4 


SUMMARY OF CANDIDATE ELIMINATION 


Technique 

Rationale for Elimination 

Autonomous TDMA 

Excessive Bandwidth (lO^ Hz) 

Excessive User Transmitter Power (5 x 10® W) 

Autonomous TDMA/ 
FDMA 

Excessive Bandwidth (1012 Hz) 

Elicited TDMA 

Requires 10 foot Position Location to Implement 
Requires Uneconomical User Transmit Power 
(5 x 10 4 W) 

Autonomous and 
Elicited CDMA 

Excessive Processing Gain (63 dB) 
Requires 105 Codes 

Elicited FDMA 

No Advantage over Autonomous FDMA Unless 
Ground-User Link Used for Advisory Messages 

Autonomous and 
Elicited FDMA/ 
CDMA 

Excessive Bandwidth 

Elicited TDMA/ 
CDMA 

Requires Ground Control to User Link, and Of- 
fers no Bandwidth Advantage 

Autonomous TDMA/ 
CDMA/FDMA 

Excessive Bandwidth 

Autonomous TDMA/ 
CDMA 

Excessive Bandwidth (6 x 10® Hz) 

Technique 

Rationale for Candidate Retention 

Elicited TDMA/ 
FDMA 

Bandwidth = 10 7 Hz T^ = 50 W 
Reasonable Implementation for Rt = 100 
Effective Bandwidth Utilization, but Requires 
Ground-User Link 

Autonomous FDMA 

Bandwidth = 107 Hz Txp = 50 W 
Best Figure-of- Merit, no Ground-User Link 
Continuous Surveillance, but Requires 10~8 
Frequency Stability and in Addition Requires a 
Methodology of Eliminating Doppler Shift 
Ground Station Must Process 10 5 Separate 
Signals 
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Table 2-5 

SUMMARY OF APPROACHES INCORPORATING SDMA 


Approach for N = 10® 

PG 


Bj 

Comments 

Eliminated Candidates 

! • — -w ‘ : 

A-TDMA/SDMA 

1 

5 x 10® 

loll 

Power and bandwidth excessive 

E-TDMA/SDMA 

1 

5 x 10 2 

10® 

Power marginally high 

A&E CDMA/SDMA 

53 dB 

0.25 

2 x 10 7 

Process gain marginally high, 
bandwidth relatively high 

A-TDMA/CDMA/ 
SDMA 
R t * 100 

38 dB 

25 W 

6 x 10 7 

Excessive bandwidth 

A-TDMA/CDMA/ 
SDMA 
R f = 1000 

1 

5 x 10 3 

ion 

Excessive bandwidth and power 

Retained Candidates 

A&E FDMA/SDMA 

1 

0.05, 

10 6 

Most efficient but requires 
ground control of frequency 
management 

E - TDM A/CDMA/ 
SDMA 

Rpp = 100 

33 dB 

25 

2 x 10 7 

Relatively high bandwidth — 
requires CPU time slot 
management 

E-TDMA/FDMA/ 

SDMA 

1 

5.0 

10® 

An efficient approach — re- 
quires CPU time /frequency 
slot management 

A&E- FDMA/ CDMA/ 
SDMA 

43 

0.25 

2 x 10 7 

Flexible approach requiring 
less ground CPU manage- 
ment, but more bandwidth 


The pure FDMA approach while being the least complex for the user and efficiently 
utilizing bandwidth has a requirement for accurate carrier stability^ on the order, of 
10" 8 . Furtherm^fh^'tnmethbd^limihating'fio^ler shift would be necessary which 
would increase the on-board complexity somewhat. Each system would require a 
separate receiver-preprocessor at the ground processor. However, these could be 
implemented on a growth basis as the system requirements grow to the 10® capacity. 

In fact, the system would not be limited to 10® users as long as bandwidth was made 
available for expansion. 
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In conclusion, these two approaches are feasible and worthy of further tradeoff 
against approaches with other configurations because of their bandwidth efficiency/ 
power efficiency, and implementation complexity. 

2.5 CONFIGURATION n ANALYSIS 

Configuration II, shown schematically in Figure 2-24, is characterized by the lack 
of a polling signaling path between the control center and the users. User signals are 
generated autonomously and sent through at least four satellites to the control center. 
Position location is done, at the control center hyperbolically, using time-of-arrival 
difference measurements. Although this configuration can inherently be used for navi- 
gation by providing suitably timed signals from satellites to user and appropriate proc- 
essing equipment on-board, only the surveillance communication problem is considered 
here. 



Figure 2-24. Configuration II 
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Many of the results obtained in the analysis of autonomous polling in Configuration 
I are directly applicable to Configuration n. There are* however, some significant 
differences arising from the different functions performed by the transmitted signals. 
In Configuration I the user transmits a long (about 100 bits) data sequence which con- 
tains the data necessary for the control center to compute his position. No time-of- 
arrival measurements are made directly on the signal. Because of this, a wider range 
of signaling possibilities exists with Configuration I than Configuration IL, 

In Configuration II, the user signals must satisfy two criteria: (1) they must per- 
mit the control center to make time-of-arrival measurements with sufficient accuracy 
to satisfy the system error requirements, and (2) they must, either implicitly by time, 
frequency, or code, or explicitly by data content, identify the transmitting user. The 
first criteria imposes the requirement that the signal be of sufficient effective band- 
width to permit range measurements. This is the essential difference between Con- 
figuration I and Configuration II signaling. In the next subsection the analysis of Sub- 
section 2. 4 is adapted as required to indicate the performance of Configuration II. 

The following basic parameter assumptions will be made where numerical results 
are computed: 

1) B s = minimum bandwidth required to transmit a ranging signal 

= 10 7 Hz 

2) T u = 1 second 

3) M = 20 bits when user ID is transmitted 

= 1 bit where ID is not needed. 

The analyses of the various multiple access-methods have been performed within 
the same framework as that used for Configuration I; the data sheets which summarize 
the results are in Appendix E, Tables E-l through E-6. 

The following paragraphs discuss the significant features of each method. 

2.5.1 TDM A 

As in Configuration I, autonomous TDMA is taken to mean that the users transmit 
autonomously with no provisions for signal separation other than the low probability of 
overlap which is assured by the short pulse duty cycle. The signals must contain ex- 
plicit identification, hence a provision for 20 bits of data is made. 

In autonomous TDMA, the addition o f use rs ca n onl y be done by shortening the 
pulse length for ail users in order to maintain the desired probability of error. An in- 
crease in N, therefore, can only be achieved by a directly proportioned increase in re- 
quired bandwidth and user power. 
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2.5.2 FDMA 


Autonomous FDMA implies that a separate frequency channel is assigned to each 
user. The results in this case are significantly different from Configuration I. Since 
identification of the user is implicit in the frequency channel assignment, the user 
sends only one bit of data per transmission. Although this is reflected in the low user 
power requirement, the bandwidth is not similarly reduced since 10? Hz is required 
per user for transmission of the ranging signal. Increase in the number of users is 
reflected in a proportional bandwidth increase with no change in user power. 

2.5.3 CDMA 

In autonomous CDMA all users transmit in a common frequency channel of width 
B s , with signal separation accomplished by PN code recognition. The transmissions 
are continuous rather than pulsed; and since user identification is implicit in the code, 
only one bit of data is sent per transmission. User power in this approach is independ- 
ent of the number of users. Since each user transmits one bit per second with a 10? Hz 
signal bandwidth, 70 dB signal processing gain is available, which is sufficient to pro- 
vide adequate clutter margin for up to 10 5 users. (Note that the achievement of 70 dB 
processing gain in practice is not feasible. ) 

2.5.4 TDMA/FDMA 

The TDMA/FDMA system arbitrarily assumes a number of separate frequency 
channels. Appendix E shows data for 10, 100, and 1000 channels. The users are 
evenly divided among the frequency channels, and within each channel the users trans- 
mit pulse ranging waveforms autonomously, with the pulse length determined by the 
desired probability of overlap of 10~ 5 . A 20-bit identification is required in each sig- 
nal. As in autonomous TDMA, increasing the number of users is accomplished by in- 
creasing bandwidth, data rate per transmission, and user power. The total bandwidth 
requirement is determined only by the number of users, and not by the channelization. 

2.5.5 FDM A/CDMA 

An FDMA/CDMA system is constructed by channelization of the bandwidth and as- 
signment of equal numbers of users to each channel. Within each channel, signal 
separation is by code division. The data rate is one bit per transmission, and the 
transmission for each user is continuous. As in the case of pure CDMA, the available 
processing gain resulting from the signal bandwidth to data bandwidth ratio far exceeds 
that required. 
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2.5.6 TDM A/CDMA 

In TDMA/CDMA all users transmit in the same frequency channel, and each is 
assigned a different pulse code. In this case, signal overlaps are acceptable, and the 
pulse lengths are determined such that the expected number of overlapping signals (as- 
suming randomly timed transmissions) is one-third of the number that can occur with- 
out exceeding the 13 dB signal-to-clutter ratio. The signal processing gain available 
is determined by the pulse length, since the signal bandwidth is fixed, and in all cases 
it exceeds the required gain for up to 10 5 users. User power is independent of the num- 
ber of users. 

2. 5. 7 COMPARISON AND ELIMINATION 

The critical comparison criteria are required system bandwidth, user power, and 
required signal processing gain. Placing realistic attainable limits on these param- 
eters leads to a first-level comparison which results in the elimination of all but one 
of the six candidates discussed in this section. The limits imposed are 

1) Bandwidth — 10 8 Hz for 10 5 users 

2) Use power — 5 kW peak 

3) Processing gain — 40 dB. 

Autonomous TDMA exceeds the bandwidth and power limits. FDMA exceeds the band- 
width limit (note that this differs from the Configuration I results). CDMA, while 
meeting bandwidth and power requirements requires excessive signal processing gain. 

In the hybrid methods, TDMA/FDMA and FDMA/CDMA both require excessive 
bandwidth. The only method with reasonable values in all three of the primary criteria 
is TDMA/CDMA, for at least 10 ^ time slots. In practical terms, however, a difficulty 
in implementing this approach arises from the fact that it appears that each user re- 
quires a separate code. An alternative to the assignment of separate codes to each 
user is to allow groups of users all having the same code to transmit at slightly different 
pulse repetition rates. (See Subsection 2. 8. 1. 2 for a description of the LIT system. ) 

A preliminary conclusion from the foregoing analyses is that there is no pure or 
hybrid multiple-access approach using only time, frequency, and code orthogonality 
that meets the practical limitations set forth above for the Configuration II communica- 
tion configuration. This conclusion is supported by the fact that the LIT (Reference 9) 
surveillance system, which is a Configuration n system, uses a form of quasi-ortho- 
gonality based upon pulse repetition rate to obtainjthe addi tional ^ si gnal^p^aMllt^ re- 
quired in a TDMA/CDMA Hybrid without using a large number of codes. The above 
analyses indicate that there is no alternative but to introduce such quasi-orthogonality 
in order to obtain a practical system. With modifications of the LIT type to solve the 
code design problem, the TDMA/CDMA hybrid is the strong Configuration II candidate. 
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2. 6 CONFIGURATION El ANALYSIS 


Figure 2-24 can be used for visualization of Configuration in by the addition of an 
interrogation signaling path passing either from the control center through one satellite 
to all users, or through more than one satellite to a specific group of users. Which 
alternative is followed is primarily a question of downlink data rate and satellite power 
limitation, and does not materially affect the results of the multiple access analysis 
presented in this subsection. 

It was seen in Subsection 2.5 that the multiple-access analysis for multisatellite 
autonomous signaling could, with certain modifications, be adapted directly from the 
single-satellite analyses of Configuration I. There is, however, an essential differ- 
ence in elicited-response multiple-access between single and multiple satellite con- 
figurations. This difference is the fundamental limitation on TDMA interrogation that 
exists in Configuration HI, which was discussed in some detail in Subsection 2. 3. This 
limitation makes the analysis framework of Configurations I and II inappropriate for 
Configuration HI. 

The multiple-access methods considered for use in Configuration in are the same 
six as have been analyzed for Configuration II — TDMA, CDMA, FDMA, TDMA/ 
FDMA, CDMA/FDMA, and TDMA/CDMA. Each of them has a somewhat different op- 
erational meaning in Configuration HI. The six candidates are described and analyzed 
in the next subsection. 

2. 6. 1 TDMA 

Elicited-response TDMA signaling in Configuration in refers specifically to the 
time- sequential polling process described in Subsection 2. 3. The limit upon the rate 
at which sequential users can be interrogated was seen to be determined primarily by 
the system geometry, and secondarily by the response message length. It was also 
found that a range-cell polling procedure can be designed to maximize the allowable 
number of users for a given poll repetition rate, T u , and total area of coverage. 

In the analyses of Configurations I and n, the number of users was taken to be the 
primary independent variable, with system bandwidth, processing gain, and user 
power determined as functions of the number of users. It is more convenient in Con- 
figuration III to take T u and the message length, Tj as independent variables and de- 
termine the number of users that can be accommodated assuming optimization of the 
range-cell polling. The details of the analysis are in Subsection 2. 3. 2 and the results 
for a number of cases are given in Table 2-1. For a more complete definition of per- 
formance, these results will be extended to include required user power. 
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In an elicited-response TDMA system, the user response pulses are uniquely 
identified by the receive time slot in which they occur. No explicit identification is 
required in the pulse. It is assumed that each response pulse consists of a short pulse- 
compression PN code which permits full utilization of the channel bandwidth for SNR 
improvement. All users can use the same code, hence the receiver consists of a single 
matched filter. The data rate sent by the user is one bit per pulse. The user pulse 
power required for the postulated user-satellite link (36 dB/Hz/W), assuming a de- 
sired E/N 0 of 13 dB is given by 



-23 + 10 log 


M_ 

T. 

i 


(2-42) 


where M is the number of bits per pulse. Equation 2-42 is plotted in Figure 2-25 for 
M = 1 and also for M = 20, which shows the power requirement for the case in which 
an explicit identification is required in each response. 


It is noted that since all users transmit in the same frequency channel, the system 
bandwidth is independent of the number of users up to the limit imposed by the polling 
restrictions. 


2.6.2 FDMA 


FDMA in the context of this analysis implies that each user responds in a separate 
frequency channel to its interrogation. Since there is no problem of signal overlap, 
the geometric polling restrictions do not apply, and the users can operate independently 
of each other. This further implies that the analysis and results of Configuration n are 
directly applicable. The only operational differences in Configuration HI FDMA are 
that the use of an interrogation signal permits the use of trilateration for position location 
and also establishes a control center-to-user data path which can be used for collision- 
avoidance advisories. 


2.6.3 CDMA 


The same argument holds for CDMA as made above for FDMA. From a signaling 
point of view there is no distinction between Configurations II and III, and the analysis 
of Subsection 2. 3 applies directly. 

2.6.4 CDMA/FDMA 

As in the preceding two cases, the analysis of Configuration II applies directly. 
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2.6.5 FDMA/TDMA 


This method is implemented by channelizing the spectrum into a number, Rjr, of 
frequency channels. Within each channel, interrogation and response proceed accord- 
ing to the geometrical restrictions discussed previously. This hybrid is a straight- 
forward extension of pure TDMA in which additional users are accommodated by in- 
creasing the bandwidth in increments of the basic 10? Hz minimum imposed by the 
ranging accuracy requirements. The total system bandwidth for a given T u and T^ is 
given by 



where [ ] denotes next highest integer. N is the number of users with which the sys- 
tem is to operate, and N]y[ is found from Table 2-1. 

The user power requirements given in Figure 2-25 apply to this case directly. 
2.6.6 TDMA/CDMA 

The basic framework of elicited-response TDMA for Configuration III can be gen- 
eralized to include CDMA in several ways. Three possibilities are considered. In all 
cases the users respond in a common frequency channel of width B^. It was noted in 
Subsection 2. 3.2 that TDMA used alone in multisatellite signaling resulted in an inef- 
ficient use of bandwidth because of the low duty cycle of channel utilization. The use 
of CDMA allows multiple responses to occupy the same time slot, and hence raises the 
bandwidth efficiency to the same order of magnitude as that achieved in Configurations 
I and n. As in the preceding analyses, we are concerned primarily with the system 
utilization of bandwidth, user power, and process gain. The following subsections 
give the pertinent analyses for each approach. These three are considered as separate 
candidates in the comparative analysis of Configuration HI approaches given in the next 
subsection. 

The three approaches are distinguished by the way in which the range cells are es- 
tablished and the different codes are assigned to groups of users. The first approach 
is characterized by a CONUS range-cell arrangement, and a corresponding TDMA in- 
terrogation method. Each interrogation triggers responses from a group of users in 
the same cell, with each user in the group being assigned a different code. In the sec- 
ond approach, the CONUS area is divided into a number of sectors. A range-cell poll- 
ing scheme is organized independently in each sector, and the interrogation in all 
sectors proceed simultaneously. All the users in each sector respond with the same 
code, and are separated by TDMA. Different sectors are assigned different codes. 

In the third approach, no range cells are used. TDMA interrogations take place at a 
rate determined by the CONUS propagation delay, with a large group of users respond- 
ing to each interrogation. 
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2.6. 6.1 CONUS Beam 

Assume a sequential polling signal being transmitted on a CONUS coverage beam 
to all the users in the system. The interrogation points are spaced by a time interval 
determined by the area of coverage, Ax, the interrogation scan time, T u , and the mes- 
sage length T(. To generalize the results we let 

T. = MT d (2-44) 

where M is the number of bits per message and Tj^ is the bit duration. 

The users in each range-cell are divided into K groups, each containing R users. 
Within each group, each user is assigned a different response code. Each polling in- 
terrogation point keys responses from all the users in one group. If N]yj is the number 
of users that can be interrogated using TDMA with an optimum range-cell procedure, 
it allows the system to interrogate x R c users per scan. The clutter which results 
from multiple responses is overcome by the signal processing gain obtained from 
pulse-compression. It is noted that since no assumption of a uniform distribution of 
users among range cells is made, the number of groups, K, will vary from cell to cell 
depending upon the level of cell occupancy. Some cells may have less than R c users. 

In general, a cell will have K - 1 groups of R c users, and one incomplete group. For 
simplicity we assume an integral number of groups. It is noted that this method re- 
quires that as a user moves from cell to cell both his response code and his interroga- 
tion point are changed by control center signaling. 

We are interested in determining the bandwidth efficiency and user power require- 
ment for this approach. Assume that the total channel bandwidth is Bj. Then the sig- 
nal processing gain (PG) achieved by pulse compression is 

PG = b t t d (2 ~ 45) 

The processing gain required is determined by the clutter level. If the R c users who 
respond to a given interrogation are randomly distributed over the cell (because of the 
relatively small cell size this is a reasonable assumption), then the response pulses 
at the receiver occur at random times in an interval t + T^, where r + Tj is the spac- 
ing between interrogations. The expected number of overlapping pulses is (with some 
simplifying assumptions) 

R T. * • 

E [overlapped pulses] = - — — (2-46) 
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For a margin of safety, we allow three times this number of simultaneous pulses to 
occur. Allowing a 13-dB signal to clutter margin, we can then write the required sig- 
nal processing gain as 


(PG) 


3 x 20 x R T. 

c 1 


req 


r + T 


(2-47) 


i 


Equating the required processing gain to the available gain, we have 

7T 


T = 


60 R T. 
c 1 


T D T + T. 


(2-48) 


Using Equation 2-44 and noting that the total number of users N = N|^R C , this becomes 
60 NM 


B T N„ (t+T.) 
M l 


The system bandwidth efficiency, expressed in users per MHz, is then 

N_ 6 _ NM < T+T i> 

B T X ° 60 M 


(2-49) 


(2-50) 


The user pulse power requirement for this method is found as follows* Assuming 
a 20-dB desired E/N 0 , we have 


P T 

E r D 

N ~ N 
o o 


= 100 


(2-51) 


where P r is the received power. Since, with the given link assumptions, we have 
P 


5=- = 4000 T 
N xp 

o 


(2-52) 


then 

T = (40 T-)^* -- --- : — - - t ' - (2-53) 

xp b 
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Figure-of-merit Equations 2-47, 2-49, and 2-53 can be used to calculate a figure-of- 
merit corresponding to that used in the analysis of Configuration I: 


Q = 


T B x PG 
xp T 


M R 


r + T. 


(2-54) 


2. 6. 6. 2 CONUS Beam, Multiple Sectors 

There are operational advantages which can be put forward in favor of dividing 
CONUS into a number of sectors, each under the control of a semiautonomous control 
center. Two advantages which are immediately clear are (1) the number of aircraft 
under the control of each sector control center is smaller than the CONUS control 
center, hence the polling signal design is simpler and the ground processing load is 
reduced over that required for a CONUS system, and (2) there is a close match with 
the existing ATC system which divides the country into semiautonomous sectors. 
Previous discussions of satellite approaches to multisector surveillance have gen- 
erally assumed a multi-beam satellite antenna, with spot beams illuminating each 
sector. In this subsection we describe and analyze an approach which achieves the op- 
erational advantages of multiple sector surveillance, while using a CONUS antenna 
beam for both uplink and downlink signaling. 

Assume that the CONUS area is divided into L equal area sectors. Within each 
sector a control center transmits an interrogation waveform to the users in that sector. 
The interrogation signal can be a different long PN code for each sector occupying the 
same frequency band. Users in each sector acquire the appropriate interrogation sig- 
nal. 


Within each sector, an optimum range-cell polling procedure is followed, with the 
range cells designed on the basis of an area of coverage equal to Ap/L, where Ap is 
the total CONUS area. (See Table 2-1 for typical parameters for L = 10. ) The users 
in each sector are thus separated by TDMA and can respond with the same code. 

There are L codes, one for the users in each sector, and all users respond through the 
same satellite repeater. The required signal processing gain to overcome clutter with 
the same assumptions as the preceding section is 


(PG) 


req 


60 L T. 
i 

r + T 

i 


(2-55) 


and the available processing gain is given by Equation 2-45. Therefore, we have in 
this case, 


60LT. 

B T T D = r+ T. 


(2-56) 
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and since N = L (where Njyj is now maximized over the reduced area), and Tj = MTj) f 


60 NM 

T = N M (T+T i ) 


(2-57) 


This is identical in form to Equation 2-49 except that the same assumptions of T u and 
T[ will result in different performance because of the difference in sedtoT area. Equa- 
tion 2-^50 for bandwidth efficiency is also valid for this case, as is Equation 2-53. 

The expression for the system figure-of-merit is 



(2-58) 


2. 6. 6. 3 CONUS Beam, No Range Cells 

The third TDMA/CDMA hybrid to be considered is constructed by dividing all the 
users over the system coverage area (CONUS) into P groups, each containing R c users. 
Within each group R c response codes are assigned. The code assignments and group 
breakdown are permanent, and are independent of user position. A single interroga- 
tion waveform contains interrogation points spaced by the two-way propagation delay 
across the entire coverage area. For CONUS, this is about 30 ms. Each interroga- 
tion keys responses from all the users in one group, hence P interrogations are re- 
quired for a complete scan. 

The analysis depends upon the assumption that all the users in each group are ran- 
domly distributed over the system coverage area. Then, the Rc response pulses from 
each group will occur at the receiver at random times over a 30-ms interval. By a 
derivation similar to that used in the preceding sections, it can be shown that 


60 R T. 

•d T c 1 

T D T+ T. 

i 


(2-59) 


whence 

60 MR 

T> = £ 

T r + T 

i 

- _ - 60 MN 

" P (T+ T.) 
' r 


(2-60) 
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since N = PR C . In this case there is a simple relationship between P and the scan 
time T u given by 

T 

= 30 x 10 3 second 


and therefore, 

1.8 MN 


= 


T T (T+T.) 
u 1 


The bandwidth efficiency is then 


T (t + T) 

N 6 u' V 6 
— x 10 = — — — — x 10 
B m 1.8 M 

T 


2. 6. 7 COMPARISON AND DISCUSSION 


(2-61) 


(2-62) 


The multiple-access methods discussed in this subsection are compared on the 
basis of bandwidth efficiency, required signal processing gain, and user transmitted • 
power; A discussion is given of implementation feasibility and operational factors. 

In Table 2-6 the comparison parameters are elevated for each method under the 
following assumptions: 

1) T u = 1 second 

2) T| = 10-5 second (where applicable) 

3) M = 1 bit per transmission 

4) Nm = 270 (from Table 2. 1). 

In constructing the table, the methods that incorporate TDMA were analyzed under the 
assumption of an optimum range cell polling procedure (except, of course, TDMA/ 
CDMA Method 3). The assumed restriction that the ranging signal bandwidth cannot be 
less than 10? Hz was also observed where applicable. 

2. 6. 7.1 Bandwidth Efficiency 

As has already been noted, the bandwidth efficiency of methods that use TDMA 
without any element of code division is low, and is fixed by the polling procedure. In 
the cases of CDMA and CDMA/FDMA, the dependence of the bandwidth efficiency upon 
N arises from the assumption of fixed bandwidth, with additional users being accom- 
modated by code division up to the limit of the available signal processing gain. The 
factor of two which separates TDMA/CDMA Method 3 from Methods 1 and 2 is a result 
of the fact that Method 3 does not use range-cell polling. The loss in efficiency in the 
use of range-cell polling results from the extra delays incurred in moving from cell to 
cell and the fact that the number of cells is approximately equal to the number of users 
(see Subsection 2.1.3). 
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Table 2-6 

EVALUATION OF COMPARISON PARAMETERS 


Method 

Bandwidth 

Efficiency 

(users/MHz) 

Processing 

Gain 

Required 

User Peak 
Power (W) 

TDMA + 

27 

1 

500 j 

FDMA* 

0.5 

1 

0.005 

CDMA* 

0.1 N 

20 N 

0.025 

FDMA/CDMA* 

0. 1 n/r f 

20 N/R F 

0.025 

FDMA/TDMA+ 

27 

1 

500 

TDMA/CDMA + 1. 

8350 

1.2 Nx 10“ 3++ 

2500 

TDMA/ CDM A + 2. 

8150 

1. 22 X 10“ 3++ 

2500 

TDMA/CDMA + 3. 

16,650 

0.60 Nx 10~ 3++ 

2500 

♦Continuous transmission and bandwidth fixed by ranging re- 
quirements 

+ Pulse transmission, Tj = 10"® seconds 
++ cf. Equations 2-47 through 2-49 


2, 6. 7. 2 Processing Gain 

To clarify the implications of the processing gain requirements, Table 2-7 shows 
the values for N = 1()5 for those methods that use code division. This table also shows 
the number of codes required to accommodate these many users, which is a significant 
parameter in CDMA methods. 


Table 2-7 

PROCESSING GAIN 


Method 

Processing 
Gain Required 
N= 10 5 (dB) 

Number of 
Codes Required 

CDMA - — — 

-56’---*- - 

io 5= ■ ■ ” : s *"■ ' 

FDMA/CDMA 

56 - 10 log R f 

10 5 /R F 

TDMA/CDMA 1. 

21 

370 

TDM A/CDMA 2. 

21 

L (= No. of sectors) 

TDMA/CDMA 3. 

18 

3000 
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2. 6 0 7. 3 User Power 


Under the above set of parameter assumptions, the required user peak power is 
the same for all methods which incorporate TDMA without CDMA. In CDMA/TDMA 
hybrids there is a 7 dB increase in the required power to reflect the fact that the sig- 
nal to thermal noise ratio is increased to 20 dB in order that the signal -to-clutter 
ratio be the determining factor in error performance. FDMA, CDMA, and their 
hybrid have continuous low power requirements. The average power is comparable in 
all cases. 

Based on the data summarized in Tables 2-6 and 2-7 some preliminary conclu- 
sions can be drawn concerning the relative merits of the methods considered. 

TDMA, FDMA, and their hybrid (TDMA/FDMA) are low in bandwidth utilization 
relative to the other methods. Pure CDMA requires unrealizable signal processing 
gain for a reasonable numbers of users. This problem is alleviated in the FDMA/ 
CDMA hybrid, but at the cost of excessive bandwidth. The best bandwidth utilization 
is achieved by the TDMA /CDMA hybrid, particularly Method 3. The processing gains 
required in all three methods are realizable. The number of codes is, however, very 
high for Method 3, and quite high for Method 1. Method 2 requires only as many 
response codes as there are sectors, typically ten. 

Based on these observations, it is concluded that the TDMA /CDMA hybrid as 
described in this subsection is the most attractive candidate, particularly in the 
Method 2 realization. Methods 1 and 3 cannot be dismissed, however, without further 
investigation of the code design and selection problem. 

2. 7 CONFIGURATION IV ANALYSIS 


The distinguishing feature of Configuration IV is that it implements a surveillance 
system using only a single satellite. Position location using only one satellite can 
only be done by the use of multiple, spaced antennas and interferometry techniques. 
With a two-dimensional antenna array an interferometer can locate a user in two 
dimensions with an accuracy of about 1 km. Altitude data must be measured by an 
on-board altimeter and the reading transmitted to the ground on the same signal that 
is used for interferometer measurements. 

The principles of interferometer position location have been described in some 
detail elsewhere (Reference 10). The significant features for the purposes of the 
present multiple access study are as follows: 

1) Position location accuracy is essentially independent of the signal fre- 
quency, and is dependent upon antenna separation. Using a maximum 
separation of 250 ft, an accuracy of 1 km requires a receiver E/N of 
about 36 dB. ° 
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2) Altimeter data and user identification (when required) must be transmitted 
on the ranging signal. 

3) The use of pseudo-random signaling for interferometer has not been 
investigated in detail. In this report we consider only FDMA, TDMA, and 
their hybrid, both with elicited and autonomous signaling. 

Since only a single satellite is used, the polling considerations discussed for 
Configuration I apply to this case directly. Range-ordered polling for the elicited- 
response case is assumed, thus negligible guard time is allowed between successive 
responses in TDMA. 

The principal operational differences between interferometer position location 
and the other methods considered are the relatively low accuracy, the two-dimensional 
limitation, and the high E/N required to achieve the 1 km accuracy. 

In the calculation of parameters for various multiple-access methods we assume 
a scan time, T u , of 1 second, a 40-bit transmission for altitude and ID, with a proba- 
bility of error of 10“ 5 . This is exceeded if a 36 dB E/N 0 is achieved, where E is the 
total energy in the 40-bit transmission. We also assume a 6 kHz minimum bandwidth 
per user to allow for Doppler effect and oscillator instabilities, regardless of the 
actual data rate. 

The results of the multiple-access analysis are given in Appendix F, Tables! F-l 
through F-5. 


As in the other configurations, it is clear that multiple-access using autonomous 
TDMA requires an excessive amount of bandwidth to achieve the desired low proba- 
bility of overlap (10~ 5 ), and also requires unworkably large amounts of user power 
to obtain satisfactory measurements with the required short transmission times . Pure 
FDMA without Doppler drift correction has a low bandwidth efficiency and requires on 
the order of 600 MHz to accommodate 10 5 users. Pure elicited-response TDMA has 
a high bandwidth efficiency, but the power requirements are very high (100 kW for 
10 5 users. 

The most promising approach lies in the elicited TDMA /FDMA hybrid which rep- 
resents a compromise solution between excessive bandwidth and excessive power. For 
example, with channelization into 100 channels 10 5 users can be accommodated in 
46 MHz of bandwidth at a transmitted power level of 1000 W. 
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2.8 SUMMARY AND CONCLUSIONS 


The use of satellite technology in future ATC systems appears certain. The 
variety of possible satellite ATC systems, although not specifically enumerated herein, 
were referred to by inference and obviously are limited only by the imagination of 
system engineers. Since each alternative provides different levels and categories of 
service and involves different assumptions, the question naturally arises as to whether 
there exists any general method of classifying and evaluating, in a total resource cost 
fashion, the performance of ATC candidates. The results of this study effort provides 
at least a partial answer to this question. Both communications and position location 
techniques have been classified so as to emphasize fundamental technique trade-offs. 

In this section the major communications design decisions have been identified. To 
assure that no significant alternatives were overlooked, general satellite configura- 
tions which would be necessary to support various position locations schemes were 
postulated and investigated on the basis of the impact on attendant communications 
requirements. Collected in this single report is a uniform evaluation of major design 
alternatives pertinent to communications performance, namely, 

1) The polling technique 

2) The multiple access technique 

3) The range measurement technique 

4) The modulation technique. 

As was indicated in the text, perhaps the most fundamental trade-off exists between 
the selection of polling/multiple access pairs. This design decision has the most 
significant effect on the utilization of the basic resources of spectrum bandwidth, uses 
transmitter power, and required processing gain. In those configurations requiring 
range measurements to support position location, it was shown that considerable 
design flexibility in supporting the surveillance function is lost, when compared with 
system design in which the functions of position location and reporting are independent. 

In the course of the study, a figure-of-merit, Q, was defined as a means of indi- 
cating the utilization of transmitter power, bandwidth, and processing gain. Addition- 
ally a uniform analytical procedure ensures that all candidates are evaluated equally. 
On this basis, general approaches to satellite ATC systems have been selected for 
further study. For comparison, the well-known concepts of LIT (location identification 
transmitter) and Woodford have also been analyzed for Q performance. 

The chief benefits of the study lie not only in the identification of promising 
candidates but in the interpretation of the results and the observation of general trends 
accompanying design alternatives illustrated, e.g. , in subsection 2.4.1 Noting that 
Q varies as N^ (N equals the total number of users to be supported by the system) for 
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all candidates excepting those that employ frequency division multiplex, which vary 
only as N, is a strong indication that this conventional approach to channelization for 
tactical systems, may yet be valid. Asa summary to this section we will conclude 
by presenting those basic system approaches which merit further consideration. A 
brief summary description of the LIT and Woodford approaches will be presented first. 

2.8.1 IDENTIFICATION AND DESCRIPTION OF PREVIOUSLY PUBLISHED 
TECHNIQUES 

2. 8. 1.1 Woodford (Reference 11) 


The Woodford system uses three satellites in inclined synchronous orbits. It is 
essentially a scheme for surveillance by two-way ranging. The ground control center 
transmits continuous data and common-timing waveform interrogation signals through 
the three satellites to the users. A separate 1-Mb/s data stream is sent through each 
satellite, thus a total of 3 Mb/s is sent. The binary data in each signal is modulated 
onto a 10-megachip/s pseudo- random sequence. Since 10^ users are serviced, 30 
bits/s is transmitted from the ground to each user. 

Each user receives and demodulates all the data from one of the satellites. All 
data except that addressed to him is discarded. In addition to data reception, the 
received signal is used to time-synchronize a 10-megachip/s pseudo-random sequence 
which is generated by the user and returned through all three satellites to ground. 

This roundtrip synchronization permits two-way range measurements to be made by 
the ground station. 

All the users transmit continuously through all three satellites. The separation 
of these signals at the ground station is done by a combination of space division (mul- 
tiple antenna beams), frequency division, and code division. Each satellite receives 
user signals on three separate antenna beams and on four frequency channels in each 
beam. With 10^ users, this arrangement provides 12 separate channels each of which 
is used simultaneously by about HP users. Within each channel signal separation is 
by code division with an integration time of 0. 1 second, giving a signal processing 
gain of 10 6 . A 20-dB signal to interference ratio is thus theoretically achieved. 

The Woodford system provides surveillance position fixed with an rms error of 
about 100 ft from all error sources. The position update rate is once per second in 
high-density areas and one per 3 seconds, otherwise. Velocity measurements are 
made by Doppler techniques to an accuracy of about! ft/s7~The bandwidth require- 
ments is about 10 MHz at L-band for the user to satellite uplink signaling. 
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This system is essentially a Configuration HI elicited-response polled system us- 
ing hybrid FDMA/CDMA multiple access and common timing waveform interrogation. 
As mentioned in Subsection 2.6.4, the multiple access analysis is the same as that 
for a Configuration n system, and is given by the data of Table E-5, Appendix E, with 
appropriate modification to account both for Woodford's use of three antenna beams 
and for the allowance for a 10 bit/s user-to-satellite data transmission rate. 

As has been discussed in Subsection 2. 7, this communication method is not a rec- 
ommended candidate because of the 60 dB of processing gain that must be achieved by 
the ground receiver. 

2. 8. 1.2 LIT (Location Identification Transmitter) Surveillance 

The LIT surveillance concept is based upon the idea of minimizing the cost of the 
user equipment. The user instrumented only for surveillance carries a transmitter, 
but no receiver. The transmitter sends a pulsed, periodic signal through four satel- 
lites to the control center, where the time-of-arrival differences are used for hyper- 
bolic position location. These signals are generated autonomously by each user. 

All the users transmit in a common frequency channel, with multiple access by 
code division and time division. LIT is a Configuration II system with autonomous re- 
porting and CDMA/TDMA multiple access. To serve 10 5 users with a signal-to-clutter 
ratio of 13 dB, the required processing gain is given by (see Table E-6, Appendix E) 

PG = (SNR)_, X3x|- 
C R T 

Since the LIT transmission uses a 51.1-^ts pulse, 

_ . .3 N 

r t = 20 * 10 - r t = 5 , 

and 

PG = 25 dB 

We have made the assumption that the number of simultaneous users to be overcome 
by code division is three times the number expected on the basis of random distribu- 
tion. 


The LIT system designers have recognized the difficulty of assigning a different 
orthogonal 511 bit code to each of 10 5 users. To solve this problem, the number of 
users is divided into a number of groups (15-20). All the users in a group transmit 
the same code, but each transmits with a slightly different pulse period. The assign- 
ment periods range in 10-/xs increments from 1. 0 to 1. 1 seconds. From the multiple 
outputs of each code matched filter, the individual user signals are identified by syn- 
chronizing with each pulse repetition rate. 


2-77 



By extrapolation of the data in Table E-6, Appendix E, it can be seen that for Rt = 
20 x 103 the required user transmitted power, Txp = 500 W for the postulated link. 

(The LIT design actually provides about 2 kW of pulse power as a result of slightly 
more stringent assumptions about the link parameters, and the provision of operating 
margin). 

2.8.2 RECOMMENDED COMMUNICATION TECHNIQUES FOR FURTHER STUDY 

The systematic approach taken to the analysis of candidate communication tech- 
niques in this section has been comprehensive in the sense that all possible approaches 
to multiple-access and polling have been considered, given the framework of the four 
basic satellite configurations. The intention has not been to select a single technique 
as having overall superiority, but to gain sufficient insight to make a first-level selec- 
tion, from all the techniques considered, of those which are clearly unreasonable in 
the demands that they make upon the basic system resources of power, bandwidth, 
and signal processing gain. Those candidates that remain are considered suitable for 
further study as part of an overall satellite-based ATC system. To summarize the 
results, these selected candidates will be listed and briefly discussed. 

2. 8. 2.1 Configuration I 

The principal feature of Configuration I is its ability to use narrowband signals. 

This is because the satellite-to-ground link carries only data, rather than signals on 
which time-of-arrival measurements are to be made. Of the 12 candidates considered, 
two are considered suitable. These are elicited-response TDMA/FDMA, (E-TDMA/ 
FDMA) and FDMA, both autonomous signaling and elicited-responses (E&A FDMA). 

2. 8. 2. 1.1 E&A FDMA 

Each user in this system is assigned a separate narrowband channel over which 
he continuously transmits binary data representing the processed satellite navigation 
signals. The total user bandwidth requirement for transmission of this data is 100 Hz. 
Thus, the total system bandwidth is 10 MHz for 10 5 users. The feasibility of this ap- 
proach depends upon the correction for Doppler frequency uncertaintly, which is in the 
6 kHz range. If Doppler correction cannot reasonably be applied, the required system 
bandwidth goes up by a factor of 60, which exceeds the bandwidth that can be considered 
available for the system. A second possible drawback is the requirement to implement 
10 5 separate ground receivers. 

2. 8.2. 1.2 E - TDMA/FDMA 

” The diff lenities encountered in implementing the pure FDMA approach in Configu- 
ration I can be alleviated by the introduction of an element of time-division into the mul- 
tiple access. Range-ordered polling using discrete address interrogation is used to 
elicit responses sequentially from groups of users in separate frequency channels. The 
time-sharing of channels permits a more coarse channelization tha n in pure FDMA, and 
removes the stability problem of very narrowband signaling. The number of ground re- 
ceivers is reduced. The user set must, however, incorporate a receiver with address- 
recognition ability so that the response can be made to the interrogation waveform. 
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A general conclusion regarding Configuration I communication techniques is that 
when narrowband signaling is possible, it results in higher bandwidth efficiency than 
code division for multiple-access, as is clear from the tables of Appendix A. 

2. 8. 2. 2 Configuration II 

The principal feature distinguishing Configurations II and III from Configuration I 
is that time-of-arrival measurements must be made on the transmitted signal, which 
must therefore be relatively wideband. This requirement eliminates the advantage 
gained by narrowband signaling in Configuration I, and CDMA becomes preferred over 
FDMA in hybrids with TDMA. 

2. 8. 2. 1.3 A-TDMA/ CDMA 

Of all techniques considered under Configuration n, this technique proves to be 
the only one which is reasonable in its demand upon system resources. It is typified 
by the LIT surveillance system, and features autonomous pulsed signaling by all users 
in the same frequency band, with signal separation by code division. 

Operationally, a significant advantage of this approach is that the user requires no 
receiver, and only an inexpensive transmitter. The autonomy of operation imposes the 
requirement for hyperbolic position location, with its GDOP deterioration relative to 
trilateration. 

2. 8. 2. 3 Configuration m 
2. 8. 2. 3.1 E-TDMA/CDMA 

The elicited responses of Configuration HI are limited in pure TDMA by the polling 
restrictions discussed in Subsection 2.3. A single frequency channel cannot be fully 
utilized with TDMA, and additional users are added by code division to increase the 
bandwidth efficiency of the system. A pure TDMA system of this general type is de- 
scribed in Reference 12, Appendix V, although the polling limitations were not at that 
time systematically explored and the potential bandwidth efficiency of the hybrid ap- 
proach was not realized. 

The range-cell polling approach, which was first described in Reference 12, gives 
this system considerable flexibility as shown by the three approaches described in Sub- 
section 2.3. Operational desirability is further enhanced by its use of trilateration for 
position location. 
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2. 8. 2. 4 Configuration IV 

Surveillance systems using a single-satellite (interferometer) position location dif- 
fer fundamentally from multisatellite systems in (1) being unable to locate the user in 
three dimensions, hence being dependent upon pressure altimeter data, and (2) being 
restricted in accuracy to the order of 1 km, rather than the 100 ft range achieved by 
multisatellite position location. 

The single-satellite signaling path bears strong similarities with Configuration I 
in allowing range-ordered polling and narrowband signaling. FDMA systems are feasi- 
ble, with the most desirable multiple access approach being TDMA/FDMA hybrid. 

CDMA signaling was not investigated for this configuration, but the same arguments 
that lead to its rejection in Configuration I can in general be applied. Further study 
of interferometer signaling methods is needed. 

2.8.3 COMPARISON OF RECOMMENDED SYSTEMS 

In this subsection we present and discuss comparative data on the five recommended 
polling and multiple-access techniques. The comparison is relative to the following 
set of primary criteria: 

1) Bandwidth efficiency (users/MHz) 

2) User transmitted power 

3) Satellite power for polling 

4) Signal processing gain required. 

The comparison data assumes 10^ simultaneous users for each system, and makes 
certain other assumptions regarding the specific configurations of each system. These 
assumptions are listed below with pertinent discussions. The basic comparison data is 
given in Table 2-8. This table also includes the same data for the surveillance portion 
of the LIT and Woodford systems. They are assumed to be operating with the same 
communication link and signal-to-noise margins in order to provide a common basis 
for comparison. 

2. 8. 3.1 Configuration I — A-FDMA 

As has been mentioned previously, the principal difficulty with realizing pure 
FDMA in Configuration I is that each user must occupy a band of only 100 Hz. This 
requires continuous Doppler correction, and may be difficult to implement economi- 
cally^ JnTable 2 -t 8 two Jiguresare, given fombandwidth efficiency. The larger num- 
ber assumes that the 100-Hz-per-user bandwidth is realized, and the smaller assumes 
that each user is allotted 6 kHz of bandwidth to allow for Doppler uncertainty and oscil- 
lator instability. The reduction in efficiency is sufficient to eliminate the second case 
from serious consideration. An attractive feature of this basic approach is the low 
peak power requirement and relatively simple user sets. 
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Table 2-8 


COMPARISON DATA FOR SURVEILLANCE COMMUNICATION TECHNIQUES 


Method 

Bandwidth 

Efficiency 

^xp 

<W) 

Rsat 

<W) 

PG 

(dB) 

Q 



Peak 

Avg. 




CI A-FDMA 

4 

10 assuming 

0.5 

0.5 

0 

0 

5 x 10 6 


Doppler 

correction 





3 x 10 8 


168 otherwise 





. 

Cl E-TDMA FDMA 

10 4 

50 

0.5 

.15 

0 

5 x 10 8 

(R F = iooo) 







CII A = TDMA CDMA 

10 4 

500 

0.025 

0 

25 

1. 5 x 10 12 

(R T = 20,000) 







LIT 

10 4 

500 

0.025 

0 

25 

1. 5 x 10 12 

CHI E-TDMA CDMA 1, 

8350 

2500 

0.025 

200 

28 

1. 8 x 10 12 

(T. = 10-5 s) 







Cm E-TDMA CDMA 2. 

8150 

2500 

0.025 

200 

28 

1. 86 x 10 12 

CIII E-TDMA CDMA 3. 

16,650 

2500 

0.025 

200 

22 

0. 9 x 10 12 

Woodford 

2. 5 x 10 3 

0.25 

0.25 

5000 

60 

2 x 10 17 

CIV E-FDMA TDMA 

io 4 

100 

10 

15 

0 

io 9 

RF = 1000) 








2. 8.3.2 Configuration I-E-TDMA/FDMA 


Channelization of the spectrum and the use of TDMA with elicited responses within 
each channel eases the difficulties associated with very narrowband signaling. The fig- 
ures shown assume 1000 frequency channels, each supporting 100 users operating in a 
TDMA mode. The satellite power given assumes that discrete-address range-ordered 
polling is used. This is a worst-case, since there are several approaches that can be 
used to simplify the polling procedure. First, the use of a coarse range-cell configu- 
ration (see Subsection 2.3) with perhaps six cells permits nonrange-ordered polling 
within each cell with about 5 ms guard time. Each user transmission takes 10 ms, 
thus the efficiency is reduced by 50%. This reduction is the cost of simplifying the 
ground processing by eliminating the need for range-ordering the users at every scan. 
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A second simplification results from noting that users in different frequency channels 
can have the, same address. Each interrogation address can therefore key responses 
from up to 1000 users, with a corresponding reduction in satellite power. 

Finally, we note that the basic channel bandwidth in this case is 10 kHz. The 
reduction in capacity incurred by allowing a margin for Doppler is much less severe 
than in pure FDMA. 

2. 8. 3. 3 CH A -T DMA /CDMA 

The basic configuration of this system is similar to that of the LIT system. 
Twenty-thousand time slots are assumed, (i.e. , a 50 /xs pulse length). The process- 
ing gain required assumes that three times the average number of overlapping pulses 
are allowed for, assuming random user distribution. 

It is noted that although the transmission time per user is reduced by a factor of 
200 from the Cl TDMA/FDMA system, the required pulse power only increases by a 
factor of ten because this system sends only one bit of data per transmission. 

2. 8. 3. 4 Cm E-T DMA /CDMA, 1, 2, and 3 

-5 

The basic assumptions for these three methods is a 10 second response pulse 
and one bit of information per transmission. These approaches have been discussed 
and compared in subsection 2. 6. 7. They differ primarily in the polling method used. 
In general they have, relative to the others, considered lower bandwidth efficiency 
and higher user pulse power requirements, because of the short pulse length. The 
pulse can be lengthened only at a reduction in bandwidth efficiency. They use higher 
satellite polling power to allow for the interrogation point updates which are required 
with range cell polling. 

2.8. 3. 5 CIV E -TDMA/FDMA 

The system evaluated assumes 1000 channels of 10 kHz bandwidth. This approach 
combines high bandwidth efficiency with moderate user power and low-satellite polling 
power. The same polling simplifications discussed above under CI-TDMA (FDMA) 
can be applied. The advantages of this system are purchased at a cost in position 
location accuracy relative to the multisatellite configuration. 

2.8.4 OPERATIONAL FACTORS 

,_ _ In_addition to J;he^strictlyd;echnical-criteria-of Table 2-8^ comparison of the rec- 
ommended approaches should also consider certain operational factors. 

An important point is the essential difference between autonomous and elicited 
systems. An elicited system has the intrinsic ability to return data to any user by 
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way of the polling waveform. In an autonomous system such a data path must be 
provided separately if collision avoidance is to be implemented for all users, and the 
additional cost to the user for this communication capability offsets the principal 
advantage of an autonomous system, which is low cost to the user. Elicited systems, 
while inherently more flexible, demand considerable total satellite-to-user data 
capacity depending upon the information content of the polling message. Hybrid 
autonomous/elicited systems may provide a viable compromise wherein elicited res- 
ponses are used only in conflict situations. 

Configuration I implies that the user is performing at least some processing and 
storage of the navigation data received from the navigation satellites, or other source 
(e.g, , OMEGA). 

This further implies some additional cost to the user although the polling and 
response function which take place through the response satellite can be implemented 
using relatively unsophisticated narrowband user hardware. The cost of processing 
and storage of time-of -arrival data could offset the advantage of user set simplicity 
obtained by narrowband signaling. 

Except for Configuration IV, the position location accuracy performance of all 
configurations is roughly equivalent. Section 3 of this report analyzes this aspect of 
system capability in detail. 

The question of the desirability of a single central control center compared to a 
decentralized system with a number of control centers controlling various sectors of 
the country is not one that is specifically addressed by this report. It seems apparent, 
however, that autonomous reporting lends itself to a central control, whereas elicited- 
response systems are more adaptible to the requirements of multiple sectors. This 
is particularly true of Configuration III, in which (see subsection 2. 4) the flexibility of 
range-cell polling, together with CDMA, permits the simulation of a multiple-beam 
sector system using a simple CONUS Beam antenna. 

2.8.5 DISCUSSION 

Table 2-8 gives the basic trade-off data for the surveillance communication 
techniques selected from the wide range of possible techniques considered. The selec- 
tion was based upon the achievement of reasonable values in each of the comparison 
criteria when implemented with 10 5 simultaneous users. Although the number of can- 
didates was thereby significantly reduced, those remaining represent a broad spectrum 
of possibilities for specific systems. 

The figure-of-merit, Q, is a convenient measure of overall use of the basic system 
resources of power, bandwidth, and signal processing gain, all of which have bounds 
which cannot reasonably be exceeded. This measure clearly shows the desirability of 
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the narrowband techniques that are possible in Configurations I and IV. The need for 
additional bandwidth to transmit ranging signals in Configurations II and HI is the 
principal factor which causes Q to increase for these techniques. The signal-process- 
ing gain required by the Woodford system results in a Q four orders of magnitude 
greater than the worst of the recommended alternatives. 

To pursue the development of satellite-based surveillance, future effort should 
concentrate upon further evaluation and comparison of the candidate systems repre- 
sented by Table 2-8. Several specific steps should be taken, such as 

1) Detailed design of a specific system representative of each candidate 

2) Evaluation of the resulting systems relative both to the criteria of 
Table 2-8 and such additional operational criteria as necessary. 

Specific additional criteria of prime importance are user hardware 
requirements, satellite design requirements, position location 
performance, relative weighting factors, and additional consideration 
of the satellite to user link design. 

3) Merging the surveillance communication system with an area 
navigation system to provide full capability. 
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Section 3 


SURVEILLANCE/NAVIGATION ACCURACY PARAMETRIC ANALYSIS 



Section 3 


SURVEILLANCE /NAVIGATION ACCURACY PARAMETRIC ANALYSIS 
3.1 INTRODUCTION 

The purpose of this section is to provide a parametric error analysis of the four 
surveillance/navigation concepts given in Table 3-1. 

Table 3-1 

SURVEILLANCE /NAVIGATION CONCEPTS 


Name 

Description 

Representative System j 

3P 

Three range difference 

621B 

2P+H 

Two range difference plus 

— 


Altitude 


3R 

Three range 

IBM (Fourth generation ATC proposed) 

2R+H 

Two range plus altitude 

PLACE 


P is defined as the range difference between two known points and an unknown (the 
navigator); i. e. , P = Rj - Rq. One possible implementation is to measure R^ and Rq 
separately via one way transmission. However, by doing so, a false geometric range 
is embedded in the measurement which is caused by a difference in the time reference 
between the transmitting and receiving systems. Since the false geometric range is 
common to both Ri and Rq, the differencing operation removes this false range. How- 
ever, it can be shown that instead of differencing R^ and Rq, a different algorithm can 
be written which not only solves for the navigator’s position but also for the false range, 
which is equivalent to solving for the reference time difference. The navigator’s position 
accuracy is not affected by the choice of algorithm. Therefore, though referred to in this 
report as the 3P concept, the reader should understand that there are two choices of im- 
plementing the 3P system depending primarily on whether there is any advantage to solv- 
ing for the time bias between transmitter and receiver. 

The sources of error affecting the accuracy of these concepts include 

1) Random errors — range, range difference, and altitude 

2) Bias errors — range, range rate, and altitude rates 
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3) Ionospheric refraction — time and space correlations modeled 

4) Satellite ephemeris — position and velocity with spatial and time 
correlations. 

The following three types of positions fix accuracy were investigated (see Figure 

1) Absolute — position fix relative to the satellite support system 

2) Relative — position fix relative to a local earth surface system; 
e. g. , an airport 

3) Resolution — distance between two aircraft in close proximity; 
i.e. , less than 10 miles. 

$ * 6 I 
\ - 1 



Figure 3-1: . Pictorial Representation of Absolute, Relative, 
and Resolution Frames of Reference 
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The accuracy of all concepts is highly geometry dependent; that is, for a fixed 
range measurement error the accuracy of the position fix varies with the location of 
the user with respect to the satellites . Some concepts are more susceptible to this 
variation in geometry than others. This geometrical effect is often referred to as 
geometrical dilution of precision (GDOP). However GDOP is but a small part of the 
total problem. GDOP gives the effect on each component of the users position, (i.e. , 
east, north, and vertical), to each and every source of error listed above. Even if 
all errors could be treated as either statistically independent or perfectly correlated 
a myriad of numbers would be needed to provide a complete error analysis. When 
compounded with the added dimension of correlated error sources, a direct assault 
on the problem would produce an unwieldy and impractial result. 

Therefore, a practical and meaningful condensation was required which did not 
compromise the primary objective of presenting the results parametrically. In par- 
ticular, the relative performance of concepts subjected to similar geometries and er- 
ror models was considered essential. 

The presentation of the results is structured to emphasize the relative perfor- 
formances among concepts. For example, the reader can compare the performance 
of the 3P and 3R systems due to an ionospheric ranging error of 1 ft and that due to 
an ephemeris error in the same figure. For example, if the ionosphere were a truly 
homogeneous stationary phenomena one might expect the ionosphere to affect the 3P 
system less than the 3R system because in the 3P system ranging errors (in this case 
the ionospheric range errors) common to all transmission paths cancel. However, the 
ionosphere is not uniform nor stationary and the GDOP is typically several times more 
for the 3P system than the 3R system. The net result is that the 3P system perfor- 
mance with respect to ionospheric effects is inferior to the 3R system and this will be 
evident from data presented in the same figure. 

All results have been normalized; no judgements or conclusions have been reached 
with regards to the magnitude of each error source because this would necessarily im- 
ply a certain bandwidth, transmission frequency, etc. Rather, the reader may select 
the proper normalized results, apply the proper magnitude of error his design war- 
rants and then combine everything for a total system error. 

3.2 DEFINITION AND GENERAL EXPLANATION OF TERMS FOUND IN RESULTS 


Due to the rather extensive amount of curves and bargraphs it is impractical to 
discuss each figure individually. Therefore, at the outset, the background informa- 
tion necessary for a proper interpretation will be given. Detailed mathematical ex- 
planations have been intentionally avoided in favor of verbal descriptions. Those 
readers interested in a complete mathematical development are referred to Appendix 
G. 
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3.2. 1 OUTPUT COORDINATES, TYPES OF ACCURACY (ABSOLUTE, RELATIVE 
RESOLUTION) 

All results are expressed in terms of the user ! s horizontal position/velocity errors 
and his vertical position/velocity errors. This reference frame seemed most appropri- 
ate in the context of an Air Traffic Control (ATC) system where, on the average, there 
is no preference for accuracy in a north/south (NS) or east/west (EW) direction. The 
horizontal error is a composite of the EW and NS components. This combination of 
EW and NS components also condenses the output; i. e. , number of graphs by 30%. 

The horizontal error could have been expressed as the rms of the EW and NS com- 
ponents, or as a circular error probable (CEP). CEP was not chosen because when the 
ratio of the major to the minor axes of the error ellipse exceeds 3 to i, the CEP indi- 
cates a value much smaller than the major axis. For example, let the ratio of the ma- 
jor to minor axes be 10 to 1, with the major axis equal to 1, a very likely condition 
with some concepts in certain geometries. Then, CEP = 0.68 and rms = 1. 005. The 
rms is always an excellent measure of the major axis (the largest error) whereas the 
CEP is not. Thus in the context of an ATC system evaluation, the rms must be con- 
sidered more suitable. 

The results may be interpreted in three reference coordinate frames; the absolute, 
the relative, and the resolution. Figure 3-1 illustrates these three frames pictorially. 

A complex of ground stations (the support system) locates, time synchronizes, and con- 
trols all of the satellites. This function may be a service to the ATC system and not 
inherently part of it, or it may be an integral part. In either case, all aircraft and 
local ground stations initially determine their positions with respect to the support 
system. This position fix is absolute and all errors made by the support system are 
passed along to all users. The distance between the support system and users may be 
on the order of several thousand miles, although it could be less. When two users 
are within several hundred miles there is a tendency to have not only the same errors 
affecting the position fix but also the same geometry between user and satellites. The 
position fixes of the two users would tend to be nearly the same under these conditions, 
A cancellation of the common sources of error is therefore anticipated. The distinc- 
tion between relative accuracy and resolution is not as well defined. The distinction 
intended is that, for resolution, all the position errors cancel except those due to 
purely random errors. When the separation distance is but a few miles; e.g. , 10; the 
resolution value probably applies. At greater distances there are incomplete cancel- 
lations because the errors are not exactly the same and the GDOP is somewhat differ- 
ent. Under these conditions the idea of relative accuracy is perhaps most appropriate. 

To avoid any confusion in the int erpretation, the results are not-referred t o.bv 
frame of reference but rather by error source and whether a benchmark calibration 
is, or is not, assumed. The final interpretation, based on the separation distance of 
interest, is left to the reader. 
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3. 2. 2 INTERPRETATION OF BENCHMARK CALIBRATION RESULTS 

No attempt was made to prejudge the results by assuming that certain error sources 
had cancelling effects. For example, the ionospheric errors were not assumed to have 
complete cancellation due to the close proximity of two users, be they two aircraft or an 
aircraft and a ground control center • i. e. , a benchmark calibration station. 

The essence of the benchmark calibration is that two users in close proximity are 
affected by the same error sources, in the same way, and in a similar geometry situa- 
tion. However, one of the users; i. e. , the control center; knows of its location by other 
means. After computing its position via radio signals any difference between the computed 
and known positions is attributed to a composite affect of the many error sources and that 

unique geometry between user and satellites. This difference (a vector) is considered valid 
over some specified area and any other user in that area is assumed to have been affected 
in exactly the same way. i.e. , the same vector error applies. Therefore the position 
fixes via radio signals are corrected by the correction vector of the benchmark station. In 
reality this correction vector is composed of many error sources but to present the re^ 
suits parametrically the composite is broken down into individual sources of error. This - 
approach permits an evaluation of the benchmark calibration technique by error source, 
and indicates that some error sources may not be as amenable to benchmark calibration 
as others. 

3.2.3 ALTERNATE CALIBRATION TECHNIQUES 

Benchmark calibration is not the only practical means of calibration. In fact, a poten- 
tial weakness to the benchmark calibration is that the benchmark correction vector must 
closely match the correction vector of any other user. Figure 3-2 j.s an example which 
illustrates the potential weakness of benchmark calibration. Here D co and D ca are of the 
same magnitude but of different orientation. The result is a sizeable residual calibration 
error. The reader may visualize a situation where -D co is colinear with D ca , in which 
case the calibration adds twice the position error instead of removing it. The possibility 
also exists that the benchmark calibration could work better with one error source than 
another. For example, the ephemeris errors may well be remqved by the benchmark 
calibration, whereas the ionospheric errors may not. If a certain error source ; e. g. , 
the ionospheric error; does not respond well to the benchmark correction, an alternative 
is for the user to make an ionospheric range correction via a model of the ionosphere for 
each and every transmission path. This correction would eliminate approximately 90% 
of the ionospheric range error. Then, when the benchmark correction is applied, the 
ionospheric error effects havebeen reduced to 10% of their original levels ; e. g. , in 
Figure 3-2, the magnitude of D co and D ca have been reduced 10 times, and hence D c is 
also 10 times less. 

This calibration technique uses a model of the ionosphere which may require consider- 
able software. For the surveillance function this would appear to be very practical since 
the model could be shared by thousands of users. However for navigation, the technique 
may be impractical. 
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(Actual Error Made by User A) 


Figure 3-2. Benchmark Calibration Techniques 

Another possibility is that employed by the TRANSIT navigation system which uses 
two carrier frequencies and the fact that ionospheric refraction is inversely propor- 
tioned to frequency squared. By measuring the range difference at these two frequen- 
cies and using the refraction-frequency relationship the refractive effect is computed 
and therefore available as a correction. The disadvantage of such an approach in an 
ATC system is the wider transmission bandwidth required. 

3.2.4 DESCRIPTION AND DISCUSSION OF EACH CONCEPT 

The discussion thus far applies to all concepts, but before discussing error 
sources, it might be well to first discuss the position location concepts. 

3. 2. 4.1 3P Concept 


Position location is derived from the mutual intersection of three hyperboloids. 
The foci of these hyperboloids are at the satellites. Figure 3-3 shows a two dimen- 
sional representation. 



Figure 3-3. Two-Dimensional Pictorial Representation of 3P System 
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Point A in Figure 3-3 is at the intersection of the two range difference contours 
that are equal and zero; i.e. , point A is equally distant from the three satellites. 

Point B is at the intersection of two unequal range difference contours. Each inter- 
section defines a unique yz position which is equivalent to a horizontal and vertical 
position fix. Position errors can occur because 

1) The measurement of a range difference yields the wrong contour. 

2) The contours are misplaced due to an error in specifying the 
location of the foci. 

Though Figure 3-3 is not drawn to scale the dependence of position error on geom- 
etry can be demonstrated. Assume that the contours shown are increments of 1 ft. 
Further, assume the true position is at A such that P^ 2 = ^23 = However, the user 
measures P 2 g = -2 and P 12 = 0. (P being the symbol for range difference.) The user 
computes his position at B. Notice that most of the error is in the z direction. As 
another example, assume the true position is at C. The true Ps and P.g = 2 and 
P 23 = 2 * If the same P error is made the same as the previous example, the measured 
Ps would be Pj 2 = 2+0 = 2 and P 22 = 2-2 = 0. The indicated position is at D. Notice 
that though the vector magnitude of the error is about the same there is considerably 
more error in the Y direction. 

Errors in the location of the foci can be considered as the combination of two effects: 

1) If a consistent or systematic error is made in specifying the location 
of all three foci, the net effect is a translation of the hyperbolas with 
respect to the earth. This does not change the shape of nor distort 
the hyperbolas. Nevertheless, this systematic error causes a user's 
position error. Ephemeris errors of this kind are said to be perfectly 
correlated and could be caused by errors in the support system which 
are common to all satellite ephemerides; e.g. , location of slave sta- 
tions with respect to the master. 

2) The other extreme, error wise, is if estimates of the satellite locations 
are randomly errored. This causes a distortion of the hyperbolas. Fig- 
ure 3-4 illustrates this possibility, The distortion in this case is caused 
by a spreading or tilting of the baselines . 

The user assumes the foci, which are derived from the ephemeris (satellites), to 
be at 1 and 2 which uniquely specifies the location of all range difference contours on 
the users map. However, the satellite is actually at 2'. Therefore, when the user 
measures P = 0 he places himself on the P = 0 contour shown by the solid vertical 
line. However, to have measured P = 0 the user must have been on the P = 0 contour 
shown by the dashed vertical line. The important point here is that the location of the 
foci is critical to the establishment of the correct range difference contours, and hence 
the users position. 

If the foci were in error such that satellite 2 was really at 2” the hyperbolas 
would be distorted and rotated as shown by the dotted contours. 
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P = 0 P = 0 


Figure 3-4. Geometrical Description of Ephemeris Error on 3P Position Location 

In general case, all three effects may be present; translation, spreading, and rota- 
tion. These three effects are statistically embodied in the covariance matrix of satel- 
lite ephemeris errors, a 3n x 3n matrix ( n number of satellites). 

The reader can visualize the complex interaction of ephemeris errors on position 
location. More specifically, if four satellites are involved, a complete’ description of 
the ephemeris position errors is a 12 x 12 matrix which has 78 entries (accounts for 
the symmetry of the covariance matrix and ig nores the velocity terms and their inter- 
action with position). The 78 terms are the position variances of each component of 
each satellite (12 terms) and the covariances of all components with each other (66 
terms). A term-by-term breakdown would yield unwieldy results of doubtful value. 
Nevertheless, the number of terms must be reduced for practical reasons, whereas 
for the completeness of the analyses all terms should be considered. The problem 
was resolved as follows: 

1) An error ellipsoid is assumed which has the same shape, size, 
and orientation for all satellites. The ellipsoid can be circular, 
pancake, or cigar shaped. The choice of shape would depend on 
how the support system determines the ephemeris. Figure 3-5 
shows some typical ellipsoids that could be specified. The 
assumption of same shape and size stems from an assumption 
that all satellite ephemerides are determined by a common 
support system. 

2) Each ellipsoid can, in a broad sense, be related to the kind of 
satellite tracking system used by the support system. For 
example, the pancake shape (T = 1, k^ = 0. 1, ^ = 1) could be 
representative of a trilateration system. 7 The radiar axis (R) 
is about equal to the range accuracy and, with reasonable plac- 
ing of the ground stations, the accuracy of the T and N axes 
would be about equal, and considerably greater than the R axis 
due to the GDOP. Furthermore the size, shape, and orientation 
of all satellites should be approximately the same because the 
satellite ephemerides are being determined by a common system. 
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Figure 3-5. Specification of Satellite Error Ellipsoids 


The T (tangential) axis is assumed to be in the plane of motion and orthogonal to the 
radius vector, R, (radial axis) which is parallel to the geocentric radius vector, and the 
N (normal axis) is assumed to be normal to the plane of motion and forms an orthogonal 
triad. Thus the error ellipsoid is specified by the size of the R, T and N axes (refer 
to Figure 3-5). However, the size of the T axis is, by definition, unity and the relative 
size of the R axis and T axis, k-. and kg units, respectively. In other words, the ellip- 
soid is normalized to the size of the T axis. It is also assumed that the principal axes 
of the ellipsoid are aligned with the RTN axes. This is equivalent to assuming zero cross 
correlations among the R, T, and N components. The shape and orientation of the ellip- 
soid is specified by k-^ and k2- 

The cross correlation (Z 3 ) among ellipsoids is accommodated by assuming the R axis 
error of satellite 1 is correlated with the R axis error of satellite 2 but not with the T and 
N axes of satellite 2; similarly* the T axis and of satellite 1 is correlated with the T axis 
of satellite 2 but not with the R and N axes, etc. The correlations among satellites is 
assumed to be equal; but the degree of correlation can be varied from -1 to +1. 

The specification of the four satellite ephemeris errors has been reduced from 78 
terms to 3 terms while retaining the essential features. For example, if P= 0, the 
ephemeris errors are uncorrelated and the random effects are emphasized. If p = 1 
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the ephemeris errors are perfectly correlated and the systematic effects are emphasized. 
By normalizing to the T axis the parametric form of the analyses is retained and by vary- 
ing kj and kg the shape and orientation is varied. 

Each range difference error is composed of two sources of error: (1) a random 
range error, and (2) an ionospheric range error. * Each source in turn consists of two 
errors introduced on the two transmission links. In mathematical terms. 


Range difference = P. = R. - R„ 

3 1 0 


j 1 


where 


j y 7 3 — 

AR q = AR 0 (NOISE) + AR q (REFRACTION) + AR BlAg 


j=l,2,3 

(3-1) 

j=l, 2,3 

(3-2) 

r bias 

(3-3a) 

r bias 

(3-3b) 


The noise term (the random error) poses no particular analytical difficulty because 
AR. is statistically independent of AR., which is independent of AR^. The assumption 
is ijiade that the variance of this random error is the same for all transmission paths. 
This assumption is reasonably valid because the dominant effect on the range variance is 
the SNR which is inversely proportional to the range between transmitter and receiver. 
The range variations among all signal paths is small. The following abbreviated analysis 
supports this assumption: 


-i" TT K ’ 

R(NOISE) = 1_ = 1 j (3-4) 

SNR SNR s 

where Kj is a constant related to the range resolution. However , from the one-way radar 
range equation 

SNR = ^ (3-5) 


♦Tropospheric effects were not considered per se. This is not necessarily to imply 
that tropospheric effects are negligible. However, thejwork-scope had.to.be restricted 
to some level due to limited funding. The ionospheric model, modified for its effective 
altitude and time and space correlation, could represent a facsimile troposphere. How- 
ever, no analyses were made based on this assumption. 
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where K 2 is another constant depending primarily on carrier frequency, power, and 
antenna gains. Thus, 


2 

^R (NOISE) 




(3-6) 


The variances for range 1 and 2 are 
<r = A R^ 

o-R 2 2 ■ A Rg 


(3-7) 


The ratio of °R / °R is 
1 2 

= Rl/R 2 


(3-8) 


The maximum range difference that can reasonably be expected for the COLM geom- 
etry is about 10 7 ft. The average range is about 10 8 ft. Thus, 


^R/ 


R„ 


~ 10 8 + 10 7 
' 10 8 


= 1.1 


(3-9) 


Generally, this ratio will be less, especially for those satellite constellations (see 
Section 4) with small inclination angles which give relatively close satellite spacing. 

The refraction error poses an analysis problem because (1) the ionospheric range 
error depends on the angle at which the ray enters the ionosphere, and (2) the range 
error for path 1 is correlated to some degree with the range error for path 2. From 
Equations 3-2, 3-3a, and 3-3b, the reader can see that if the two ionospheric range 
errors were perfectly correlated and had the same entry angle ARj = ARq, the result 
would be a cancellation of the ionospheric range error* 

Rather than assume this condition the ionospheric error model provided for 

1) Time and space correlation 

2) Entry angle effect 

3) Latitude dependency (the ionosphere is known to produce less refraction 
at higher latitudes). 

This model was derived directly from Reference 1. 
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The variance of the ionospheric error is defined as a iqnO equator, and 

°^IONO e ~ a< ^ a * latitude j > . ’’a” is constant which makes e~ a <f> = 0.5 when <f> = 45° 
The range variance of any ray passing through the ionosphere at any angle other than 
perpendicular iso- foNO e “ a ^ Q( E ) equals unity when E = 90° , and about 2. 5 

when E = 20°. Q is sometimes referred to as the obliquity factor. 

The covariance of ionospheric range errors 1 and 2 is a ^jqno ea ^l.Q( E i) e a< ^2 
Q(E 2 ) P (S) P{ t) where P (S) is the space correlation and P (t) is the time correlation. 
The correlation models used (taken from reference 1) are shown by Figure 3-6 and 
3-7. 



Time Separation (h) A June A 


O September • 

O December ■ 

Figure 3-6. Ionospheric Time Figure 3-7. Ionospheric Space 

Correlation Correlation 

The separation distance S of Figure 3-7 is illustrated by Figure 3-8. Since S de- 
pends on (Reference 1 gives Hjq^q ~ 200 nmi miles) and P (S) in turn depends 

on S, H iono should be treated as a parameter, because for any given geometry of 
satellites and aircraft lowering Hjq N q should make S smaller, thereby increasing the 
space correlation P (S) . It can be shown that increasing the space correlation causes 
more of the ionospheric range error to cancel. 

Since the ionosphere constitutes a relatively thick shell about the earth ranging from 
an altitude of 50 mi (D layer) to an altitude of 500 mi (F 2 layers), was varied 

between 50 and 400 nmi. 


Reference 1. C. C. Chen, Range Difference Error Due to the Presence of 
Ionospheric Electrons, TRW Systems Group. 
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Figure 3-8, Pictorial Model of Ionosphere 


3. 2.4.2 2P-KH Concept 

The significant difference between the 3P and 2P+H concepts is that the 2P+H 
concept uses the earth radius whose locus is a sphere in lieu of the third hyperboloid. 
Figure 3-3 may again be used to visualize the concept. Point A is the intersection of 
Pl2 (or P 23 ) with a sphere-of-radius (r e + h). An r e is assumed, and h is provided by 
some form of altimeter. Point B is the intersection of P^ with a radius of h plus A h. 
If the user makes an altitude error, Ah, it may produce a horizontal error, but the 
vertical error is always equal to the altitude error. Therefore, the 2P+H results do 
not show any vertical error results. The sensitivity of horizontal position errors to 
both P and H errors is also geometry sensitive. 

In general, the comments on the 3P system with regards to random errors, ephe- 
meris errors, and ionosphere errors are applicable to the 2P+H system. One less 
satellite is involved which means one less signal path, but this is properly accounted 
for in the results. 

3 . 2. 4. 3 3R Concept 

A position fix with the 3R concept involves the intersection of three spheres whose 
centers are located at each satellite respectively. A two dimensional representation is 
shown by Figure 3-9. 
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Figure 3-9. Two-Dimensional Pictorial Representation of 3R System 


Point A is the intersection of R^ = 11 and R 2 = 11. If, in the measurement of R 2 , 
an error of minus 1 ft were made, then R 2 (measured) = 10 and the estimated position 
would be at B. Note that both horizontal and vertical errors result. 


The sensitivity of position errors due to range errors is again geometry depend- 
ent, as will be shown, but not to the extent of the 3P concept. This is because the 
range loci of the 3R concept are spheres which have a uniform curvature whereas the 
3P concept loci are hyperboloids which have a variable curvature. 

Furthermore, the sensitivity is proportional to the angle of intersection of the two 
loci. A right angle intersection is best; i.e. , least error magnification or GDOP. 
Though the reader should not draw too quantitative a conclusion from Figures 3-3 and 
3-9, it is fair to say that, in general, the 3R concept will generally have intersections 
closer to right angles than the 3P concept. Based on these simple observations one 
would expect the 3R concept to be more accurate on the average, and also to have less 
variation in accuracy than the 3P concept. 


In other words, if all other things were equal, and a histogram were prepared of 
all the sensitivities for all possible geometries, one could conceivably generate a plot 
similar to that shown in Figure 3-10. 

Though the preparation of the complete histogram would have been prohibitively 
expensive for this study; there are two parineters which "can be used to quantitatively 
characterize the histogram. These two parameters are mean and standard deviation. 
The mean deviation gives the average sensitivity (or accuracy) over the geographic re- 
gion of interest while the standard deviation shows how much the sensitivity might vary 
within the region. A small mean deviation implies an accurate system and a small 
standard deviation implies a dependable accuracy. 
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Magnitude of. Sensitivity ► 


Figure 3-10. Typical Histogram of Error Sensitivities for 3P and 3R Concepts 

Therefore the mean and standard deviation are provided for each concept and for 
each source of error. The reader should find them useful in making comparisons of 
concepts and how those concepts behave to the different sources of error. In addition, 
results are provided which show in detail how the concepts behave as a function of geom- 
etry and source of error. 

The 3R concept is sensitive to all error sources identified for the 3P system plus 
one — the range bias. It may be seen from Figure 3-9 that a range bias causes, pri- 
marily, an altitude error with no horizontal error when directly under the satellites. 

As the user moves away from a point directly under the satellites the horizontal error 
due to a range bias increase. 

3.2.4. 4 2R+H Concept 

The 2R+H concept is very similar to the 3R concept with the one exception, that, 
in effect, the third satellite is located at the earth T s center, and in lieu of the third 
range an altitude measurement is substituted. In general, the 2R+H concept is sensi- 
tive to all error sources identified for the 3R concept. However, the vertical error 
always equals the altimeter error. 

3.2. 5 SOME ADDITIONAL EFFECTS RELATING TO EPHEMERIS AND IONOSPHERIC 
ERRORS WITH BENCHMARK CALIBRATION 

This subsection has been reserved until the background information related to the 
ephemeris and ionospheric errors and the benchmark calibration was presented. 

As discussed in Subsection 3.2.1 the relative accuracy is based on a correction 
vector from a local ATC center. If the total system were time stationary the correc- 
tion would be valid over an indefinite time period. However, there are three things 
which negate the idea of a time stationary system: 

1) The satellites are not truly stationary with respect to the earth (see Section 
4), which causes the GDOP to vary with time. 
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2) The satellite error ellipsoids are not time stationary. It can be shown 
theoretically, that in an energy conservation system the uncertainty of 
the total state vector is time invariant; but the position and velocity un- 
certainty grows. Thus while the position ellipsoid and velocity ellipsoids 
constantly change their shape and orientation the correlations between 
position and velocity change accordingly so that the total state uncertain- 
ty is time invariant. 

3) The ionosphere as shown by Figure 3-7 is nonstationary. 

The first item above is treated by providing results for different satellite positions. 

The different satellite positions are dictated from an initial choice of constellations 
(see Section 4) and the time of day. The second item is treated in the results as fol- 
lows: the user does not move after the calibration; however, the calibration vector 
should be recomputed continually. However, if only updated periodically, the effect 
of time correlation on the ephemeris error should be considered. Due to the limited 
funding it was practical to consider only two cases which are the extremes of correla- 
tion- In the first case, the satellite ephemeris error is assumed to be uncorrelated 
after some time (unspecified in this report), and in the second, the ephemeris is as- 
sumed to be perfectly correlated. The actual degree of correlation depends on how 
the support system estimates the ephemeris. This aspect was beyond the scope of 
this study. In a practical sense the updating interval could be sufficiently short so 
that the perfectly correlated results are probably closer to realism. 

The third item has been treated much like the second except that the time correla- 
tion of the ionosphere has been varied from zero to one. The reader can interpret these 
results by considering a user initially calibrated while waiting at the end of the runway 
prior to takeoff. After takeoff the aircraft circles the airfield (so as not to affect the 
GDOP) using the satellite radio signals for navigation and the correction vector regis- 
tered at takeoff. However, because the ionosphere is not stationary, the correction 
vector changes but the user is unaware of the change, or at least incapable of compen- 
sating for the change. The error analysis then shows how the relative accuracy varies 
with time variations in the ionosphere. 


3.2. 6 ANALYSIS OF DETERMINATION OF AIRCRAFT VELOCITY FROM DOPPLER 
MEASUREMENTS 


Appendix G gives a mathematical development showing how the error in measuring 
Doppler affects the determination of aircraft velocity. For the 3P and 2P+H concepts 
the Doppler measurement is actually a Doppler difference measurement: 



o 

where R. 

) 


^true + ^BIAS 



(3-10) 
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RBIAS is, in a way, one transmission, due to unknown unequal LO frequencies in the 
transmitter and receiver. As seen from Equation 10, the bias in the 3P and 2P+H con- 
cept cancel. This bias does not cancel in the 3R and 2R+H concepts. (Like the position 
algorithm if desired, the RbiAS ma y solved for which would then allow the user to 
synchronize his LO with the satellite LO. ) 

The derivations in Appendix G show that the propagation of Doppler or Doppler dif- 
ference errors into velocity errors is exactly the same as for the propagation of range 
or range difference errors into position errors, with one addition. The error in range 
(or range difference) which produces a position error also produces a velocity error. 

This effect in a well engineered design is usually not dominant, but if precise Dop- 
pler is to be measured and very inaccurate ranges are measured because the primary 
intent was to have a highly accurate estimate of velocity, the velocity error could be dom- 
inated by the position error. Thus, the effect of range (or range difference) errors on 
velocity have been investigated and are on the order of 10“3 to 10-5 ft/s/ft of range error. 

3.3* RESULTS 

3.3.1 COMPARISON OF 3P, 3R, 2P+H and 2R+H CONCEPTS VERSUS ERROR SOURCES 
AND BENCHMARK CALIBRATION 

Figures 3-11, 3-12, 3-13, and 3-14 provide a summary comparison of the four sur- 
veillance/navigation concepts employing multiple satellites. The reader is referred to 
subsection 3.2.4 for a description of each source of errors. The solid bars give the 
mean (average) accuracies over the surveillance region and the white bars give the 
standard deviations. 

Three satellite constellations are given for each concept corresponding to the mid- 
night Y constellation (see Figure 6-a) with inclinations of 10°, 20°, and 30°. Other 
figures will give comparable results for the Y, 4 A.M. constellation. 

The ordinates of all graphs are in dimensionless units (ft/ft). If a satellite to user 
radio link were designed to provide a range standard deviation of x ft, the users posi- 
tion error would be the value of the ordinate times x. Note also that the normalization 
has been reduced to one radio link. In some concepts 2, 3 or 4 satellites may be in- 
volved; however, the curves properly account for the number of satellites by a statis- 
tical combination of the errors on all of the necessary radio links. 

For example, consider the 3P case, i = 10°. Four satellite radio links are re- 
quired. Assume each link was contaminated with 10 ft of range jitter (a random type 
error). Then, on the average, the horizontal position error will be 460 ft and the var- 
iation from this average (mean) depending on where the user is in region of interest, 
will be 108 ft. (For Figures 3-11 through 3-14, the region of interest is defined by 
Region 2, Figure 4-1.) In other words, the position error could vary from 460 - 2(108) 
ft, = 244 ft to 460 + 2(108) ft = 676 ft. Since the complete histogram has not been com- 
puted, there is no quantitative correlation of probability with standard deviation; i.e. , 

3 <T does not imply 99.7%. The standard deviation given should be considered more in 
the context of how much the position error is likely to vary from the mean at any given 
location. 
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Figure 3-11. Comparison of 3P, 3R, 2P+H arid 2R+H Concepts vs Error Source 
without Calibration/Horizontal Position, CONUS, Midnight Y, ki = . 1, ko = 1 
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without Calibration; Vertical Position, CONUS, Midnight Y, k* = .1, k 2 = 1 
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Assuming the same radio link design, if applied to the 3R concept, the average hori- 
zontal position error would be 52 ft with a 2. 2-ft variation. The ionospheric range error 
results have been normalized to 1 ft. The actual ionospheric range error is inversely 
proportional to the carrier frequency squared. At 1. 6 GHz, the ionospheric ranging error 
is between 10 and 100 ft. 

3.3.2 COMPARISON OF 3P, 3R, 2P+H, AND 2R+H CONCEPTS VERSUS VALUE OF 
BENCHMARK CALIBRATION AND SIZE OF AREA TO BE SERVICED BY 
CALIBRATION 

Based on these first four figures the benchmark calibration appears to be of little 
value. However, only one benchmark was assumed for the entire region (region 2, 

Figure 4-1) which may be too large an area. Therefore the region was reduced to a 300 
x 300-nmi square (a 5° x 5° grid). The center of the 5° x 5° grid is at the center of the 
larger grid (CONUS) and approximates the area about Chicago. The results are shown in 
Figure 3-15. The solid bars are for the 5° x 5° grid and the white bars are for the area 
given by region 2. In most cases the value of benchmark calibration is doubtful even for 
the 5° x 5° region. The increase for random errors closely approximates V3F because 
of the RSSing of the user and benchmark errors. However, the results for the ephemeris 
and refraction are somewhat surprising and the reasons for the rather poor improve- 
ments will be investigated in subsequent figures and discussion. 

Ini the case of the ephemeris error the results have been presented for a worst-case. 
The results assume that the correlation of the current correction vector is uncorrelated 
with the one made at the time of the benchmark calibration. In practice the correlation 
is determined from the type of satellite location system employed, the length of the filter 
span in the orbit determination algorithm, and the time from calibration. Hence these 
results could be loosely interpreted as a long time (i. e. , several hours) between bench- 
mark calibrations as opposed to a continuous type of benchmark calibration. 

3.3.3 EFFECTS OF EPHEMERIS ERROR ON HORIZONTAL POSITION 

Figure 3-16 shows the effects of ephemeris error on horizontal position assuming 
(1) no calibration, (2) calibration with zero correlation, and (3) calibration with per- 
fect correlation. Figure 3-16 thus presents the effects of the ephemeris error three 
different ways: (1) the top bar graph represents the effects of the ephemeris error 
assuming no calibration; i. e. , the user position error is in an absolute frame; (2) the 
middle bar graph represents the effects of the ephemeris error for the previously des- 
cribed situation; i.e. , calibrated but the correlation has gone to zero; and (3) the bot- 
tom bar graph represents the effects of the ephemeris error with calibration and 
where the correlation is unity; i. e. , equivalent to a frequent updating of the bench- 
mark calibration. * " 
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Figure 3-16. Effects of Ephemeris Error on Horizontal Position, Region 2 


The bottom bar graph indicates that the benchmark calibration is almost always 
beneficial when the correlation is high or equivalently when the benchmark calibration 
is performed frequently. Figure 3-16 also assumes a single calibration station for 
the entire Region -2 (see Figure 4-1)^^ - ^ " r - — - 

3.3.4 HORIZONTAL ERROR DUE TO EPHEMERIS 

Figure 3-17 shows the horizontal error due to ephemeris within a 300 x 300 nmi 
grid versus time from calibration and presents the effects of the ephemeris error for 
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a smaller region (5° x 5°) with correction correlation equal to one. A somewhat dif- 
ferent aspect of the calibration update interval is apparent. In Figure 3-16, the cor- 
relation of the benchmark calibration is based on the correlation of the ephemeris er- 
rors themselves, but the movement of the satellites was not factored into the results. 

In Figure 3-17, the movement of the satellites is emphasized. The movement of the 
satellites is shown by Figure 4-2. As the satellites move, even if the ephemeris 
errors are perfectly correlated, the correction vector of the benchmark calibration 
will change. This is due strictly to a change in the relative geometry among satellites 
and users. If this correction vector is not updated frequently, even with perfectly 
correlated ephemeris errors, there will be an error introduced in the presently com- 
puted position. Figure 3-17 shows the value of the benchmark calibration when used 
over a 5° x 5® region. The degradation of the correction is only moderate over a 3- 
hour interval. The time span was chosen where the relative satellite motions would 
be greatest; i.e. , the midnight constellation. 

3.3.5 EFFECT OF CHOICE OF SATELLITE ERROR ELLIPSOID ON THE HORIZONTAL 
ERROR 

Another consideration affecting the ephemeris error results is the choice of error 
ellipsoid (see Figure 3-5). The results shown thus far assume the error ellipsoid sec- 
ond from the left in Figure 3-5. It is the pancake shape with the radial error 10 times 
smaller than the tangential and normal errors. Figure 3-18 compares the results of 
the pancake shape ellipsoid with the cigar shape (first figure on the left in Figure 3-5) 
ellipsoid. In Figure 3-18, the solid bar represents the pancake shape and the white 
bar represents the cigar shape. The differences are seen to be small even though the 
normal error has been reduced by 10 to 1; i.e. , the volume of uncertainty of the pan- 
cake shape is greater than for the cigar shape. 

3.3.6 EFFECTS OF IONOSPHERE RANGE ERROR AND MODEL ON POSITION 
LOCATION 

Figure 3-19 shows the effects of the ionospheric refraction on the four position lo- 
cation concepts. The results are presented in the context of mean and standard devia- 
tions discussed previously. The region of calibration is a 5° x 5® grid (300 x 300 nmi). 
Figure 3-19 indicates that the 3P system accuracy is actually degraded when the cali- 
bration is applied, and that the improvement for the 3R, 2P+H and 2R+H concepts is 
marginal. The variation in accuracy is due to the time correlation of the ionosphere 
and the movement of the satellites with respect to the earth. 

A simplified mathematical explanation of why the calibration does not improve the 
accuracy is given below. Consider first the range difference error due to refraction: 

AP 2 = AB 1 -AR q (3-11) 

( 3 - 12 ) 



P(B) 


(3-13) 
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Figure 3-18. Effect of Satellite Error Ellipsoid 
on the Horizontal Error 


Equation 3-11 is a simple statement of the range difference error; Equation 3-12 is 
the ensemble average of AP 2 and Equation 3-13 is but a redefinition of terms in Equa- 
tion 3-12. cr and <tr q are the variances due to the ionospheric refraction and P( S) 
is the correlation distance (see Figures 3-7 and 3-8). Neglect the effect of latitude on 
the ionosphere and assume the obliquity factor is the same for paths 1 and 0. Then, 


o R 


1 


a R 0 
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Time-After Calibration (h) Time-After Calibration (h) 

Figure 3-19. Comparison of 3P, 3R, 2P+H and 2R+H Concepts due to Refraction 
with and without Calibration/Vertical Position, 

300 x 300-nmi Region, Midnight Y, k^ - - . 1 
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Under this assumption Equation 13 becomes 
2 2 

O’p = 2 o- r [1 - ?(S)] (3-14) 

P(S) is determined from the relative positions of the user, the satellites and the iono- 
sphere. Figure 3-20 gives a simplified geometrical picture. 



Figure 3-20. Simplified Diagram Illustrating Correlation Distance 


Applying a straight line fit to Figure 3-7 gives 
= 1 " 2400 

and substituting Equation 15 into 14 gives 
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(3-15) 


(3-16) 


(3-17) 


whence 

- o- R (0. 17) 

2 

In other words about 17% of the refraction error remains in crp. However, the GDOP 
(defined in Subsection 3. 1) is very nearly the same as the sensitivity of the random 
error. From Figure 3-11, we have GDOP = 46. Therefore, the horizontal position 
error would be 46 (0. 17) = 7. 8 which agrees quite well with the value given in Figure 
3-11; i.e., 7.6. 
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When i = 30° , crp = [0. 09] and a p = 0. 3 ctr. The GDOP for i = 30 is seen to 

be 4. 7 and hence the position error due to refraction is (4.7) (0. 3) = 1.4 which again 
agrees closely with the result given in Figure 3-11; i.e. , 1. 77. 

10 

When i = 10 c , S ~ 200 (5774), which is about 35 miles. Apparently the error can- 
cellation property of range differencing does not counterbalance the loss in accuracy 
due to GDOP. In other words, to have good ionospheric error cancellation characteristics 
the satellite spacing must be very close which imposes a severe GDOP penalty. A similar 
argument can be applied to explain why the calibration is not too effective even for short 
distances. The calibration depends on a very high correlation of the range difference 
error to be effective because, due to the large GDOP, any residual refraction errors 
are multiplied by the large GDOP to yield a large position error. The severity of the 
effect is shown in Figure 3-21. This figure shows the growth of the position error as 
the user moves away from the calibration point. 



Figure 3-21. Effect of Refraction Error on Horizontal Position 
with and without Calibration as a Function of Distance from Calibration Point 

Figures 3-22, 3-23, and 3-24 explore another aspect of the refraction problem. 
From the previous discussions and results, the impact on error sensitivities due to S 
has been established. One of the parameters directly influencing S is the height of the 
ionosphere, hjoNO* Since there is no clearly defined altitude at which the refraction 
takes place, hjQ N Q was varied. 

The results of this parametric analysis are summarized in Figures 3-22, 3-23, 
and 3-24 with good agreement with the previous results; the 3P concept is very sensi- 
tive to anything that changes P( S). Furthermore, the lower hjQNQ Is, the less change 
there is in ^(S) and the 3P system should indicate better performance. The figures 
are in agreement. 
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Figure 3-22. Effect of Height of Ionosphere on Horizontal Error 

without Calibration 
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Figure 3-23. Effect of Height of Ionosphere on 
Horizontal Error with Calibration 
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3.3.7 EFFECT ON HORIZONTAL ERROR DUE TO ROTATION OF SATELLITES IN 
THE Y CONSTELLATION 

Figures 3-11, 3-12, 3-13 and 3-14 assumed a midnight Y constellation. Figure 
3-25 and 3-26 assume exactly the same conditions except that four hours have passed 
and the Y constellation has moved from its midnight to its 4 AM configuration, (see 
Figure 4-2a). The results are very nearly the same for i = 10° but show some differ- 
ences for i = 30° . The explanation lies in noting from Figure 4-2a that for i = 10° the 
Y shape is maintained throughout a 24 hour period whereas for i = 30° , the Y shape is 
transformed into an inverted T after 4 hours. The change from the Y to the T shape 
is very likely the major reason for the differences. If the horizontal errors were 
broken into two parts; i. e. , EW and NS; there would be a noticeable shift in the accu- 
racy of the EW and NS components with time. But since the EW and NS components 
are RMSed, the rotation of the Y shape (i = 10°) does not have a major influence. The 
shift would be most noticeable in areas of low elevation angle where there is a tendency 
for the EW and NS errors to be unequal. 

3. 3. 8 HORIZONTAL ERROR DUE TO ALTIMETER ERROR WITH AND WITHOUT 
CALIBRATION 

The 2P+H and 2R+H concepts use an altitude measurement in lieu of a P or R mea- 
surement. An error in the altitude measurement produces an error in the horizontal 
measurement (see Subsection 3. 2. 4. 2). Figure 3-27 addresses this aspect consid- 
ering both the 2P+H and 2R+H concepts with and without calibration. The results 
show a nearly one to one effect; that is, 100 ft of altitude error causes 100 ft of hori- 
zontal error. The value of the calibration is marginal, but there was only one cali- 
bration station assumed for the entire Region 2. Considerably better results are pos- 
sible when the serviced region is substantially reduced. The black bars represent 
the mean deviation and the white bars the standard deviation. 

3.3.9 MISCELLANY 

3. 3. 9.1 3P Performance With Box Constellations 


The 3P concept can, for certain satellite constellations, yield accuracies which 
are much worse than those indicated thus far. In particular, the box constellation 
as its name implies arranges the satellite into a box pattern. In Figure 4-26, the 3 
AM or 3 PM constellation with the center satellite removed is one example of a box 
constellation. In Figure 4-8, 12 midnight case, if satellite 4 (or 5) were located at 
the bottom of the figure 8, a box constellation would result. 

Figure 3-28 illustrates quantitatively how disastrous the consequences are when 
a box constellation occurs. The top figure should be compared with the two bottom 
figures. All things are equal except for the arrangement of the satellites. 
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Figure 3-26. Comparison of 3P, 3R, 2P+H and 2R+H Concepts vs 
Error Source with Calibration, 4 AM Y Constellation 
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It is important to realize that, though the satellites are arranged as the corners of 
a rectangle, the baselines upon which the 3P concept derives its range differences and 
hence its hyperboloids are as, shown in Figure 3-29. 



Figure 3-29. Baselines Formed frpm Box Constellation for 3P Concept 

When e is about 120° , we have the midnight Y configuration and when e is about 
90° we have the 3 AM configuration (see Figure 4-20). From this, one could infer 
that the overall performance of the 3P concept is strongly related to e . Figure 3-30 
has been constructed to show this relationship. 



Figure 3-30. Effect of 3P Performance vs 
the Angle between the Three Baselines 


At zero and 180° the baselines are parallel, which explains the shape of the curve 
there. The flatness of the curve between 90° and 170° is somewhat surprising but in- 
dicates a very uniform performance of the 3P concept for that range of e . Certainly 
e <60° should be avoided. There is no known simplified explanation for the perfor- 
mance singularity at e = 45° . An isolated case can be constructed which shows graph- 
ically that when €~ 45° and the user is at the center of the box there is no vertical in- 
formation, i.e. , the intersection of three hyperboloids is a vertical line. This iso- 
lated case also helps to explain why the 2P+H concept does not have a singularity, the H 
provides the vertical information. 
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SATELLITE COVERAGE 


4.1 INTRODUCTION 

All the surveillance/navigation systems that have been investigated, require si- 
multaneous visibility to two, three, or four satellites (more than four is always re- 
dundant). A problem exists in how to deploy the least number of satellites to provide 
visibility over the largest desired area of the globe. There are two classic approaches 
to this problem: (1) combinations of 24-hour period satellites in geostationary and in- 
clined/elliptical orbits which form stationary constellations over a fixed geographic 
region, and (2) less than 24-hour period satellites in inclined, near circular orbits with 
uniformly spaced orbit planes. The latter produces a nongeostationary pattern of sat- 
ellites over any geographic region and, in fact, tends to provide global coverage 
though not perhaps on a continuous basis. 

Many investigations have been performed for Defense Navigation Satellite Systems 
(DNSS) with regards to the satellite coverage problem and both the stationary and nonsta- 
tionary cases have been investigated in depth'*'. However, these studies have had global 
coverage as the principal requirement. In this study we were interested primarily in 
coverage for a specified and limited region of the globe: i.e. , the 49 continental United 
States, Puerto Rico, and Hawaii. 

Three test cases, or geographic regions were established (see Figure 4. 1). Region 
1 covers CONUS and southern Canada. Region 2 extends the boundaries of Region 1 to 
include most of Alaska, Puerto Rico, and Hawaii. Region 3 covers CONUS, Hawaii, 
Puerto Rico and a very narrow strip of Canada along our border, but not Alaska. 

Initially, it was throught that the stationary constellations would prove superior 
for the limited coverage problem primarily because (from the DNSS studies) it was known 
that global coverage requires about 20 satellites and the total number of satellites was 
about the same for either the stationary or nonstationary constellations. The stationary 
constellation, on the other hand, requires only four to six satellites. Therefore, based 
solely on the number of satellites in orbit, the stationary constellation appeared more 
attractive. 

However, since the objective of this study was to perform a parametric treatment 
of all aspects of the COLM/ATC problem, the nonstationary constellations were also 
investigated. 


^"621B Geometric Performance for a 20 Satellite Global Deployment, !f Aerospace 
Report No. TR-0172 (2521-01)-!; November 5, 1971. 
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4.2 RESULTS 


4.2.1 STATIONARY CONSTELLATIONS 

Coverage for two types of stationary constellations are presented in this subsection. 
The first type is the so called ft pin wheel" constellation which has a geostationary sat- 
ellite in the center and three or four 24-hour period, inclined, and eccentric satellites 
which appear to an earth observer to produce a nearly circular pattern about the geo- 
stationary satellite (see Figures 4-2, 4-3 and 4-4). The satellite movement in 1-hour 
increments shows the effect of the eccentricity of the orbit which varies as e = tan i/2 
in order to produce the circular ground trace. Notice that the Y or X formations are 
retained for all time for i = 10°, whereas there is considerable change for i = 30°. The 
change is caused by the eccentricity which causes a speeding up in the southern hemi- 
sphere and a slowing down in the northern hemisphere. 

The second type is a variation of the pinwheel, which uses two geostationary sat- 
ellites and three or four elliptical/inclined satellites, the latter of which appear to an 
earth observer as an egg-shaped pattern about the two geostationary satellites. 

The pin wheel constellation has two parameters to be varied: (1) the longitude of the 
geostationary satellite, and (2) the inclination of the satellites in the ring. The longitude 
of the geostationary satellite was chosen to provide maximum coverage of CONUS. The 
inclination was then varied and the contours of visibility were established. The contours 
are defined as the locus of points where the earth observer can see the required number 
of satellites, and where at least one of the satellites is at an elevation of 10°, which is 
generally accepted as a practical radio navigation horizon. These contours establish the 
boundaries of multiple satellite visibility. The results for the four satellite pinwheel con- 
stellations are shown in Figure 4-5. The contours are shown for the worst case; i. e. , - 
when one of the satellites is at perigee or at the point of maximum southern latitude. 

This condition occurs three times daily. Otherwise, the visibility contours would ex- 
tend farther north. However, since an ATC system must function continuously, the 
extended coverage is only of secondary value. 

From Figure 4-5, it can be concluded that for CONUS coverage, with any position 
location system using four satellites, the inclination must be less than 20° , whereas to 
provide coverage for CONUS, Hawaii, and Puerto Rico the inclination should be about 
10°. For position location systems using only three satellites there is considerably 
more coverage in the northern latitudes, but still not adequate for coverage of both CONUS 
and Alaska (see Figure 4-6). Figure 4-6 is shown for the 12 noon configuration (Figure 
4-3) using satellites 2,3, and 4, which provide the best accuracy but least coverage. 

For all other time configurations, the three-satellite visibility contours are essen- 
tially that for the geostationary satellites because, as seen from Figure 4-3, the satellite 
south of the equator would not necessarily be used. 
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Figure 4-3. Four-Satellite Pinwheel Constellation 



Figure 4-4. Five-Satellite Pinwheel Constellation 
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Figure 4-5. Four-Satellite Pinwheel Constellation, Four-Satellite Visibility Contours 



Figure 4-6. Four-Satellite Pinwheel Constellation, Three-Satellite Visibility Contours 
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If a fifth satellite is added to the pinwheel constellation, there is an improvement in 
the coverage, particularly in the northern latitudes. This configuration is shown in 
Figure 4-7. 


Figure 4-7 assumes the 3 AM configuration (see Figure 4-4), and the position fix 
would be derived from satellites 1, 3, 4, and 5 when west of satellite 1, and from satellites 
1, 2, 3, and 4 when east of satellite 1. The midnight and 4 AM, 6 AM, and 8 AM config- 
urations will provide slightly more coverage at higher latitudes. However, Figure 4-7 
is based on the premise that 100% visibility is required, at all times. 

Three satellite visibility contours are essentially the same as for the geostationary 
satellite since the satellite in the southern hemisphere is not visible at the northern lat- 
itudes. Note, however, that at 4 AM and 4 PM, satellites 3 and 4 are quite close, which 
affects the position accuracy. 



Figure 4-7. Five-Satellite Pinwheel Constellation, Four-Satellite Visibility Contours 

Figure 4-8 shows another variation of the stationary satellite constellation wherein 
the inclination is very high (60°) and two geostationary satellites are placed inside the 
ring. The eccentricity for this constellation exceeds 0.5. Consequently, perigee alti- 
tude is only about 7000 miles which results in some loss of visibility and higher Doppler 
rates when viewing those satellites during their passage in the southern hemisphere. 
However, as Figure 4-8 shows, it is never necessary to use a satellite south of the equator 
for position fixing. 
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Figure 4-8. Six-Satellite Constellation with Large Inclination and Eccentricity 

Figure 4-9 shows the coverage for six-satellite constellations. The 3 AM con- 
figuration was selected because it has the largest spread, thereby providing the least 
amount of coverage but the best position accuracy. Figure 4-9 also assumes that air- 
craft west of the center of the geographic coverage use satellites 1, 2, 5, and 6, whereas 
those east of center use satellites 1, 2, 3, and 4. Combinations like 1,2, 4, and 5 
provide very poor accuracy (see Figure 3-30). 

Another stationary constellation is shown by Figure 4-10. This constellation uses 
two geostationary satellites and three others in inclined circular 24-hour orbits. The 
equatorial spacing of the geostationary satellites should be about equal to the inclination 
angle so as to provide an equal position accuracy in the EW and NS directions. 

The augmented Figure 8 five- satellite constellation has a noticeable disadvantage in 
that the 4 AM and 4 PM configurations have no satellite well north of the equator. There- 
fore, during these two times- each dayr the=NS positiomaccuracy will be from two to three 
times worse than the EW position error. At other times, the EW and NS position errors 
will be about equal. 
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Figure 4-9. Six-Satellite Constellation with 
Large Inclination and Eccentricity, Visibility Contours 
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The coverage and effect of varying the inclination angle of the Figure 8 five-satellite 
constellation is shown in Figure 4-11. 



Figure 4-11, Figure 8 Five-Satellite Constellation Visibility Contours 

Comparison with Figure 4-7 (five-satellite pinwheel constellations) shows very com- 
parable coverage characteristics. 

In summary, all of the stationary constellations considered exhibit similar coverage 
characteristics for the same number of satellites. The fifth and sixth satellites provide 
additional coverage especially in the northern latitudes. However, no single stationary 
constellation provided coverage for CONUS, Puerto Rico, Hawaii, and Alaska simultane- 
ously. 

The reader can visualize sliding the coverage plots in a longitudal direction from 
whence it can be seen that Alaska and Puerto Rico can never be covered with a single 
stationary constellation. 

4.2.2 NONSTATIONARY CONSTELLATIONS 

Since these constellations are not stationary, they tend to provide coverage of the 
entire globe. Therefore, to provide 100% coverage over any given area is almost equiva- 
lent to providing coverage for the entire globe. Since this requires more than five or six 
satellites, why consider the nonstationary cases? 
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Besides the reason of performing a parametric study, there are such reasons as 
(1) a single stationary constellation could not provide coverage for CONUS, Puerto Rico, 
Hawaii, and Alaska; (2) the nonstationary constellation will have the satellites at less 
than synchronous (24 hours) altitudes thereby having potential benefits in terms of launch 
costs and less free space radio propagation loss; (3) the nonstationary case will have a 
very soft degradation as the failure of any given satellite produces only intermittent cover- 
age problems; (4) each satellite can be independently commanded ON or OFF which implies 
the potential of increasing coverage when and where necessary or desired without a marked 
increase in system costs. The satellite would be commanded ON when entering the desired 
coverage area and OFF when leaving. This feature is particularly easily implemented with 
polar orbits where a single station at Fairbanks or Thule can "see" each satellite once every 
orbit. Simple timing mechanisms reset at each north polar passage would provide the OFF 
command for the south going passes. Increased coverage would be provided by altering the 
ON/OFF commands; e.g. , a longer On time on the south-going pass would provide cover- 
age for many Latin American countries. 

As was shown in Subsection 4.2.1, the stationary constellation coverage deficiency was 
decidedly in the northern latitudes. Therefore, the choice of nonstationary constellations 
focused on correcting this deficiency. This was done by concentrating on polar orbits, 
though orbits with 40° inclinations were also briefly investigated. The polar orbit grossly 
"over covers" the polar regions at some expense to the equatorial regions. The 40° incli- 
nation restores some of the imbalance. However, the inclined orbits are not visible from 
the same ground station; e.g. Fairbanks; on every pass, which would complicate the 
Satellite control by requiring several control stations. 

A computer program was developed to study the nonstationary constellations and 
a description may be found in Appendix G. Basically, the program distributes the satel- 
lites about the world in a uniform manner, restricted only by the number of orbit planes 
and the total number of satellites. The orbital altitude and inclination are free variables 
but must be the same for all satellites. A region is specified and the program computes 
the number of satellites visible in that region. For large regions; e.g. , like those shown 
in Figure 4-1; the program generates a raster within the region; and at each point in the 
raster, computes the number of satellites visible. The program then determines the number 
of raster points where either three or four satellites were visible. Time is then advanced 
and new positions of the satellites relative to the earth are computed, and the sequence of 
determining the number of raster points where with three or four satellites were visible is 
repeated and at the completion of some time (24 hours) the program computes the total num- 
ber of raster points satisfying the visibility criterion and divides by the total number of 
raster points. This result is defined as the percent visibility. 

The results of this investigation are summarized by Figures 4-12 and 4-13. Figure 
4-12 is for 100% coverage, and Figure 4-13 is for 90% coverage. Comparison of Figures 
4-12 and 4-13 tend to show the soft degradation of the nonstationary constellations. For 
example, in Region 3, with a period of 8 hours, four satellite visibility requires 23 satel- 
lites for 100% coverage and 18 satellites for 90% coverage. Therefore, if 23 satellites 
were initially deployed and five failed, there would be long term outage of only 10%. 
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Figure 4-13. Number of Satellites for 90% Coverage 
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Notice the knee of the curves is around 6000 nmi which corresponds closely to the 
knee of the single satellite visibility curve (see Figure 4-14). 



Figure 4-14. Single Satellite Earth Coverage 

Figure 4-12 and 4-13 show that very attractive constellations exist with periods 
between 6 and 12 hours. Consider Tables 4-1 and 4-2, which were constructed from 
Figures 4-12 and 4-13. 


Table 4-1 

SUMMARY COMPARISON OF NUMBER OF SATELLITES REQUIRED VS 
ORBITAL PERIOD FOR REGION 2 AND FOUR- SATELLITE VISIBILITY 


Period (h) J 

— 
100% Coverage 

90% Coverage 

6 

26 

20 

8 

23 

17 

10 

21 

16 

12 

20 

15 

24 

18 

14 
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TABLE 4-2 


SUMMARY COMPARISON OF NUMBER OF SATELLITES REQUIRED VS 
ORBITAL ALTITUDE FOR REGION 2 AND THREE-SATELLITE VISIBILITY 


; Period (h) 

100% Coverage 

90% Coverage 

6 

22 

14 

8 

19 

12 

10 

18 

12 

12 

17 

11 

24 

14 

9 


The 10- hour period constellation with its 21 satellites could be deployed in three 
equally spaced orbit planes with seven equally spaced satellites in each plane. The 24- 
hour period with its 18 satellites could be deployed in three equally spaced orbit planes 
with six satellites in each plane. The difference between six and seven satellites per 
launch vehicle could be more than offset by the almost 2 to 1 reduction in orbital altitude; 
i.e. , the seven- satellite launch may require a smaller booster than the six-satellite launch. 
A failure of one satellite in each launch would mean 21 - 3 = 18 operating satellites for the 
10-hour period, and 18*- 3 = 15 operating satellites for the 24-hour period. Either would 
provide more than 90% coverage on the whole. 

Table 4-3 summarizes the effect on the number of required satellites vs region 1, 2, 
and 3. 

Table 4-3 

EFFECT ON THE NUMBER OF SATELLITES REQUIRED VS REGIONS 1, 2, OR 3 


Period (h) 

Region 

1 

2 

3 

6 

24 

26 

26 

8 

21 

23 

23 

10 

20 

21 

21 

12 

18 

20 

20 

24 

16 

18 

18 
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Percent of Time Satellites are Visible 


By comparing Table 4-3 and Figure 4-1, it may be concluded that the reason for 
the additional satellites is because Regions 2 and 3 extend into more southernly latitudes. 
In other words, expansion of Region 1 in a northward or east/west direction does not 
create the need for additional satellites, whereas expansion southward does. 

Figures 4-15 through 4-25 provide additional details. In particular Figure 4-16 
indicates that with orbital periods of less than 4 hours the increase in the required 
number of satellites is very rapid. For example, a 2-hour period requires about 100 
satellites for four satellite visibility which is about 4 to 5 times that required for 6- 
hour orbits. 
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Figure 4-15. Percent Coverage vs Number of Satellites, Region 1 
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Figure 4-17. Four-Satellite Visibility in Region 1 
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Figure 4-18. 


Three-Satellite Visibility in Region 1 
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Figure 4-20. Percent Coverage vs Number of Satellites, Region 2 
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It is possible that there are several combinations of the number of orbit planes, and 
satellites within a plane, which will yield the same coverage. Therefore, though the co- 
verage curves are plotted as if all total number of satellites were possible, the reader 
should recognize that some totals are impractical. Therefore, the next higher number, 
which provides a rational choice of orbit planes and satellites in a plane, should be se- 
lected. 

Certain total numbers have several possible combinations. For example, 24 satellites 
can be arranged into eight planes, three satellites in each plane, or three planes and eight 
satellites in each plane; or four planes and six satellites in each plane, or six planes and 
four satellites in each plane. The percent of coverage depends to some extent on which 
combination is used. However, reasonable selections; e.g. , three planes and eight satel- 
lites in each plane, or four planes and six satellites in each plane; provide about the same 
percentage coverage. Less reasonable selections, as 12 satellites in each plane and two 
planes, provide less coverage. 

Basically uniform coverage in an EW direction is obtained by a uniform spacing of 
orbit planes, whereas uniform coverage in a NS direction is achieved by a uniform spacing 
within an orbit plane. These comments are made for guidance and should not be taken too 
literally. However, in the equatorial regions, the reader can see the logic of the state- 
ment. The reader is cautioned that the curves generated herein are for trade-off purposes, 
and for optimization detailed investigations should be undertaken with regards to booster 
costs, in orbit maintenance, and booster failure rates. 

As mentioned previously, the lower altitude satellites have less free space loss which 
may have advantages relative to the accuracy of the position fix; i. e. , for a given EBP, 
the accuracy of the range or range difference improves as the range decreases. However, 
the free space loss gain is offset somewhat by the fact that as the altitude increases the 
satellite antenna gain increases because the beamwidth necessary for earth coverage 
(horizon to horizon) is less. 

The net loss is really the significant factor. Figures 4-28, 4-29, and 4-30 show the 
free space loss, the earth coverage antenna gain, and the net loss (free space loss minus 
antenna gain) respectively for three frequencies, 400 MHz, 1GHz, and 6GHz. Above 6000 
nmi, the net loss is about 1/7 dB per 1000 miles. Below 6000, nmi the net loss decreases 
rapidly with altitude. 
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Figure 4-27 . Earth Coverage Antenna Gain vs Altitude 
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DERIVATION OF RESOURCE PARAMETERS 
FOR MULTIPLE ACCESS TECHNIQUES 


The ensuing subsections briefly describe the prime performance derivations of 
resource parameters for multiple access techniques. 


A. 1 AUTONOMOUS TDM 


PG = 1, since all system users utilize the same code. Then the transmit time Tj 
available for each message is Tj = T u /N. However, since the transmissions are ran- 
domly generated, this transmit time must either be shortened or more time slots allo- 
cated to ensure an acceptable probability of a clear channel. The probability of a clear 
channel can be calculated as follows. Assuming the probability of a user occupying any 
interval, Tj, is uniformly distributed, then the probability of a user selecting a particu- 
lar slot is T{/T u , where T u is a 1-second frame interval to provide adequate data refresh. 


Then, probability of clear channel, P c , is 



for large N and for — «< 1 
u 


P = exp - (N - 1) 

C i 

u 


or 


p = 1 - ( N - l) — + .... 

c ' • 1 T 

u 


Now probability of error, P e , 
slot; is 

T. 

P e = (N - 1) ^ and P g = 
u 


caused by message overlap, e. g. , two users in same 


-5 5 

10 is assumed which yields one overlap per 10 messages 
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Then, 


, T 10 

T = — — 
i N 


-5 


The total bandwidth required can now be expressed as 


M MN 7 XT 

B = — = - — = 10 N 
T T. TP 
i u e 


where M = number of bits/message = 100 


The peak user transmitter power is expressed in terms of the desired output or de- 
tected signal-to-noise ratio. This is given as 


P T 
r 

^r 


D 


o 


T 

xp 


L_ G 
fs ac 


KT B. 

i 


G 

sc 


where 

Te> = one bit interval 

B. = detection bandwidth 
i 

= free space loss 

G = gain of aircraft antenna 
ac 

G gc = gain of spacecraft antenna 
KT = Johnson noise in 1-Hz bandwidth 
From this relation 
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This expression can be restated as 

36 + 10 log T - 10 log B. = 13 
xp a 

for the assumed system parameters and link power budget. Solving for Tjjp yields 

10 log T = -23 + 10 log B/ 
xp i 


or 


T = antilog 
xp 


-23 + 10 log — 
i 

10 


A. 2 ELICITED TDM 
For elicited TDM, 
T 

_ u 
T. = — 
i Rm 
T 


where 


R^, = the number of time slot assignments = N 


As with autonomous TDM, PG = 1; therefore, becomes 


_ M _ MN ,.2 XT _ 

b t “ tT c t ' t - ° 10 nc t 

1 U 


where 


C^, = the guard spaces between time slots and is set equal to one for purposes 
of system comparison 


Similarly, 


T - antilog 
xp 


' • , M 

-23 + 10 log — 

i 

10 



A. 3 AUTONOMOUS OR ELICITED FDM 

For autonomous or elicited FDM, PG = 1 , and = T u = 1. Then, each channel 
bandwidth equals B { = M/T if and the total bandwidth becomes 

B_ = N B. = N ^- = NMC_ 

T i T. F 


where 


= guard space between channels = 1 . 


and 


T = antilog 

xp 


-23 + 10 log = 7 - 
i 

10 


A. 4 AUTONOMOUS AND ELECITED CDM 


For autonomous and elicited CDM* 

1 ) The number of time slots, Rt = 1 

2) The number of frequency slots, Ry = 1 

3) The number of code users, Rq = N 

4) The transmit time per message, Tj = T u 

5) A signal-to-clutter ratio, SNR C , of 20 is assumed necessary to achieve 
P e =10- 5 . 

Then the processing gain PG may be expressed as 

PG = SNR N 
c 


and the bandwidth expression is 


B = PG ^ 
T T. 


___ NM 

SNR — 
- c T 

u 


N > 2 


An additional expression is introduced into the user transmitter power expression to 
ensure that the thermal noise is much less than the signal-to-clutter factor. A signal- 
to-thermal noise ratio, SNRt = 20 x 5 is assumed for CDM to ensure that thermal 
noise is negligible when referenced to the above clutter noise. 



The corrective factor of 5 shows up in the formula as follows: 


T = antilog 
xp 


M 

-23 + 10 log =- 
i 


10 ' +log 5 
A. 5 AUTONOMOUS TDM/FDM 


For A-TDM/FDM, the processing gain = 1, and the message transmit time is 


T. = 

1 r t 

where 

R = number of time slots per frequency slot 

However, to ensure a probability of an error of 10~ 5 for random overlapping messages, 
the expression must be modified as follows: 


T = 


T P 
u e 


i R„ 


Then, 


B = — • R^ 

T T F 
i 


xtt^ M ■ ' MN _ 

> R F ~ T P ‘ R T R F ~ T P N > R F 


u e 


u e 


or when N < R^, T^ = T u , i. e, , all users are accommodated via FDM; 
M 


B = ^ N 
T 


= — N 
N < R_ T 
F u 


N < R, 


Then, 


T = antilog 

xp 


-23+10 log^ 
i 

■ ~ 10 - 


= antilog 
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-23 + 10 log ^ 
i 
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or 


T = antilog 
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-23 + 10 log 
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N > R, 
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A. 6 AUTONOMOUS AND ELICITED FDM/CDM 

For A- FDM/CDM and E-FDM/CDM, 

T. = T 
1 u 

and the number of coded subscribers per FDM channel is 


Then, 

PG = SNR N = SNR ~ 
c s c R f 

and 


B = PG R = SNR N 
T T. F c 

1 


r - sm c 

u 


NM 

T 

u 


User transmitter power is identical to FDM when 2 > N < Rj and identical to CDM 
when N > R F . 


T =• antilog 
xp 


M 

-23 + 10 log — 
i 

lo ~ 


2 > N < R_, 
F 


T = antilog 
xp 


-23 + 10 log Y* 
i 

10 


+ log 5 


N >R f 


A. 7 AUTONOMOUS TDM/CDM 

The conditions for autonomous TDM/CDM are similar to elicited TDM/CDM ex- 
cept that the random nature of the technique causes overlap between time slots. This 
random overlap distribution has been assumed to be binomial, therefore, to ensure suf- 
ficient signal process gain for simultaneous users of time slots N g is redefined as 

N g = Three times the average numbers of simultaneous users per time slot. 


A -6 



In a binomial distribution, the average number users per slot is given as 


N = N — — 
s R T 

where 

= number of equally probable slots. 

Since the variance is also given as N/Rx letting N s - 3 N/Rx means that the probabil- 
ity of a greater number of simultaneous users is equal to that probability beyond the 3cr 
point on the binomial distribution. Then, 

N NT 

PG = SNR • 3 • -s- = SNR 3 -=-l 
R T C T „ 

and 


r, M 

B_ = PG — 
T T. 

l 


T.N , 

- „ l / 

, m \ 

„ NM 

= SNR 3 — 

t7 ) 

v 1 / 

- SNR 3 ?— 

N > 2 u ' 

N > 2 T u 


N > 2 


Similar to elicited TDM/CDM 


T = antilog 
• xp 


-23 + 10 log 


M 


10 


+ log 5 


A. 8 ELICITED TDM/CDM 

For elicited TDM/CDM, the number of code users within a time slot may be ex- 
pressed as 

N =iL 

s R t 

where Rx the number of time slots and 


T = 
i R, 


u 


As before, the processing gain may be expressed as 


PG = SNR N = SNR ~ 
c s cR t 
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Then, the bandwidth becomes 



SNR 


N 


c R„ 



= SNR 
N > 2 ° 


NM 

T 

u 


N > 2 


Now, the expression for the user transmitter power is identical to the pure TDM case 
when the number of users N < Rt and identical to the CDM case when N- R T : 


T = antilog 
xp 


-23 + 10 log f- 
i 

10 


2 > N < R t 


T = antilog 
xp 


M 

-23 + 10 log — 
i 

io — 


+ log 5 


N > R t 
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Appendix B 

DETAILED DATA SHEETS FOR 
MULTIPLE ACCESS TECHNIQUES 




Parallel Receiver Signal Processor 



CM 

I 




■8 

H 
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AUTONOMOUS TDM - EACH USER ASYNCHRONOUSLY 
BURST TRANSMITS HIS POSITION MESSAGE 

Concept: Autonomous TDM — Each User Asynchronously Burst Transmits His Position Message- 
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ELICITED TDM - USERS BURST TRANSMIT WITHIN A 



Parallel Receiver Signal Processor 




ELICITED TDM/CDM - GROUPS OR ORTHOGONALLY CODED 
USERS RESPOND SIMULTANEOUSLY ON COMMAND 

Concept: Elicited TDM/CDM — Groups of Orthogonally Coded Users Respond Simultaneously on Command 



Parallel Receiver Signal Processor 




Concept: Autonomous TDM/FDM 



Parallel Receiver Signal Processor 
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AUTONOMOUS OR ELICITED FDM/CDM 
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Concept: Autonomous TDM/CDM 
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PQ 
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Appendix C 


GRAPHICAL PRESENTATION OF NORMALIZED BANDWIDTH 
VERSUS NUMBER OF USERS 



Normalized Bandwidth 



1 10 10 2 10 3 10 4 10 5 


Figure C-l. Number of Users vs Normalized Bandwidth for 
Autonomous and Elicited FDM, TDM, and CDM 
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Number of Users (N) 

Figure C-4. Number of Users vs Normalized Bandwidth for 
Elicited TDM /CD M 


C- 


Normalized Bandwidth 



Appendix D 


MISCELLANEOUS POWER AND CODE GRAPHS 



Peak Power <W) 



Figure D-l. Number of Users vs User 
Transmitter Power for Competitive Techniques 
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Code Processing Gain (PG-p) 




Appendix E 

DETAILED DATA SHEETS FOR 
MULTIPLE ACCESS TECHNIQUES, C n 



AUTONOMOUS TDM CONFIGURATION H 




AUTONOMOUS FDM CONFIGURATION II 







AUTONOMOUS CDM CONFIGURATION H 















AUTONOMOUS TDM/FDM CONFIGURATION H 



Users/MHz 










AUTONOMOUS FDM/CDM CONFIGURATION H 
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AUTONOMOUS TDM/CDM CONFIGURATION II 
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Users/MHz 




















Appendix F 

DETAILED DATA SHEETS FOR 
MULTIPLE ACCESS TECHNIQUES, C IV 



E-TDMA, CONFIGURATION IV 


in 

o 
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Table F-5 
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Appendix G 


PROGRAMMING CONSIDERATIONS 



Appendix G 


PROGRAMMING CONSIDERATIONS 
G. 1 ERROR ANALYSIS PROGRAM 


The program first computes the satellite position, velocity, and the RTN to ECFC 
transformation matrix. The satellite computations are done in an inertial frame 
(assummed to be coincident with the ECFG frame at t = 0). 

The input parameters (for each satellite) are 

0^ Inertial longitude of ascending node 

i Inclination of orbit plane (0 < i <180°) 

£ Argument of perigee 

t Time of perigee 
P 

e Eccentricity 


T Period. 

The computations (for the jth satellite) are (reference Figure G-l) 

n. = 2 tt/T., mean rate 
J J 


a j = r /n i 


1/3 


, semimajor axis 


where 


(G-l) 

(G-2) 


M = 1.40765 x 10 16 ft 3 /s 2 


M. = n. (t - t mean anomaly 
J J PJ 


E. = M. = e. sin E., eccentric anomaly 
J J J J 


C . =2 tan 

J 


1R4 


true anomaly 


(G-3) 

(G-4) 

(G-5) 


G-l 







(G-6) 


R . = a. (1 - e.)/(l + e. cos £.), geocentric radius 
gj ) 3 3 3 


V. = 


.2 1 > 

M (r . ” a. 
g] J, 


11/2 


, speed 


fi. = 


cos 


e. sin i . 
3 3 


,. 1/2 


(1 + 2 e. cos X . + e. ) 
J J J 


, flight path angle 


<t> = 


- sin 


-1 


sin i; sin ( £ . - £ .) 
J 3 P3 


, latitude 


_1 /sin - C )c °s i\ 

£ . = 0 . + sin 7= — l| f inertial longitude 

J N] V cos 0 5 ^ 


where 


0. = £ - ^t, ECF longitude (Earth Centered Fixed) 

J 3 e 


-5 


oj = 7. 27221 x 10 rad/s 
e 


Define W. + W. so that 
3 3 


P ~ W p 

ECFC j ECRTj 


and 


p = W P + W P 

ECFC j ECHTNj j ECRTNj 


Then. 


w. = 

j 


cos (0 XT - <0 t) - sin (0- T - t) 0 
N. e N. e 

J J 

sin ^0 -a>t) cos -tot) 0 

N. e N. e 

J J 


0 cos i. - sin i 
3 J 


0 sin i. cos i. 

J 3j 


cos (£. - £ .) - sin ({. - £ .) 
) PJ 3 PJ 


sin (£. - £ .) cos ({. - £ ,) 
J b PJ J PJ 


and 


W = to 
j e 


sin (0 XT -tot) cos(0 XT -tot) 0 
N e N. e 

j 3 

- cos (0 -to t) sin ( 0 XT -tot) 0 
N. e N. e 

J J 


0 cos i. - sin i. 

J J 


0 sin i. cos i. 

3 3. 


cos({. - £ .) -sin({. - £ .) 
*J PJ J PJ 


sln(£.-£ .) cos(l - £ .) 
J PJ j '’PJ 
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(G-7) 

(G-8) 

(G-9) 

(G-10) 

(G-ll) 

(G-12) 

(G-13) 

0 

0 (G-14) 

1 

0 

0 (G-15) 



The position (P 


P . = 
SJ 


P . = 
S3 


s; 


r ~i 

sy 

X. 

J 


R . 
gJ 



Y 

3 

= w 

J 

a 



Z. 

3 


0 



X. 

3 


R . 
g] 


V. cos B. 
3 3 

Y. 

3 

= 

3 

, 0 

■+ w. 

3 

V. sin B. 
3 3 

A 


0 


0 


(G-16) 
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The Navigator's position, P^, and velocity, P^, in ECFC are computed as 


P N = 


and 


P N 5 


x„ 



cos L cos L 





a o 


y n 

" + 

h) 

cos L sin L 
a o 


Z^ T 



sin L 


L n j 



a 


X 


-sin n sin L - cos n sin L cos L 1 

N 



o 

a o 

Y__ 

= V^ T - 

sin n cos L - cos n sin L sin L 

N 

N 


o 

a o 

Z 


cos n cos L 


Nj 


L 

a 



(G-18) 


(G-19) 


where 

L , L is the ECFC latitude and longitude of the navigator 
a o 

R = 2.092567 x 10 7 ft 
e 

h, V XT , n are input quantities 
N 

h = altitude of the navigator 

= navigators speed relative to earth 
n = the ECFC heading angle from north. 

(Note: The navigators flight path angle is defined to be zero which is 
equivalent to flying horizontal) 
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The range, R., and the range rate, Rj, to the jth satellite are 

|l/2 

j = 1, 2, 3 


B J ' [^ - X / + - ¥ )» 2 + (Z N - ¥ 

S ~i. K-¥^-¥ + < V¥ ( V y j> 


+ ( Z N - z j) < z 


ij - Z.)1(G- 
N j _ 


(G-20) 
21 ) 


For the range difference navigation schemes, the jth range difference, Pj, and range 


rate difference, Pj, are 


P. = R - R, 
] j+1 1 


P. = R - R 
3 3+1 1 


The error equations (by navigation type) are 






AP 

SI 

APf' 




AP 

S2 

AP 2 




AP 

S3 

AP 3 

sf- 

ii 

" ap n 

+ n t 

AP 

S4 

AP 1 

ap n 

AP 

SI 

AP 2 

1 


aP 

S2 

AP 3. 




AP 

S3 





AP S4 
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(G-23) 


(G-24 3 p) 


' Ap i' 




" ap si‘ 

AP 2 




AP S2 , 

Ah 

“ m t 

AP 

N 

+ n t 

AP 

S3 



AP 1 


AP 

N 


AP 

SI 

AP 2 




AP 

S2 

< 




aP 

S3_ 
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v 2P+H ; 
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AR 1 




AP S1 

AR 2 




A? S2 

AR 3 

' ll 

AP 

N 

+ n t 

AP 

S3 

ARi 

AP 

N 

AP 

SI 





AP S2 

_ A \ 




_ AP S3_ 


AR 2 




AP S1 

Ah 


ap n 

+ 

T 

AP S2 

A R 1 

_ ap n_ 

ap si 

AR 2 




AP S2 


Ah 


(G-24 3R ) 


< G - 24 2 R + H> 


A condensed notation for the error equation for ail navigation schemes is 


AP 

AP 


“t 


AP 


N 


AP. 


N 


+ N„ 


AP, 


AP, 


(G 24 general ) 


The matrix, M T , is 




M 


M 


O 


M 


smce 


d R. dR. 


, etc. 


(G-25) 
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Where 



G-7 



or 



(cont) 


G-8 




9R 9R 
• 3 3 

when evaluating R , R , — — , — ■— 
3 3 3 3 

The matrix is 



where 



(G-27) 



G-9 







The error Equation 0-24^^, can be solved for navigation errors: 


. -1 AP -1 

■ M T ~Z - m t h i 

__AP_ 


where 




-1 

r -l 

M 

! o" 
1 _ 


M 

, 



-1 . 

M 

M_ 


-M MM 


0 [“A J 0 

"IT ” ” i 

1 L A ' A J 


A more relevant coordinate system for the navigation errors is an earth surfaced 
tangential (ESTC) frame located at the navigator. Define Q such that 


Q AP X 


ECFC 


East component of AP 


North component of AP^. (at the navigator) 


Up component of AP 


ESTC 


Then, 


- sin L 


cos L 


Q = - cos L sin L - sin L sin L cos L 

o a o a a 


cos L cos L sin L cos L sin L 


cos L 


- sin L 


Q = L sin L sin L - cos L sin L 0 


o a 


o a 


- sin L 

o 


cos L cos L cos L 0 

a o a 


+ L - cos L cos L 


cos L sin L 
o a 


sin L cos L 
o a 


sin L 


sin L sin L cos L 

o a 2 
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(G-35) 


(G-36) 


(G-37) 



where the subscript g indicates ESTC frame, and the covariance matrix of navigator 
errors is 



G-13 



and 



(G-41) 


(G-42) 


(G-43) 


for the 3P navigation mode, is suitably reduced for navigation schemes using two 
or three satellites. 


If QA and QANW are defined as 



G-14 



QANW = QANW + QANW + QANW + QANW 


(G-45) 


then the equation producing the covariance matrix of the navigator errors, can be written 
as implemented : 


NAV 


QAt 2 (QA) T i QAt 2 (QA) T 

_pp I PP 

• ? T T • ? • T 

QArl(QA) : QAt (QA) 
PP | PP 


r i 
0 ! 


o 


0 | Q A t 2 (QA) T 
1 PP 

I 


QANWr 2 (QANW) T ! QANWt 2 (QANW) T 

. : SS l ' SS ^ 


T ! 


QANWt (QANW) | QANWt (QANW) 

SS i SS 


(G-46) 


0 ! QANWt (QANW) 

i SS 


For the case where a reference station calibration is assumed, the navigator error 
Equation 38 becomes 


AP M,rl 


AP 


AP 

Ng 

= Q A 

rp 

— 

i 

<o 

> 

1 

5 

s 

1 

i 1 

AP 

i X X 

1 


- Q A 
T T 
o o 


AP 


AP 
L oj 


+ Q t a t n t 

o o o 


AP 


AP 


oJ 


(G-47) 


where the subscript, o, indicates that the quantity is evaluated at the time of the calibra- 
tion using the satellite/navigator geometry defined by placing the navigator at the reference 
station at the time of calibration. 

Then the covariance matrix of the navigator errors (with reference station calibration) 
becomes 


REF 


T NAV + \ T NAV ) ~ T CROSS 




2 ) 

CROSS ) 


(G-48) 
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where 


CROSS 


< QA >o T p p (QA)Tl 
o 


K t p o p< qa » T 


« A >o T p o p «' A > T 


_i_ 2 — T 

(QA)oTp p (QA) 
o 


0 j 0 

i * 

0 ! «Vp p ,QA)T 

U i o 


(QANW) T 2 (QANW) T 
o s s 
o 


(QANW) r (QANW) 
o 


(QANW) t (QANW) 
0 o 


(QANW) t (QANW) 
o 


'(G-49) 


i 


0 ' (QANW) t 2 QANW) T 
I' oss 

L I o 


2 2 2 *2 

The terms r , f , T , f (defined in Subsection G. 2) relate the source 
ppppssss 
o o o o 

errors at calibration to the current source errors. 

For example 
2 


T = E 

»o P 


{ap o ap t } 


(G-50) 


The program computes the diagonal elements of 


2 , 2 
t nav 31111 T REF 


for each of the error sources defined in Subsection G. 2 and prepares the following 
quantities for output: , : 


< 


2 2 

r... tr + T XTATr = Horizontal Position Error 

NAV 1,1 NAV 2,2 




NAV 


= Vertical Position Error 


3,3 


i 

I' 


2 2 

T NAV + T NAV = Horizonta ^ Rate Error 
4,4 5,5 


NAV. 


= Vertical Rate Error 


6,6 


S (G-51) 


Similar expressions using are formed and output for the four components of error 

with reference station calibration. 
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G. 2 ERROR SOURCES 


G.2.1 3 P NAVIGATION 

AR^ARj 

AR - AR 
4 1 

AR„ - AR, 
2 1 

AR - AR, 

3 1 

AR - AR, 

4 1 


AP 

AP 


The i j element of T is equal to 
PP 

E (AP t AP.) = E (AR. + t AR. + 1 ) + E (AR 1 AR^ 
-E(AR. + 1 AR 1 )-E(AR 1 AR j+i) 

2 

the i i element of T is 

V 

E (AP. AP.) = E (AR. + 1 AR. + 1 )+E(AR 1 AR^) 

O O' 1 o 

-E(AR. + 1 AR i )-E(AR l AR ) 
o o 3 

• 2 

the i j element of T is 
PP 

E (AP. AP.) = E (AR. AR. ) + E (AR AR ) 
l 3 1 J * 1 11 

-E(AR. + 1 AR 1 )-E(AR 1 AR. + 1 ) 


(G-52) 


(G-53) 


(G-54) 


(G-55) 
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2 

and the i i element of f is 

P o P 


E (AP. AP.) = E (AR. J , AR. ,) + E (AR ARJ 
1 j' ' l+l j+r ' 11 


-E(AR i+1 AR 1 )-E(AR i AR. + 1 > 


(G-56) 


G. 2. 1. 1 3P Navigation, Range (and Range Rate) Noise Errors 

When normalized to a unit error, the expected values needed to evaluate Equations G-53 
through G-56 are 


E(AR i+1 AR j+1 ) = E(AR i+1 AR j + 1 )=8 ij 


(G-57) 


E (AR X ARj) = E (AR 1 AR J ) = 1 


(G-58) 


where 8 . . is the Kronecker delta function. 

The remainder of the terms have the value zero. Therefore, 
2 1 1 


2 . 2 
T “ T 
PP PP 


12 1 
112 


(G-59) 


and 


2 • 2 r A 1 

T - T = [ 0] 

P 0 P p o P 


(G-60) 


G. 2. 1. 2 3P Navigation, Range (and Range Rate) Bias Errors 

When normalized to a unit, error, the expected values needed to evaluate Equations G-53 
through G-56 are all equal to the value one. 


Therefore, 


2 .2 2 -2 r „ , 

T = T =T = T = L 0 ] 
PP PP P P P P 

o o 


(G-61) 
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G. 2. 1.3 3P Navigation, Refraction Errors 

When normalized to a unty error, 

E (AR. AR.) = p . 
i 1 i] 

and 

E (AR. AR.) = p. . 


for all where p . . and p . . are defined in Subsection G. 3. 

13 V 

Therefore, 




(G-62) 


(G-63) 
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G. 2.1.4 3P Navigation, Satellite Position (and Rate) Errors 

The ephemeris errors are normalized to the tangential axis. The radial and normal 
errors are related to the tangential error by the input factors kj and k 2 . The ephemeris 
errors are assumed to be correlated over the four satellites by the input factor p s . The 
source error covariance matrix is 


2 . 2 
T = T = 
ss ss 


K 

’s K 


p s K 

P S K 

K 

"s' 5 
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p s K 

:p s k " 
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P s K 

"D 
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° 8 k 

K 
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where 


K = 


“1 

0 


0 

1 

0 


0 

k 2 

s 


(G-67) 


2.2 

Two values of r and t are used. These represent the limiting cases of 
oo 

temporal correlation, as follows: 


2 2 

r = t = [ 0 ] (uncorrelated) 
s s s s 

O O V 


(G-68) 


and 


t = r = t (correlated) 
s s s s ss 
oo 


(G-69J 
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G.2.2 2P+H NAVIGATION 


AR - AR, 
2 ' 1 

AR - AR, 
3 1 


AR - AR, 
2 1 

AR - AR, 
3 1 


(G-70) 


2 .2 2 *2 

For i and j both less than 3 , the i j elements of t , t , t > and t are given 

PP PP P Q P P 0 P 

by Equations G-53 through G-56, respectively. 

For i or j equal to 3, the i j elements are 


T (ij) - E (AP. AP.) = 8. . E (Ah Ah) 
PP i J i] 


(G-71) 


T 2 (i j) = E (AP. AP.) = 8. . E (Ah Ah) 
PP i ] iJ 


(G-72) 


T (i j) = E (AP. AP.) = 8. . E (Ah Ah) 
P 0 p 1 Q 3 1J o 


(G-73) 


T (ii) = E (AP. AP.) = 8 . E (Ah Ah) 
P 0 P » 0 3 iJ o 


(G-74) 


where 8. . is the Kroenecker delta function. 
iJ 

G. 2. 2. 1 2P+H Navigation, Range (and Range Rate) Noise Errors 
Using Equations G-57 and G-58 produces 


2 1 0 

T 2 =T 2 = 1 2 0 

PP PP 

0 0 0 
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and 


2 - • 2 r a! 

T “ T = [ 0 ] 

P o P V 


(G-76) 


G. 2. 2. 2 2P+H Navigation, Range (and Range Rate) Bias Errors 
Similarly to the 3 P case (G. 2. 1. 2), 

2 - 22 . 2 - - 

T = T =T=T=lO] (G-77) 

PP pp p p p p 
o o 

G.2.2.3 2P+H Navigation, Altitude (and Altitude Rate) Noise Errors 

Normalizing to a unit error and using Equations G-71 through G-74, gives 

*0 0 (f 


2 . 2 
T = T 
PP PP 


0 0 
0 0 


(G-78) 


and 


T 2 — T 2 = [0] 

P o P P o P 


(G-79) 


G. 2. 2.4 2 p+ H Navigation, Altitude (and Altitude Rate) Bias Errors 
Normalized to a unit error, 
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0 1 
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G. 2. 2. 5 2P+H Navigation, Refraction Errors 
Similarly to the 3P case (G. 2. 1.3), 
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.2 

and t are not applicable for this error source. 
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the i i element of t is 
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G. 2.3.1 3R Navigation, Range (and Range Rate) Noise Errors 

When normalized to a unit error, 
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G. 2. 3.2 3R Navigation, Range (and Range Rate) Bias Errors 

When normalized to a unit error 

E (AR. AR.) = E (AR. AR.) = E (AR. AR.) = E (AR. AR.) = 1 
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G. 2. 3.3 3R Navigation, Refraction Errors 
Equations G-62 and G-63 hold; therefore, 
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G. 2. 3. 4 3R Navigation, Satellite Position (and Rate) Errors 


2 2 2 • 2 
T and t are givenby Equation G-83. The two values of t and r are given 

OC3 dCJ S S S S 

o o 


ss ss 
by Equations G-68 and G-69. 


(G-97) 


G. 2. 4 2R+ H NAVIGATION 
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For i and j both less than 3, the i j elements of r > T , T , T are given by 
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Equations G-85 through G-88, respectively. 

For i or j equal to 3, the i j elements are given by Equations G-71 through G-74. 

G. 2.4.1 2R+H Navigation, Range (and Range Rate) Noise Errors 
Using Equations G~89 and G-90 produces 
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G. 2.4. 2 2R+H Navigation, Range (and Range Rate) Bias Errors 
Equation G-93 holds, and 
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r = t = r = r 
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G. 2.4.3 2R+H Navigation, Altitude (and Altitude Rate) Noise Errors 
2 . 22*2 

r , r , T , f are given by Equations G-78 and G-79. 
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G. 2.4.4 2R+ H Navigation, Altitude (and Altitude Rate) Bias Errors 

2 . 22.2 

T , t , T , f are given by Equation G-80 
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G. 2.4. 5 2R+H Navigation, Refraction Errors 
Using Equations G-62 and G-63. 
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G. 3 REFRACTION CORRELATION FACTORS 
The term 



represents the correlation factor between the range measurement errors due to refrac- 
tion effects along the two paths 


R. andR., 

1 3 

the vectors between the navigator and the satellites i and j at the current time. From 
Chen T s paper, "Range Difference Error due to the Presence of Ionospheric Electrons", 


O. . = G. G. (1 - 0. 665 x 10 -7 S. .) 
i] i J ij 
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where S. . is the distance, at ionosphere altitude, between R. and R. in feet, and 
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with 


L , the latitude of the navigator in radians 
a 

E^, the elevation angle of R^ in radians. 
The term, 



represents the correlation factor between the measurement errors due to refraction 
effects along the two paths 
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and R., 
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where 


R. is the vector between the reference station and satellite i at t = t . 

1 o 

o 

R . is the vector between the navigator and satellite j at the current time. 
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Then, 
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with 

L , the latitude of the reference station in radians 
a 

o 

E. , the elevation angle of R. in radians, 
o o 
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MULTIPLE SATELLITE COVERAGE 


This program is designed to provide coverage patterns for a specific number of 
satellites in a specified number of orbit planes. Distribution of the satellites takes 
place internal to the program and is based on the following input data: 

1) Number of satellites - Maximum of 153 

2) Number of orbits - Must be an integer factor of the above 

3) Satellite period - T (minutes) 

4) Eccentricity - e 

5) Time of perigee for first satellite - T p (1) (minutes) 

6) Longitude of ascending node for first satellite - (1) (degrees) 

7) Angle of inclination - i (degrees) 

8) Argument of perigee for first satellite - Cn (1) (degrees) 

The intent is to provide uniform spacing of the satellites about the globe, initially 
and as the satellites move within their orbits. 

. The secular effects of oblateness, i. e. , regression of the line of nodes and argu- 
ment of perigee, are accounted for in the satellite motions. The following criteria 
have been used to properly space the satellites. 

The space between orbital planes is equal to 180° /number of orbits. Start the 
first orbital plane off at the equator and Greenwich. Space the satellites uniformly 
within the orbit by changing the time of perigee (T p ). Delta T p is equal to period/ 
number of satellites per orbit. 

Example 1 (Figure H-l) 

8 satellites t = 90° 

2 orbital planes 
4 satellites per orbit 
2 hour satellite period 
Delta T p = -T/4 = -30 minutes 
Delta £ = 180°/2 = 90" 
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Example 2 (Figure H-2) 

15 satellites 

l 

3 orbital planes 
5 satellites per plane 
2 hour period 
i = 90° 

Delta Tp = 120/5 = 24 minutes 
Delta £ = 180°/3 = 60° 

Thus, the second and third orbital planes will be spaced 60° and 120° away from the 
first. 

Delta £ p = 60% = 20° 

Thus, the argument of perigee for the second and third orbital planes will be 20° and 
40°, respectively. 

Plane 1 

£ n = 0° 

•Cp/o.o°‘. 

T p = 0.0, -24.0, -48.0, -72.0, -96.0 

Plane 2 

£ N = 60.0° 

Cp = 20.0° 

T p = 0.0, -24.0, -48.0, -72.0, -96.0 

Plane 3 

4 N = 120.0° 

C p = 40. 0° 

T p = 0.0, -24.0, -48.0, -72.0, -96.0 

Once the satellites have been spaced properly, it is possible, at any given time, 
to display a graphic representation of where all the satellites are located and how 
many can be seen from any point on the earth. At time t, the vector for each satellite 
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Figure H-2 

is computed and stored. At every 5° latitude and longitude, a test is made to deter- 
- mine how many satellites are visible, i. e. , above 10° elevation angle to meet naviga- 
tional requirements. 


_ Altitude (R s ), latitude (4> B ), and longitude ( 6 S ) of a satellite produces the vector 
R g (see Figure H-3). See Program flow for equations deriving R g , «#> s , and 8 S . 

cos <f > g cos 0 g ' 


Hg Rs 


cos 4> b sin 0 S 
sin <*> s 


Altitude (Rp = radius of earth), latitude (<£p), and longitude (0^) of a position on 
the earth's surface gives the vector Rp (see Figure H-3). 


R p “ Rp 


cos <t > p 

cos 6^ 

COS </> p 

sin ^ 

^sin <p p 
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Figure H-3 

A line (Rel) drawn from a satellite tangent to a point on the earth gives the maxi- 
mum distance that satellite can be from the point on the earth (see Figure H-4). Thus, 

o O 

Eel = Rg - Eg where, R e = radius of earth 
The vector R^ (point on the earth to satellite) is computed by: 

^ps ” R s “ ^p 
(See Figure H-4. ) 


If Ep S > Rel, the satellite is not acceptable since it is beyond the horizon. 
If not, the elevation angle of the satellite ( € ) is computed as follows : 


R • 

R 

ps 

P 

|r I 

R 

1 Ps| 

P 


e = 90° - J3 

- H the elevation angle is above 10°, the satellite is therefore acceptable for naviga- 
tional purposes. 
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Program Flow 



ACELV 

STIM 

ETIM 

DTIM 

DLAT 

DLONG 

NSTOT 
NORB 
PERIOD = T 
!e 



fp 
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SYMBOLS 


R - Radius of earth 
e 

R - Distance from center of earth to observation point 
P 

■ 'to - Earth's rotational rate 
e 

a - Semimajor axis 

e - Eccentricity 

T - Time of perigee 
P 

4^ - Longitude of ascending node 

i - Inclination 

4 - Argument of perigee 

P 

T - Satellite period 

n - Mean motion 

p - Semilatus rectum 

n - Mean value of mean motion 

R - Distance from center of earth to satellite 
s 

R - Distance from satellite to observation point 
ps 

in - Mean value of mean anomaly 

- Mean value of perturbed true anomaly 

4^ - Mean value of longitude of ascending node 

<f> - Geodetic latitude of satellite 

s 

4 - Inertial longitude 

0 - Longitude of satellite 

s 

d> - Latitude of observation point 
P 

0 - Longitude of observation point 

P 
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INPUT 


Card 1 

ACELV 

i 

F10.5 

Acceptance value for angle of elevation (degrees) 


STIM 

F10.5 

Start time for observations (minutes) 


ETIM 

F10.5 

End time for observations (minutes) 


DTIM 

F10. 5 

Time increment (minutes) 


DLAT 

no 

Delta latitude (degrees) 


DLONG 

no 

Delta longitude (degrees) 

Card 2 

NS TOT 

15 

Number of satellites 


NORB 

15 

Number of orbits 

T 

PER 

F10.5 

Period of satellites (minutes) 

e 

E 

F10.5 

Eccentricity 

T 

P 

TP(1) 

F10.5 

Time of perigee for first satellite (minutes) 


XIN(l) 

F10.5 

Longitude of ascending node for first satellite 
(degrees) 

i 

AIN 

F10.5 

Angle of inclination (degree) 


ARGPER 

F10.5 

Argument of perigee for first satellite (degrees) 
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Appendix I 

MARITIME USE OF SATELLITES 

LI INTRODUCTION 


During the past 25 years or more, electronic aids to navigation and collision avoid- 
ance have been used aboard ship. While certainly improving matters, each of the var- 
ious systems has its problems. For example, Loran A has a range limitation of 600 
to 900 miles by day and 1400 miles by night, insufficient for a truly world-wide system. 
Loran C has a longer range, but its lower reliability and higher cost make it unsuitable 
for commercial users. 1 Radar has proven so unsatisfactory with respect to collision- 
avoidance that the term "radar-assisted collisions" has come into use. This is not to 
say that the equipment causes collisions, but that ship’s personnel either have difficulty 
interpreting the information provided by the radar or do not take the trouble to deter- 
mine the proper action. The interpretation problem can arise from interference, such 
as rain and sea clutter, traffic density, which stresses the ability of even a competent 
operator, or simply lack of training. Recognition of the density problem has led to the 
development of a variety of automated systems capable of tracking multiple targets. 2, 
These systems certainly ease the burden of the individual ship’s officers. However, it 
remains to be seen whether the systems prove accurate enough for all types of users and 
suitable for ever-increasing traffic densities. The latter in particular will highlight, if 
it hasn’t already, the absence of a marine equivalent of the current air traffic control 
arrangements among nations. Satellites appear an attractive approach both to increased 
accuracy and cooperative traffic control. 


1.2 USER COMMUNITY AND ITS REQUIREMENTS 


Table 1-1 summarizes the results of a study performed for the US Coast Guard in 
1969. 5 This study provided estimates of current quantities, and forecasts for 1980, of 
the world’s maritime vessels exceeding 100 gross tons. Not surprisingly, those vessels 
engaged in the transportation of raw materials and the search and acquiring of food had 
the fullest utilization. In all, there were in 1969 some 21,000 vessels at sea, of which 
nearly half were fishing vessels. The projections for 1980 were for an almost 50% in- 
crease in ships at sea, most of the increase accounted for by yet more fishing vessels. 

A geographic distribution appears in Table 1-2. 


While the overall traffic does indicate the magnitude of the navigation/collision- 
avoidance problem, the predominance of fishing vessels in itself has certain consequences. 
The traffic patterns of such vessels are much less predictable than for carriers of raw 
materials or manufactured products, as fishing boats are engaged in continuing search 
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Table 1-1 


FORECAST OF VESSEL POPULATION GROWTH 


Type of Ship 

1969 

1980 

Total 

At-Sea 

Total 

At-Sea 

Tankers 

5,869 

4,637 

7,039 

5,559 

Ore and bulk carriers 

2,378 

1,474 

5,370 

3,519 

General cargo and passenger 

26,100 

5,269 

19,565 

4,935 

Other 

3,980 

270 

2,984 ! 

254 

Total exclusive of fishing vessels 

38,327 

11,650 

34,958 

14,267 

Fishing vessels 

11,949 

9,323 

19,156 

14,946 

Overall total (over 100 gt) 

50,276 

20,973 

54,114 

29,213 


Table 1-2 

SUMMARY OF SHIP DISTRIBUTION TRENDS 


Location 

1969 

1980 

Atlantic 

9,507 

11,130 

Pacific 

7,431 

10,790 

Indian 

2,465 

5,776 

Miscellaneous 

1,570 

1,517 

Total (over 100 gt) 

20,973 

29,213 


-2 







operations. Moreover, such vessels ever- increasingly operate in fleets, so that the 
collision- avoidance problem is aggravated for ships encountering them on the high seas 
and perhaps within the fleets themselves. Still further, fishing boats require an accu- 
racy of less than 0. 1 mile in all oceans and under virtually all weather conditions, and 
for reasons of business security would prefer a navigation system which does not require 
position reporting. 6 

Table 1-3 presents accuracy and other user requirements for a variety of ship types. 7 
Considering the usual estimates of satellite lifetimes, (5 to 7 years), the 1975-1985 col- 
umn is of the most interest here. Note that the 0. 1-mile or better accuracy is a require- 
ment for several types of vessels, with 0.5-mile or better specified for all but recrea- 
tional craft. The requirement for "instantaneous" determination of position mandates 
some form of automated system. Coverage capability of-world-wide or near-world-wide 
extent underscores the need for a multination approach to traffic control, weather report- 
ing, and search and rescue missions. 


1.3 SATELLITE TECHNIQUES AND THEIR ADVANTAGES 

A variety of experiments 6 have demonstrated that satellites can provide satisfactory 
marine surveillance and communications under a wide variety of conditions, some of 
which render HF and lower frequencies unusable. In fact, satellite communications are 
about as good as bridge-to-bridge communications in a harbor, which themselves are 
very useful. Obviously, the flexibility of such a system would be far greater than a line- 
of-sight system using such media as HF radio and cable. One of the questions to be 
answered (but not here) is that of orbit altitude. The Navy ! s TRANSIT system, which is 
the only fully operational satellite navigation system for marine use, has satellites in 
600-mile polar orbits. One analysis 6 alleges that this type of system is less costly 
than one using synchronous orbits, and states that a synchronous system is superior only 
when service over a limited portion of the globe is required. However, the cost figures 
in that analysis are quite close, and one can argue that limited coverage per satellite or 
constellation is desirable as it minimizes tracking and handover problems. Another 
analysis 10 indicates a preference for mid-altitude satellites (6,500 to 14,000 miles). 

Still another study 11 refers to "preferred deployments" of geostationary orbits. How 
a marine satellite system eventually would be configured has yet to be resolved, but the 
TRANSIT system has already demonstrated an accuracy capability meeting the require- 
ments set forth in Table 1-3. 12 

A cost-benefit analysis for six possible satellite configurations was performed by the 
National Academy of Sciences. 13 These configurations were 

1) Not appropriate for marine application 

2) Two geostationary satellites using the VHF band and having limited 
capacity for marine service, primarily traffic control 


1-3 



Table 1-3 

MAR TIME USER REQUIREMENTS 



1-4 


/ 








































3) Two geostationary satellites using both the VHF and L bands, and pro- 
viding position fixing and communications for both air and marine traffic 
surveillance 

4) A similar system using only L band 

5) One geostationary satellite with communications in L band and position 
fixing at L or C band 

6) A worldwide system (Systems 1-5 are transoceanic) using six geo- 
stationary satellites. 

The costs of these systems are summarized in Table 1-4. The R&D and the launch costs 
are assumptions, as is a satellite lifetime of 5 years and a launch malfunction factor of 
1.2. The nature and amount of the benefits is summarized in Table 1-5. It turns out 
that there is an unquestioned advantage (defined by the intersection of the lowest estimate 
of benefits with the highest estimates of cost) for Systems 3, 4, 5, and 6 if the systems 
are available to both US and foreigh users. The points of advantage are reached from 2 
to 4 years following the date of availability of the satellite if aircraft as well as maritime 
applications are included. System 2 takes longer than ten years to reach the same point. 
Obviously, other configurations would have different cost-benefit tradeoffs, but the 
NAS study illustrates the kind of analysis that can be made to justify a continuing effort. 
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Table 1-4 

CQSTS ASSOCIATED WITH SATELLITE SYSTEMS 


Capital Equipment Costs ($ Million) 

System 

1 

2 

3 

H 

. 5 

6 

R&D 

10-20 

10-20 

20-40 

20-40 

30-50 

20-40 

No. of Satellites 

1 

2 

2 

2 

1 

6 

Cost of Flight Model Satellites 

2-4 

2-4 

6-9 

6-9 

7-10 

6-9 

Cost of Launch (per satellite) 

5 

5 

15 

15 

15 

15 

Total Cost of Satellite Portion of 
System (R&D plus Satellites plus 
Launch, including 1.2 Malfunc- 
tion Factor) 

18-31 

27-42 

70-98 

70-98 

56-80 

171-213 

Ground Complex Costs (including 
Marine Control for systems 3-6) 

3-4 

3-4 

6-8 

6-8 

6-8 

. 18-24 

Total Capital Costs 

21-35 

30-46 

76-106 

76-106 

62-88 

189-237 

Operation and Maintenance of 
Ground Complex (per year) 

2-3 

2-3 

3-5 

3-5 

3-5 

9-15 

Replacement of Satellites (per 
year) 

1.6-2 . 

3.2-4 

5-11.5 

5-11.5 

5.3-6 

15-35 

Total Annual Costs 

3.6-5 

m 

8-16.5 

8-16.5 

8.3-11 

24-50 

User Equipment Costs ($ Thousand) 

Capital Costs (Equipment + Instal- 
lation) 







Aircraft Users (each) 

20-40 

20-40 

20-40 

60-120 

60-120 

60-120 

Ship Users (each) 

* 

7-15 

25-50 

25-50 

25-50 

25-50 

Operation and Maintenance (per 
year) 







Aircraft Users (each) 

2-4 

2-4 

2-4 

6-12 

6-12 

6-12 

Ship Users (each) 

* 

0.7-1. 5 

2.5-5 




*Not appropriate 
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Table 1-5 


BENEFITS OF VARIOUS SATELLITE SYSTEMS 




/ 

Average Monetary Return 
per Year per Aircraft or Ship 

Benefit 

Bene- 

ficiary 

Reason for Benefit 

Average 

Net (10 year 
Amortization of 
Costs of User 





Equipment 

Improvement 

Owners 

Savings in fuel and 

Reduction in delay 

$2-4k rising to 

in Air 

of 

increased utiliza- 

for 700 flights per 

$44-8 8k 

Traffic 

Aircraft 

tion of aircraft 

year at $20-40 


Control 


by reduction in 

savings per flight 




en route delays 

in an early time 
frame, rising to 
$80-160 per flight 
by 1980, $14k- 
28k to $56-102k 


Marine 

Owners 

Reduction in insur- 

Lowering of negoti- 


Traffic 

of 

ance premiums 

ated rate of 


Control 

Ships 

through increased 

premium if ex- 




safety 

perience warrants 
it by 0. 5%-l. 0% 
of cost of vessel, 

(A ship would re- 




($10M-20M aver- 

ceive benefits 




age cost) $50k- 

from both with a 




200k 

single installa- 

En route 
Navigation 

Owners 

of 

Savings in fuel and 
increased utiliza- 

1% of fuel cost 
(300k-500k); 1% 

tion of equipment) 
$51k-$202k 

for Ships 

Ships 

tion of ships by 

gain in days avail- 




improved course 

able for revenue 




keeping 

carrying purposes 
($300k-$700k per 
year), $6k-$12k 


Aid to 

Owners 

Savings in man- 

10% of manpower 

Zero-$15k 

Commercial 

of 

power and fuel by 

and fuel costs, 


Fishing 

Ships 

enabling ship to 

($50k-$250k per 




return directly to 
proven fishing 
grounds 

year) $5k-$25k 



* Column 4 less user costs from Table 4 


k = Thousand 
M= Million 
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